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abstract 


Several  materials  undergo  transitions  from  a  semiconducting 
to  a  metallic  state  at  a  critical  temperature.  Previous  theoretical 
attempts  to  understand  such  transitions  have  been  generally  quali¬ 
tative  and  have  not  been  able  to  account  for  all  the  specific  experi¬ 
mental  results. 


In  this  work,  an  explanation  of  semiconductor-to-metai  transi¬ 
tions  is  presented  using  a  band  model.  It  is  shown  thermodynamically 
that  the  energy  gap  of  a  semiconductor  closes  down  significantly  with 
the  number  of  excited  carriers  if  the  gap  has  a  large  pressure  coef¬ 
ficient,  as  is  found  in  several  of  these  materials.  This  shrinkage  of 
the  energy  gap  is  due  to  explicit  variation  of  the  crystalline  volume. 
There  may  also  be  a  constant  volume  carrier  concentration  dependence 


of  the  gap,  which  cannot  be  evaluated  thermodynamically. 

Two  specific  models  are  discussed.  If  the  energy  gap  arises 
from  the  splitting  of  the  first  Briliouln  zone  by  an  antiferi  omagnetic 
exchange  interaction,  the  gap  will  decrease  linearly  with  the  number 
of  free  carriers.  The  relationship  is  demonstrated  by  means  of  a 
one-dimensional  model,  a  three-dimensional  tight  binding  model,  and 
a  virtual  crystal  approximation.  In  all  three  cases,  the  same  result 
is  obtained.  An  analogous  situation  occurs  if  the  energy  gap  is  due 
to  a  crystalline  structure  distortion  to  lower  symmetry.  In  par¬ 
ticular,  the  pairing  of  ions  in  a  one— dimensional  crystal  is  analyzed. 
The  relationship  between  gap  and  f*ee  carrier  concentration  is  de¬ 
rived  for  the  two  cases  of  a  contact  interaction  and  a  Mathiea  inter- 


action.  The  results  found  are  quite  similar  to  those  when  the  gap 
is  caused  by  antiferromagnetism. 

These  relationships,  together  with  the  application  of  Fermi- 
Dirae  statistics  to  the  conduction  electrons,  enable  us  to  calcxilate 
the  electrical  conductivity  as  a  function  of  temperature  in  two  opposing 
limits,  the  effective  mass  approximation  and  the  limit  of  narrow  bands 
It  is  found  in  both  cases  that  a  singular  increase  in  the  carrier  con» 
conti'ation,  and  tlua  the  conductivity,  occurs  at  a  giv^en  temperature. 
This  leads  to  the  disappearance  of  the  energy  gap,  and  therefore  to 
a  semiconductor-to-metal  transition. 

The  transition  temperature  is  evaluated  in  terms  of  the  zero 
tempi.raiure  energy  gap,  A  number  of  results,  relating  experimentally 
measurable  quanlit  es  such  as  the  pressure  cceflicientf  u£  the  transi- 
t  on  temperature  and  energy  gap,  are  derived. 

The  experimental  results  dealing  with  the  crystals  which 
exhibit  scmiconductor-to-metal  transitions  are  presented,  and  the 
predictions  of  the  theory  are  tested.  Very  good  agi  eement  is  obtained. 

The  theory,  taken  together  with  the  available  experimental 
information,  is  used  to  suggest  possible  band  schemes  for  each  of 
the  materials  under  investigation.  It  is  found  that  the  crystalline 
symmetry  in  every  case  is  compatible  with  an  energy  gap  which  arises 
from  either  a  pairing  of  cations  or  antiferromagnetism. 

Finally,  the  effects  on  the  theory  of  spin-disorder  scattering, 
formation,  and  non— stoichiometry  are  discussed  quantitatively. 


. . 
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Chapter  I 
INTRODUCTION 

A.  Transition  Metal  Oxides  which  are  Insulators  at  all 
Temperatures 

The  transition  metal  oxides  provide  a  striking  example  of  the 
inadequacy  of  simple  band  theory  when  an  attempt  Is  made  to  predict 
the  electrical  transport  properties  of  crystals.  Mcst  of  these  oxides 
are  msulators  [  ij  ,  despite  the  apparent  presence  of  a  partially  filled 
3d  band.  The  materials  which  remain  non-conducting  at  all  tempera¬ 
tures  shall  be  referred  to  as  the  NiO  class.  These  include  NiO,  MnO, 
and  The  experimental  situation  has  been  reviewed  by  Morin[2] 

and  by  jcnker  and  van  Houten  [3];  the  only  Information  of  direct  inter¬ 
est  tc  us  here  is  that,  where  the  magnetic  structure  has  been  deter¬ 
mined,  these  oxides  ate  all  antiferromagnetic,  and  they  are  insulating 
both  below  and  above  the  Neel  temperature.  We  are,  however,  inter¬ 
ested  'n  discu&smg  the  proposed  theoretical  models. 

Many  attempts  h^ave  been  made  to  explain  the  non-conducting 
nature  of  these  materials.  DeBoer  and  Verwey  [  1]  assumed  that  the 
movement  of  electrons  between  neighboring  cations  is  impeded  by  a 
high  potential  barrier.  Thus  the  3d  electrons  are  essentially  localized, 
conductivity  being  possible  only  via  tuiineling  of  electrons  through  the 
ha^rriers.  Verwey  later  extended  the  model,  showing  that  the 
activation  energy  was  due  to  the  ionization  of  the  cations. 

Heikes  and  Johnston  j^5j  noted  that  the  observed  ionization  energies 


are  not  large  enough  to  account  for  the  measured  activation  energies  of 
the  NiO  class,  which  are  of  the  order  of  0.  5  eV.  They  suggested  that 
the  activation  energy  should  be  associated  with  the  mobility,  rather 
than  with  the  carrier  concentration.  The  activation  energy  in  this  theory 
arises  from  Landau  trapping  [6j,  which  occurs  when  an  electron  is 
localised  around  a  lattice  point  in  a  polar  crystal  long  enough  to  cause 
a  deformation  of  the  crystal.  The  electron  is  then  able  to  form  a  bound 
state  with  the  potential  well  due  to  the  polarization.  Heikes  and  Johnston 
ascribe  the  conductivity  to  a  hopping  or  diffusion  of  electrons  from  one 
site  to  a  neighboring  site,  which  will  occur  when  the  lattice  polarization 
around  both  is  momentarily  the  same. 

Yamashita  and  Kurosawa  [7]  attempted  to  work  out  the  theory, 
starting  with  a  HelUer- London  approach.  They  found  it  necessary  to 
assume  that  the  electrons  arc  iccalized  around  the  cations,  and  were 
then  able  to  show  that  the  electrons  could  be  trapped  by  the  resulting 
lattice  polarization,  and  that  conductivity  could  occur  by  means  of  equal- 
site  hopping.  However,  no  reason  was  given  for  the  localization  in 
the  first  place. 

Anderson  [8]  ,  iJiking  up  a  suggestion  of  Mott  [9]  ,  put  the  theory 
on  a  somewhat  firmer  footing  ty  pointing  out  that  the  localization  could 
be  due  to  the  large  correlation  energy  we  should  expect  in  a  low  density 
system.  Anderson  called  U  the  energy  necessary  to  remove  one 
electron  from  a  given  cation  and  put  it  on  another  cation  wh  h  is  far 
away.  U  is  thus  an  intra-ionic  Coulomb  integral,  and  can  be  estimated 
from  the  free  cation  ionization  energi  Anderson  approximates  U  as 
of  the  order  of  10  eV,  although  there  are  reasons  to  believe  thic  is  too 


high.  If  U  is  largei*  than  half  the  bandwidth,  it  outweighs  the  kinetic  _ 
energy  loss  which  co  ild  be  accomplished  by  spreading  through  the 
crystal,  and  the  electrons  would  remain  localized  Although  this  argu¬ 
ment  appears  to  justify  the  procedure  of  Yamashita  and  Kurosawa,  it, 
too, is  only  a  self-Consistent  hypothesis,  since  the  presence  of  many 
free  electrons  would  reduce  U  significantly  by  means  of  screening  Thus, 
the  large  order  of  magnitude  given  by  Anderson  for  U  is  accurate  only 
when  the  electrons  are  known  to  be  localized,  and  it  is  possible  for  a 
metallic  state  to  be  self-consistently  lower  in  total  energy. 

Holstein  [lO]  has  applied  the  theor-,'  of  polaron  formation  of 
Frohlich  [11]  to  the  situation  in  the  transition  metal  oxides.  A  polaron 
is  e  bound  state  of  an  ectron  and  its  associated  lattice  deformation 
which  can  be  treated  as  a  quasiparticle  moving  througl  the  crystal  . 
Holstein  considered  in  particular  the  "small"  polaron,  whose  dimen- 
sic.  s  arc  of  the  order  of  a  lattice  constant.  The  polaron  will  be  "small" 
if  the  electronic  bandwidth  is  much  leas  than  the  maximum  binding 
energy  of  the  polaron,  a  constant  which  is  proportional  to  the  square  of 
the  strength  of  the  electron-phonon  interaction.  Holstein  found  that  at 
low  temperatures  there  is  sufficient  overlap  to  form  a  polaron  band 
through  which  the  electron,  together  with  its  lattice  deformation,  can 
move.  An  unfilled  polaron  band  carries  a  current  which  decreases 
exponentially  with  increasing  temperature.  Since  the  width  of  the 
polaron  band  is  also  an  exponentially  decreasing  function  of  tempera¬ 
ture,  and  since  the  polaron  states  have  a  finite  lifetime,  eventually 
the  point  is  reached  when  the  bandwidth  is  smaller  than  the  uncertainty 
in  energy  associated  with  the  lifetime,  and  the  bands  lose  their  physical 
meaning.  Above  this  temperature,  estimated  by^  Holstein  as  half  the 
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Debye  temperature,  conduction  occurs  when  the  interaction  oi  phonons 
with  the  deformation  potential  causes  ionic  rearrangements  at  neigh-  & 
boring  sites,  enabling  the  polaror  to  jump  from  one  to  the  other.  This  ^ 
mechanism  is  just  thermally  activated  hopping  of  polarons,  * 

An  entirely  different  approach  was  suggested  by  Slater  [12]  .  t 
He  retained  the  band  picture  and  suggested  that  the  antiferromagnetic 
ordering  can  bring  about  an  extra  band  splitting.  Slater’s  argument 
goes  as  follows;  According  to  the  Hartree-Fock  approximation,  an 
electron  will  be  repelled  much  less  strongly  by  electrons  with  p>arallel 
spins  than  by  those  with  antiparallel  spins.  Since  the  crystal  is  anti- 
ferromagnetic,  each  electron  feels  a  potential  which  has  the  periodicity 
of  the  superlattice  rather  than  that  of  the  lattice.  Such  a  potential  will 
spilt  the  first  Briliouin  zone  in  half,  each  half  containing  1/2  N  V  states 
for  each  direction  of  spin,  where  N  is  the  concentration  of  cations  per 
unit  volume  and  V  is  the  volume  of  the  crystai.  Thus  even  when  there 
is  one  electron  per  cation,  the  first  zone  could  in  principle  be  filled, 
the  second  empty.  Two  difficulties  with  this  theory  immediately  ^ 

emerge.  The  3d  band  should  be  highly  degenerate,  particularly  for  g 
cvb’c  lattices,  such  as  MnO  and  CoO,  and  Slater  does  not  attempt  to  ^ 
show  how  such  crystals,  containing  5  and  7  3d  electrons  per  cation, 
respectively,  can  have  anything  but  a  partially  filled  valence  band,  g 
even  with  the  antiferromagnetic  splitting.  Furthermore,  Slater  gives  ^ 
no  explanation  of  the  non-conductivity  of  the  NiO  class  of  materials 
above  the  Neel  temperature. 

Recently,  Hubbard  [13]  considered  another  band  approach  whiclf 
removed  the  latter  objection.  Hubbard  showed  that  electronic  corre-  3 


cuiuu  i-esiiit  tn  the  splitting  of  a  iiarrow  band  into  two  sub-bands, 
.A-di-gyus  to  the  splitting  due  to  exchange  energy  in  Slater's  model, 
bat  independent  of  the  magnetic  ordering  of  the  system.  Thus  the 
crystals,  if  insulating  at  T  =  0,  should  remain  insulating,  even  above 
the  Neel  temperature.  This  correlation  splitting  is  essentially  a 
manifestation  of  the  Anderson  effect  [8]  ,  in  that  it  represents  a  partial 
localization. 


B.  Transition  Metal  Oxides  which  Exhib't  Insulator-to-Matal  Transitions 

The  lack  of  electrical  conductivity  is  not  the  only  enigma  pro¬ 
vided  by  the  transition  metal  oxides.  Another  subclass,  referred  to 
here  as  the  class,  consists  of  several  oxides  of  vanadium  and 

titanium,  which  are  non-conducting  at  low  temperatures,  but  undergo 
a  transition  to  a  metallic  state  above  a  critical  temperature.  These 
are  the  materials  with  w'hich  we  are  concerned  in  this  w'ork 

The  first  hint  of  unusual  behavior  in  this  class  of  oxides  came 
when  Perakis  [I4j  discovered  that  the  magnetic  Susceptibilities  of  VO^ 
and  each  went  through  a  sharp  transition  at  a  given  temperature. 

A.nderson  [15]  found  a  specific  heat  anomaly  in  in  the  vicinity  of 

150  K,  with  a  latent  heat  of  approximately  700  calories  per  mole 
measured  over  a  20^K  temperature  inter^'al.  Foex  [16]  investigated 
the  electrical  properties  of  and  discovered  a  sharp  increase  of 

a  factor  of  10  in  conductivity  at  the  temperature  of  the  specific  heat 
anomaly.  He  also  reported  a  volume  contraction  at  the  same  paint. 

Fohx  and  Loners  [17]  found  a  transition  in  Ti^O^  at  about  475^K,  the 
electrical  conductivity  jumping  by  a  factor  of  10,  Jaffray  and  Dumas  [18| 
verified  the  transition  in  and  also  reported  an  insulator-to-metal 
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transition  as  well  as  a  small  contraction  in  volume  in  VO_  at  the  tem- 

z 

perature  of  the  magnetic  susceptibility  anomaly,  Pearson  [19]  noted 
an  absorption  edge  of  0.  15  eV  in  a  sample  of  powdered  Ti,jOj  in  KBr, 
but  did  not  obtain  a  break  in  the  conductivity  curve.  He  did  find  that 
the  conductivity  as  well  as  the  lattice  constmts  were  rapidly  changing 
in  the  vicinity  of  500°K. 

The  situation  with  regard  to  the  electrical  properties  of  the 

lower  oxides  of  titanium  and  vanadium  was  resolved  by  the  work  of 

Morin  [20],  who  performed  measurements  of  conductivity  as  a  function 

of  temperature  on  single  crystals  of  Ti^O^,  TiO,  VO,  and  VO^. 

His  results  are  given  in  Fig.  I-l.  As  can  be  seen,  all  these  materials 

with  the  exception  of  TiO  are  semiconducting  at  low  temper  .tures,  but 

undergo  a  transition  to  the  metallic  state  at  a  given  temperature,  which 

we  shall  call  T  . 

o 

The  experimental  situation  since  Morin's  work  has  been  reviewed 
in  great  detail  by  Fcinleib  [21] ,  and  we  shall  at  this  point  mention 
only  three  important  discoveries  of  recent  years.  Yahia  and  Fred- 
erikse  [22]  were  able  to  measure  the  Ifell  coefficient  in  Ti^O^,  and 
found  that  the  Hall  activation  energy  was  virtually  the  same  as  the  con¬ 
ductivity  activation  energy.  Abrahams  [23]  performed  neutron  dif¬ 
fraction  experiments  on  Ti^O^,  and  determined  that  it  was  antiferro¬ 
magnetic  with  a  small  moment  of  0.  2  Bohr  magnetons  per  cation, 
the  Neel  temperature  being  in  the  vicinity  of  the  semiconductor-to-metal 
transition.  Most  recently,  Feinleib  [24]  performed  optical  measure¬ 
ments  on  single  crystals  of  stoichiometric,  semiconducting  and 
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observed  an  absorption  edge  at  approximately  0.  1  eV. 

The  first  attempt  at  an  explanation  of  the  sharp  discontinuities 
in  conductivity  was  given  by  Morin  [20]  ,  who  adapted  the  theory  of 
Slater  [12] .  If  these  materials  were  semiconducting  due  to  a  band 
splitting  arising  from  antiferrornagnetismi  then  a  transition  to  the 
metallic  state  would  be  expected  at  the  Neel  temperature. 

Callaway  [25]  made  the  Slater-Morin  model  a  little  more  quant 
tative  by  studying  the  energy  band  structure  of  a  body-centered  cubic 
antiferromagnet.  Considering  only  the  first  Fourier  component  of 
the  exchange  potential,  a  major  simplification,  he  found  that  an  in¬ 
sulating  state  exists  whenever  an  interaction  parameter  (proportional 
to  the  effective  mass,  to  the  strength  of  the  exchange  potential,  and 
to  the  square  of  the  lattice  constant)  is  sufficiently  large.  Neither 
Callaway  nor  Morin  discussed  the  nature  of  the  semiconductor-to- 
metal  transition  beyond  noting  that  the  band  gap  should  disappear  at 
the  Neel  temperature. 

There  are  three  difficulties  with  the  Slater-Morin  theory, 
outside  of  the  fact  that  it  has  never  been  quantitatively  applied  to  the 
oxides  of  titanium  and  vanadium.  FirsUy,  the  existence  of  antiferro¬ 
magnetism  has  been  demonstrated  only  in  Ti^O^,  and  even  In  that 
material  the  antiferromagnetic  moment  is  extremely  small.  Secondly, 
no  model  for  the  structure  of  the  degenerate  3d  bands  has  been  pre¬ 
sented  which  explains  how  Ti^O,.  and  VO^  can  be  semiconduct¬ 

ing  even  with  the  antiferromagnetic  splitting.  Finally,  it  has  not 
been  proven  that  the  energy  gap  must  vanish  with  the  disappearance 
of  long-range  order;  it  is  conceivable  that  the  large  amount  of  short 


g 

range  order  present  until  two  or  three  times  the  Neel  temperature  is 

sufficient  to  maintain  the  effective  double  periodicity  seen  by  the 

slowly  moving  3d  electrons,  and  thus  maintain  an  energy  gap 

A  different  approach  to  the  problem  was  presented  by  Mott  [26], 

who  proposed  that  there  would  be  a  sharp  transition  from  an  insulatmg 

to  a  conducting  state  at  a  critical  value  of  the  lattice  parameter  R 

■  c 

Mott  argued  that  a  free  electron  and  a  free  hole  attract  each  other  with 

a  Coulomb  potential,  V  (r)  =  -e^/er  ,  and  form  a  bound  state,  similar 

to  an  exciton,  from  which  neither  can  participate  in  conduction.  However, 

if  a  large  number  of  free  carriers  exist,  the  attraction  between  an  electron 

and  a  hole  is  a  screened  Coulomb  potential,  V(r)  -  -e^e*^^/cr.  When  q 

becomes  sufficiently  large,  this  potential  becomes  too  weak  for  the 

formation  of  bound  states,  and  a  discontinuous  increase  in  the  number 

of  free  carriers  brings  about  a  transition  to  the  metallic  state, 

Mctt  [27]  assumed  Fermi-Thomas  screening,  and  roughly  estimated 

R  as  : 
c 

o 

R  Z  — J - 2.  1  A 

C  i 

m  /m 

where  z  is  the  number  of  3d  electrons  per  cation,  e  is  the  dielectric 
constant,  and  m  is  the  effective  mass  of  an  electron  in  the  conduction 
band.  Mott  did  not  attempt  to  apply  quantitatively  his  theory  to  the 
vanadium  and  titanium  oxides,  but  it  is  possible  that  a  crystalline  con¬ 
traction  at  a  given  temperature  could  indeed  lower  the  lattice  param¬ 
eters  in  Such  a  way  as  to  bring  about  a  transition  to  the  metallic  state. 

Since  there  is  evidence  of  phase  transformations  at  or  near  T  in  all 

o 

four  materials  which  undergo  semiconductor-to-metai  transitions, 


it  ii  possible  that  this  inod^l  applies*  Hewev^ef *  the  teniperature 
variation  of  the  lattice  constants  often  produces  changes  greater  than 
those  which  occur  during  the  phase  transformation,  but  In  the  “pposfte 

direction.  Furthermore,  the  expression  derived  for  R  should  favor 

c 

metallic  conductivity  in  VO.  where  a  =  3,  over  that  in  TiO,  where 
a  the  lattice  parameters  and  ionic  radii  of  the  two  crystals  being 
roughly  the  same*  Actually,  the  reverse  is  the  case,  TiO  being 
the  metal  at  all  temperatures.  Other  objections  to  this  theory  wilj 
be  presented  in  Chapter  VI. 

Goodenough  [28]  accepted  Mott's  hypothesis  of  a  critical  cation 
separation,  and  proposed  that,  due  to  the  anisotropy  of  3d  orbitals 
and  complex  crystalline  structures,  localized  and  collective  electrons 
could  be  present  simultaneously  in  a  given  crystal.  By  considering 
the  entire  class  of  transition  metal  oxides,  he  was  able  to  estimate 
R^  ^  2*  9  A.  Noting  that  the  lattice  parameters  in  the  V  O,  class 

j 

are  smaller  than  R^,  Goodenough  [29]  suggested  that,  due  to  direct 
cation-cation  interaction,  all  would-be  conduction  electrons  could  be 
trapped  in  homopolar  bonds  at  low  temperatures.  Above  a  critical 
temperature,  the  free  energy  associated  with  a  metallic  state  may 
become  lower  than  that  of  the  bonding  state,  because  of  the  extra 
entropy  associated  wtth  the  continuum  of  electronic  states  in  the  metal¬ 
lic  phase,  and  a  transition  would  then  occur,  Goodenough  [30]  did 

hypothesis  to  the  oxides  of  vanadium  and  titanium,  and  was 
able  to  account  for  many  of  the  previously  unexplained  symmetry 
changes.  However,  the  theory  does  not  lend  Itself  to  quantitative  in- 
vestigation,  and  has  other  difficulties  which  will  be  detailed  in 


Chapter  VL 


Two  recent  theoreUcal  discussions  are  quite  relevant,  alihough 
neither  has  been  applied  to  the  class  of  materials.  Kohn  [31] 

eveloped  a  characterization  of  the  insulating  state,  which  essentiaxjy 
makes  quantitative  part  of  Mott's  original  proposal  [9]  .  Kohn  con¬ 
sidered  a  ring  of  hydrogen  atoms*  and  shewed  that  if  the  separai  on 
between  nearest  neighbors  is  sufficiently  large,  one  can  express  each 
wave  function  as  a  sum  of  localized  functions  which  do  not  overlap, 
and  such  wave  functions  correspond  to  an  insulating  system,  Kchn  also 
investigated  the  Mott  transition  [27]  by  considering  a  ferromagnetic 
simple  cubic  array  of  hydrogen  atoms.  If  the  potential  between  a 
spm-up  hole  and  a  spin-down  electron  were  a  Coulomb  attraction,  an 
insulating  spin-wave  state  is  lowest,  independent  of  the  wtrakness  of 
the  interaction.  However,  for  a  delta  function  interacticn,  a  cr  ’ica 
strength  exists,  below  which  only  a  continuum  of  states,  characteriz¬ 
ing  a  metal,  is  present.  Thus,  an  insulator -to -metal  trans't.on  can 
occur  at  a  critical  value  of  the  lattice  parameter. 

Hubbard  [32]  extended  his  own  theory  of  correlations  in  narrow 
bands,  and  found  that  at  a  critical  ratio  of  the  bandwidth  ,  to  the 
intraionic  correlation  energy,  U- 

^  -  2 
/3 

the  energy  gap  due  to  correlation  has  shrunk  to  zero,  and  a  semicon- 
ductor-to-metal  transition  occurs.  However,  the  energy  gap  shrinks 
continuously  as  the  ratio  of  tc  U  increases,  so  that  slightly  below 
the  critical  value,  the  material  has  an  nfmltesimal  band  gap,  while 
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just  above  the  transition  point,  the  density  of  states  at  the  Kermi 
surface  is  negligib^.y  small,  so  that  the  material  is  a  poor  metal. 

This  type  of  transition  is  unlikely  to  be  responsible  for  a  change  in 
conductivity  of  a  factor  of  10^  at  T  ,  as  is  found  in  V,0,  and  VO. 

It  IS  possible  that  Hubbard's  neglect  of  interactions  of  electrons  on 
different  ions  camouflageU  the  sharpness  of  the  transition,  but  until 
worked  out  explicitly,  this  remains  a  speculation.  If  the  nature  of 
the  disappearance  of  the  energy  gap  is  as  Ihibbard  has  found,  then 
Mott's  hypothesis  of  a  sharp  increase  in  the  number  of  carriers  is 
incorrect. 

ft  is  clear  that  all  theories  of  insulator-to-metal  transitions 
have  serious  deficiencies  when  applied  to  the  class  of  materials. 

In  this  work,  we  shall  attempt  to  present  a  model  which  overcomes 
the  objections  discussed  above.  We  shall  retain  the  band  picture, 
and  assume  the  non-condu  ting  state  of  these  materials  is  that  of  a 
normal  semiconductor,  having  a  filled  valence  band  separated  from 
an  empty  conduction  band  by  an  energy  gap.  In  Chapter  II,  we  shall 
show  how  such  a  gap  can  arise  in  the  transition  metal  oxides  from 
antiferromagnetism  and  also  from  a  crystalline  structure  distortion 
to  lower  symmetry.  It  will  be  demonstrated  that  in  these  two  cases 
the  energy  gap  will  shrink  as  carriers  are  excited  across  it,  and  the 
decrease  in  the  gap  will  be  quantitatively  estimated.  A  thermodynamic 
argument  will  be  presented  to  relate  this  decrease  to  the  observed 
pressure  coefficient  of  the  gap.  In  Chapter  HI,  the  theory  of  con¬ 
ductivity  will  be  worked  out  in  two  limits,  the  effective  mass  approxi¬ 
mation  and  the  limit  of  narrow  bands.  We  shall  demonstrate  that  a 


sharp  semiconductor-tQ-nietaJ  traRs^uon  occurSj  and  caiculate  the 
transition  teriperature  in  terms  of  obser%'able  quantities*  In  Chap¬ 
ter  IV,  the  theory  of  Chapters  II  and  III  will  be  applied  to  the  exp,ri  - 
mental  data  on  the  oxides  of  titanium  and  vanadium.  Certain  rela¬ 
tionships  predicted  by  the  theory  will  be  critically  ter/fed  quantitdit-: ve  y 
Chapter  V  will  consist  of  suggested  energy  band  models  for  the 
materials  investigated,  based  on  a  combination  of  the  experimental 
data  available  and  the  predictions  of  the  theory.  Finally,  in  Chap¬ 
ter  VI,  we  shall  discuss  the  effects  of  spin-disorder  scattering  and 
broadening,  polaron  format’on,  and  non-s^oichiometry,  as  well  as 
demonstrate  additional  inconsistencies  with  some  of  the  other  models 
v>  hich  have  been  presented. 
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Chapter  II 

DEPENDENCE  OF  ENERGY  GAP  ON 
CARRIER  CONCENTRATION 

A.  General  Hypothesis 

Consider  an  intrinsic  semiconiuctor  for  which  the  top  of  the 
valence  band  is  separated  from  the  bottom  of  the  conduction  band  by  an 
energy  gap,  E  .  In  general,  E  depends  on  the  concentrat-cn  of  carriers 

D  S 

in  the  conduction  band,  n,  and  on  the  temperature,  T. 

E  =  E  (n,  T) . 

g  g 

At  low  temperatures,  the  concentratjon  of  carriers  is  aisc  small,  and 
we  can  write: 

E  "  E  -  cT  -  8n 

g  -TO 

3  E  3  E 

where  is  the  gap  at  T  =  0,  a  )^,  and  p  •  Although 

the  term  linear  in  T  is  respon  jie  for  the  major  part  of  the  decrease 
in  band  gap  at  very  low  temperatures,  it  does  not  contr:ibute  to  the 
semiconductor-to-metal  transition  and  therefore  will  be  dropped.  In 
Appendix  A,  we  shall  show  that  only  a  small  error  is  introduced  into 
the  calculation  by  ignoring  this  term.  We  are  left  with: 

E  -  E  -  pn  {2.  1) 

g  go  1  *  '  . 

which  is  our  fundamental  relation.  The  remainder  of  this  chapter  will 
be  devoted  to  demonstrating  the  applicability  of  Eq.  (2.1)  In  two 
particular  situations,  where  the  energy  gap  is  due  to  antiferromag¬ 
netism  and  where  the  gap  arises  from  a  crystalline  structure  distortion 
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to  lower  symmetry.  We  shall  evaluate  p  for  both  of  these  cases,  and 
show  that  Eq.  (Z.  1)  remains  valid  as  n  becomes  relatively  large. 


B.  Thermodynamic  Argument 

In  this  section,  we  shall  present  a  thermodynamic  crlculation 
of  the  change  in  energy  gap  of  a  semiconductor  with  the  concentration 
of  excited  carriers.  This  will  provide  us  with  a  general  expression 
for  p  in  Eq.  1).  The  first  part  of  this  argument  follows  closely  a 
recent  paper  by  Figielski  [33]  . 

The  differential  form  of  the  Gibbs*  free  energy  for  a  system 
where  the  number  of  particles  may  vary  is; 

dG  =  -S  dT  +  V  dP  +  E  p.  dN.  (2.21 

}  J  } 

th 

where  is  the  number  of  particles  in  the  j  phase  and  p ^  is  the 
chemical  potential  of  the  pliase.  Treating  electrons  in  the  valence 
and  conduction  bands  as  different  phases  and  ignoring  inner  orbitals, 
we  can  write; 


Ep.dN.  =p  dN  +p  dN  .  (2.3) 

J  j  ^c  c  V 

For  an  intrinsic  semiconductor,  N  +  N  =  constant,  and  (2.  3)  becomes: 

VC 

Ep.  dN.  =  (p  -  p  )dN  (2.4) 

where  N  is  the  number  of  carriers. 


a 


2 

a  G 


Since  =  3  .  (2.2)  and  (2.4)  give: 


(1^) 

'FN 'P,  T 


V  'SIT'P,  T  ~  V  T 


ap 

V 

^TP^n.  T 


(2.5) 


where  n  =  N/V  is  the  concentration  of  carriers.  The  chemical 


potential,  must  bt  calculated  from  the  relation. 


ae  E  -  ji 

=  /  p  JE)  f(-r^)  dE 


where  is  the  concentration  of  electrons  in  the  valence  band,  p  (E)  is 

the  density  of  states  m  the  valence  band,  and  f  (x)  =  [e^+  ij‘^  The  chenncai 

potential  ,  is  determined  from  an  entirely  analogous  relation.  It  is 

clear  that  at  T  =  0,  and  =  E^,  where  E^  is  the  energy  of  the  top 

of  the  valence  band  and  E^  is  the  energy  of  the  bottom  of  the  conduction 

band.  In  general,  is  the  change  in  free  eneigy  when  an  electron  is 

removed  from  the  valence  band  and  placed  in  the  conduction  band.  It  can 

be  shown  that  differs  .rom  by  small  terms  proportional  to  kT, 

which  we  can  essentially  neglect  Substituting  this  result  in  {2.  5): 

18V  ^  ^  ®  ^ 

V  'ffr’p.T  "  <TF*n,T  -  <-rP  >n.T  =  'Tli'n.T 

where  E^  *  E^  •  E^_  is  the  energy  gap. 

Using  the  thermodynamic  relations; 


T  ■  ■  W 


O  Zs 


where  K  is  the  isothermal  compressibility,  we  find: 

,1^,  .  1  ,8V, 

'tt’p.t  -  ■  tv  '-Tpin.T  *Fir>p,T  <Tn^’v,  T 

Thus,  (2*6)  and  (2,  7)  yield: 


(2,7) 


5E  1  8E  2  8E 

*-3Tf'p,T  " -t  *<TK^>V.T- 


(2.8) 


This  is  the  general  thermodynamic  expression  for  p.  The  first  term 
on  the  right  represents  the  contribution  to  0  resulting  from  changes 
in  the  volume  of  the  crystal.  This  term  can  be  evaluated  easily  from 
the  experimentally  measurable  quantit  es,  y,  K,  and  is  always  posi¬ 
tive.  The  other  contribution  to  0  is  an  explicit  dependence  of  the  gap 
on  carrier  concentration  at  constant  volume,  and  can  have  either  sign 
This  thermodynamic  argument  shows  that  whenever  there  is  a 
pressure  dependence  of  the  energy  gap  of  a  semiconductor,  the  gap 
must  also  depend  on  carrier  concentration.  Since  the  pressure  co¬ 
efficients  of  the  gap  in  and  VO,  two  of  the  materials  with  which 

we  are  especially  concerned,  are  anomalously  large,  we  expect  a 
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relatively  large  decrease  of  gap  with  carrier  concentration  in  these 

materials*  But  this  argument  demonstrates  only  that  E  depends 

on  the  number  of  excited  carriers.  It  does  not  say  anything  about  the 

validity  of  Eq.  {Z.  1)  as  n  becomes  fairly  large,  nor  does  it  indicate 

the  microscopic  reasons  for  such  a  variation  in  energy  gap.  Therefore, 

we  now  turn  to  specific  models  for  which  we  can  calculate  express  ons 

for  E  as  a  function  of  n.  These  calculations  will  provide  us  also 
g 

with  expressions  for  p  which  can  be  tested  experimentally. 


C.  Antiferromagnetism 

Consider  an  antiferromagnetic  crystal  which  can  be  described 
by  Bloch  wave  functions.  Assume  that  the  crystal  is  an  insulator  at 
T  -  0  because  of  the  splitting  of  the  first  Brillouin  zone  by  the  doubly- 
periodic  exchange  potential.  In  other  words,  we  have  an  empty  con¬ 
duction  band  which  begins  a  distance  E  above  the  filled  valence  band, 
with  E  being  a  measure  of  the  exchange  energy.  The  lower  band 
refers  to  wave  functions  whose  amplitudes  are  large  at  the  sublattice 
positions  of  the  electron  under  consideration,  whereas  the  upper  band 
wave  functions  have  large  amplitudes  at  the  positions  ot  the  sublattice 
of  opposite  spin.  As  the  temperature  is  increased  from  T  =  0,  the 
upper  band  becomes  thermally  populated.  When  an  electron  is 
excited  across  the  energy  gap,  the  net  magnetization  on  either  sub¬ 
lattice  decreases,  and  thus  the  gap  decreases  with  increasing  con¬ 
centration  of  carriers.  Thus  a  relationship  like  Eq.  (2.  1)  can  be  ex¬ 
pected  to  hold.  In  this  section,  we  shall  determine  for  how  large  n 
Eq.  (2.1)  remains  valid  and  also  calculate  the  value  of  p.  We  shall 


I 
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begin  with  a  one -dimensional  model,  since  it  can  be  solved  in  terms 
of  simple  band  parameters  and  the  solution  exhibits  quantitatively 
the  points  discussed  above.  We  shall  then  repeat  the  calculation  for 
the  case  of  a  three-dimensional  face-centered  cubic  crystal  (such  as 
VO)  using  the  approximation  of  tight  binding.  Finally,  we  shall 
evaluate  P  in  the  limit  opposite  to  that  of  elementary  Bloch  theory, 
employing  a  virtual  crystal  approximation. 


1.  One -Dimensional  Model 

Consider  a  one-dimensional  antiferromagnetic  crystal,  with  N 
ions  spaced  a/2  apart.  Assume  each  ion  contains  one  electron  outside 
the  core  and  take  the  ion  at  the  origin  to  have  spin  up.  The  exchange 
potential  has  periodicity  a,  and  thus  the  first  Brillouin  zone  is 
-ir/a  <  ir/a.  In  accordance  with  L.  C.A.O.  theory,  the  Bloch  wave 
functions  are: 


xl  -  (|^)  c,,(k)  Z  e^^"^  a(x  -  na)  a 

1  Z  II  n 


x)  =  {^)  ^  a(x  -  na)  a 


(2.  141 


,  ,  N  .  '  ^  /I  I  V  l/2)a  r  i  .  1  1  a 

+  {^)  ^22^^^  f  ®  a[x  -  (n  +  2  )a]  p  . 


Assuming  the  ionic  functions  do  not  overlap,  the  secular  equation  is: 


det  !H  (k)  -  E  6  1=0 

mn  mn 


(2.  15) 


where: 


H|.(k)=  Z  coskna  /  dx  a(x  -  na)  a{x) 

*  *  n 

H22W  =  £  coskna  /  dx  a[x  -  t  i)a]  H(x)  a(x  *  ^a) 

^  CO*  k  (n  -  ^  )  a  /  dx  a{x  -  na)  H{x)  a{x  -^a)  * 

Let: 

«^  =  /  dx  a(x  -  na)  H(x)  a(x) 

=  /  a[x  -  (n  +  i)a]  H(x)  a(x  -  ia) 
s/  dx  a(x  -  na)  H(x)  a(x  -  ^a)  . 


With  these  definitions,  and  assuming  nearest-neighbor  interactions 
only,  (2.  15)  becomes; 


Zp 


,  coa^ka 
{e  +  c  )  E 


+ 


2^^  cos  ka  i 

%-E  I 

a  c  - 4p  cos  ^  ka 


where  we  have  dropped  the  subscript  o.  Solving  the  secular  equation: 

etc 


®  ^  2  1, 
V' — 2 — 1  *  ^ka 


E{k)  =  < 


(2.  16) 


+  a  /  /  «  -  fl  .^2  2  1, 

~Z —  "  V  ^ — 2 — ^  ^  2 


Consider  an  electron  with  spin  down;  thus  a  >  c  .  Also, 
consider  the  narrow  band  case,  where  p  «  la  -  a  I  .  Then  (2.  16)  may 


be  written; 


E^(k)  i  e  + 


4p  cos  ^  KA 


E  {k)  =  c 

V 


2  1  u 

cos  ”  k  a 


a  ~  it 


where  E  (k)  refera  to  the  conduction  band^  E  (k)  to  the  valence  band, 
c  V 

From  {2.  17),  we  can  immediately  write  down  the  important  band 
parameters: 

Energy  Gap  s  E  ~  a  -  a 
S 

Bandwidth  h  4p^/(c  -  a)  =  4p^/E  (2.  18) 

8 

*  2  2 

Effective  Mass:  m  =m,  =m  =2li  /E,  a  , 

eh  b 


In  the  one-dimensional  case,  the  energy  gap  is  exactly  the  difference 
in  exchange  energy  between  a  spin  up  and  a  spin  down  electron. 

Evaluating  the  wave  functions  corresponding  to  electrons  in 
the  valence  and  conduction  bands,  we  obtain: 


1  - 


■5a2  2  1, 

2p  cos  ^  ka 
(c  -  c  ) 


+  = 
c 


2S  cos  ^  k  a 


= 


-a  i  1  \ 

-p  cos  2  ^  ^ 

■  "31 

a  -  o 


1  - 


?  2  1 
Zr  cos  y  ka 


Equations  (2.  19)  lead  to: 


c  \  ^  ^  2  1, 

^z>c  =  2  -  ^  2 

g 


(2.  19) 


c  V  1  b  2  i  , 

S  >  -OS  irka 

z/  V  2  E  2 

g 
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EcjiiiitloTiS  (2.  20)  Can  Dt;  liscd  to  give  iis  srsmf?  insight  into 
antiferromagnetism  li  the  bandwidth  is  small  compared  to  the  band 
gap,  then  the  average  spin  on  each  subiattice  will  be  close  to  4  i/2, 
and  neutron  diffraction  experiments  should  show  an  antiierromagnetic 
moment  near  one  Bohr  magneton  per  catiton  How^ever.  as  the  band- 
w'ldth  gets  to  be  of  the  order  of  magnitude  of  the  energy  gap,  much 
lower  values  for  the  sublattice  magnetization  can  be  obtained. 

In  connection  with  this,  it  is  useful  t  investigate  the  exact 
solution,  {Z.  16),  rather  than  the  narrow  band  approximation,  {2  17). 
Although  the  energy  gap  remains  exactly  E  =  a  -  a  ,  independent 
of  the  bandwidth,  the  general  solution  for 


E 

g 


1  1  . 


Hence,  as  p  approaches  E  ,  the  real  bandwidth  is  smaller  than  the 

g 

value  given  by  (2.  18)  When  p  =  E  /4  ,  E,  =  0.  21  E  ;  when  8  =  E 

g  ^  i  ‘  g  * 

E.  =  1. 6  E  , 
b  g 

The  exact  solution,  (2,  16),  can  be  most  simply  discussed 
in  terms  of  the  parameter  $  ,  defined  bv: 


tan  6 


4p  cos  ka 


Q  -  a 


The  exact  solutions  for  <  S  >  ,  and  <S  >  can  be 

C  V  •  2  c  '  2  V 

expressed; 
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huSj  the  average  value  of  at  the  bottom  of  the  Conduction 


IS  ^iV/ciV? 


exactly  l/2.  However,  at  the  top  of  the  conduction 


band; 
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As  the  ratio  of  bandwidth  to  band  gap  increases,  the  average  value 

of  at  the  top  of  the  conducting  band  decreases  monotonically . 

When  E.  =  0.  21  E  ,  <S  >  =  0.  34;  when  E,  =  1 .  b  E  ,  <S  >  =  0.  09 
b  g  2  b  g  ’  z 

It  IS  clear  that  the  average  values  of  in  the  valence  band  are 

just  the  negatives  of  those  in  the  conduction  band  in  the  one-dimensional 

case. 

At  T  =  0,  the  valence  band  is  exactly  filled,  the  conduction 
band  completely  empty.  Thus,  the  magnetization  of  the  sublattice 
where  the  spins  are  down  at  T  0  is  given  by^  the  sum  of  <  S  >  _ 

Z  r 

over  all  points  in  k-space: 


\  -  ^  ^ 

^  =  /.  i  - 

Clearly,  as  the  bands  get  wider,  the  total  magnetization  on 
each  sublattice  decreases  considerably  In  the  limit  where  the 
bands  are  much  wider  than  the  gap,  the  average  sublattice  magnet¬ 
ization  goes  to  zero.  In  other  words,  the  antiierromagnetism 


p 


becQtneg  undetectabie  when  the  exchange  energy  is  iniich  smaller 
than  the  overlap  between  nearest  neighbors.  We  shall  return  to 
thin  point  later. 

Calculating  from  the  change  in  spin  as  an  electron  is 

excited  across  the  energy  gap; 

?  I 

AS  (k)  -  1  — = —  cos“^  ka 
Z  d. 


=  1  -  (i  +  cos  ka) 


(2 


Integrating  (2.21)  to  find  the  net  change  in  spin  when  n  electrons 
have  been  excited  across  the  gap; 

IT /a 

A  S^(n)  =  ^  /  dk[  1  -  { 1  +  cos  k  a)] 
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=  £  (1 _ ^ 

N  '  ^  E  '  E 
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since  n/N  =  i  -  ka/w.  Expanding  (2.22)  in  powers  of  n: 


A  S  (n) 
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TT  n  . 
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The  first  term  is  the  major  contribution,  not  only  because  n/U  is 
usually  small,  but  also  because  we  are  considering  the  situation 
where  E^  «  E^, 

The  assumption  is  now  made,  in  view  of  the  fact  that  the 


energy  gap  is  a  direct  measure  of  the  exchange  energy  see  (2.  18)]  , 


{^.241 


E  “ 


g  ^  W  '"e 


=  1)U 


where  is  the  number  of  down  spins  and  U  is  an  average  exchange 
integral.  What  (2.24)  says  is  that  the  exchange  energy  is  porpoitlonal 
to  the  number  of  pairs  of  parallel  spins,  the  assumption  here  is  that 
each  electron  sees  the  spins  of  all  the  other  electrons  la  the  band,  a 
hypothesis  consistent  with  the  spirit  of  band  theory.  We  can  evalua-e 
U  by  noting  that  at  T  -  0: 

U  =  2N  Eg^/(N/2  (N/2  -  1)] 

»  8  E  /N. 

go 

Thus  (2.  24)  becomet 
4E 


:  = - (n^  -  n.)  . 


(2.25^ 


Since  <^S  >  =  (n  -  nrt)/2N  and  N  =  n  +  n^,  where  n  is  the  numbei  oj 

z  ■  '  c  p  '  r  p’  a 

up  spins,  (2.23)  and  (2.25)  show: 
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since  N  "  0(10  )  »  1,  Equation  (2.26)  indicates  that  the  energy  gap 

will  close  down  linearly  with  rumber  of  excited  carriers  for  quite  a 
large  range  of  n.  In  the  limit  of  narrow  bands,  all  terms  of  higher 
order  in  n  are  negligibly  small.  Furthermore,  even  when  the  band 
width  is  of  the  order  of  the  gap,  the  first  correction  term  is  an  order 
of  magnitude  smaller  tl^an  the  linear  term  up  to  values  of  n/N  “  0.  3. 


This  result  tends  to  justify  the  use  of  «  q^ation  (Z.  '  #  In  later  chapters. 
From  (2.26),  the  value  of  p  for  a  one-dimensional  antiferromagnet 
is: 

P  =  4  E  /N  .  (2.27) 

go 

2  Three-dimensional  Model 

Mach  of  the  calculatior>  of  the  three-dimensional  model  is 
entirely  analogous  to  that  of  the  one-dLmensional  problem.  We  return 
to  the  use  of  n  and  N  as  concentrations  of  excited  carriers  and  cations, 
respectively,  and  consider  the  case  of  a  face -centered  cub.  c  lattice 
of  magnetic  ions  with  spins  antiparallel  in  alternating  (H  ll-planes. 

We  shall  assume  that  there  is  one  3d  electron  per  cation  outside  the 
core  and  that  this  eleciron  can  be  descr-bed  by  a  non-degenerate 
Wannier  function  due  to  a  slight  distortion  o»‘  the  lattice.  Th^s  is  a 
good  representation  of  the  situation  in  VO  However,  the  calculation 
is  virtually  unchanged  by  taking  another  structure,  such  as  corundum, 
or  by  choosing  a  larger  number  of  3d  electrons  per  cation  in  a  degenerate 
band.  In  the  case  described,  energy  gaps  due  to  the  exchange  potential 
appear  at  the  surfaces  of  the  reduced  part  of  the  first  Brillouin  zone. 
Thus,  if  we  take  a  cation  of  spin  up  at  the  origin  and  one  of  spin  down 
at  ^  =  (110)  a,  energy  gaps  appear  at  the  planes  k  ■  ^  =  fff/2.  In 
three  dimensions,  it  does  not  necessai  Jy  follow  that  there  is  a  real 
energy  gap  in  the  density  of  states,  and  these  materials  in  principle 
could  exhibit  semimetallic  rather  than  sem  conducting  behavior,  ft 
ir  clear  that  when  overlap  is  small  compared  to  the  exchange  energy. 


the  bands  will  be  narrow  and  a  real  gap  will  occur.  But  this  is  not  a 
necessary  condition  for  the  existence  of  <•  gap.  Consider  a  non¬ 
magnetic  crystal  with  no  exchange  potential,  and  assume  that  E(k) 
has  been  calculated  for  the  fdouble-sized)  first  Brillouin  eone.  When 
the  exchange  potential  is  turned  on,  this  zone  is  split  into  tw-o.  If 
the  original  energy  levels  in  the  outer  half  of  the  large  zone  were 
higher  than  those  in  the  inner  half  c£  the  zone,  then  a  gap  exists  even 
for  a  relatively  small  exchange  energy.  Such  a  situation  could  occur 
for  the  asymmetric  3d  electrons  in  a  low  symmetry  lattice,  as  well 
as  in  special  highly  symmetrical  lattices,  such  as  body  centered  cubic. 
We  shall  discuss  this  furtli-;r  in  Appendix  B,  where  attempts  at  energy 
band  calculations  will  be  summarized. 

The  L.  C.A.  O.  calculation  proceeds  au  in  the  one-dimensional 

case.  As  before,  we  take  one  3d  electron  per  cation.  Let  P  be  the 

a 

positions  of  the  ions  whose  :\d  electrons  have  spin  up  at  T  =  0,  and 
let  Rp  be  the  positions  where  the  3d  ele  .rons  have  spin  down.  Then 
the  Bloch  wave  functions  are: 


^  iK  ’  R 

4*i(^,r)  =  c..  ®  p(f  -  R^) 
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(2.28) 


Then: 
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taking  an  ion  of  type  n  at  the  or*gin.  Let  the  origin  be  the  Site  ui  a 
spin  up  ion,  and  let  f  =  p.  be  the  positions  of  the  nearest-neighbor 
spin  down  ions.  Then,  in  the  nearest -ne  ghbor  approximation; 

H,,(ki  =  pOB  k  ■  P' 

Hj,(k)  =  H^jCk)  =  Pt 


where  p.'  are  the  positions  of  the  nearest  “neighbor  sp-tn  up  ions,  and 

c  =  /  dFuCT)  H  (“r)p(r 

o  c  c 

“o  =  /  drp(r  -~)  -7) 

c  ,  =  /  dr  p  (  r  I  H  =  r '  p  (  r  P  ) 

c  j  =  /  d?p(T  -  p  i  t  -  P  -  R  ) 

Pj  /  dFpCT  -  F)  *13.,  ’  ^  e  -  -  p  -  p‘)  . 

We  do  not  know  where  ^n  k  =  spaf  »^  the  top  ot  the  v'alenct  band 
and  the  bottom  of  the  conduction  band  are.  From  our  previous  dis¬ 
cussion  It  is  reasonable  to  assume  tha’  ^hesi*  points  art-  both  somt  - 
where  on  the  planes  k  •  "p  -  ±ir/Z.  However,  this  is  not  necessat  Jy 
the  {.ase.  For  simplicity,  we  shall  assume  the  simplest  situa*;cn: 
the  conduction  band  minimum  and  the  valence  band  maximum  both 


at  k  =  0,  Since  we  are  using  this  calculatien  primarily  to  evaiaate 
the  average  change  in  spin  as  electrons  are  excited  across  the  gap= 
and  since  only  the  zeroth  order  change  will  turn  out  to  be  significant 
in  the  evaluation  of  p,  no  error  is  introduced  at  this  point.  The 
calculation  has  been  repeated  using  more  physical  assumptions  about 
the  band  extrema  [conduction  band  minimum  at  k  =  ^  (110),  valence 
band  maximum  at  k  -  ^  (HI)] ,  and  it  was  found  that  only  the  higher 
order  term#  are  affected.  Near  k  =  0,  in  the  (111)  direction,  the 
secular  equation  (2.29)  can  be  written: 


.  8 , 2  2.  _ 
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(2.30) 


pj(6  -|kV) 


C:  j.  ^  tL  8,22.  _ 


As  before,  we  consider  an  electron  with  spin  down  and  assume  that  the 
bands  are  narrow  compared  to  the  gap.  Then  the  solution  of  (2.  30)  is: 
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Equations  (2.31)  show  that  the  assumptions  about  the  band  extrema  are 
self-consistent  only  if  e  ^  is  negative  and  c  j  is  positive.  The  band  gap 
is  then: 


The  corresponding  wave  functions  are: 
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Equations  (2.  32)  lead  to: 
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(2,  34 


Using  the  same  assumption  about  the  form  of  E  as  in  the 

g 

dimensional  calculation: 


n  -T 

E  =  (n- _ 


(2.  35 


where  V  is  the  total  volume  of  the  crystal.  From  {2.3-4)  and  (2.  35L  wt 


obtain  the  reeultj 


(2,36) 

Ikjuation  {2.  36)  ahowt  that  P  in  (2.  1)  is  a  constant,  independent  of  n, 
to  within  terms  of  the  order  of  the  square  of  the  small  quantity, 

"  *o^'  coefficients  of  these  terms  should  not  be 

disturbing,  since  in  Appendix  B  we  show  that  the  bandwidth  is  greater 

than  12  ,  while  the  band  gap  is  considerably  smaller  than  (a  -'a  ) 

oo 

Thus,  the  value  ol  72  Pj  /  is  a  good  deal  less  than  half  the 

sq*i4re  of  the  ratio  of  bandwidth  to  band  gap,  which  we  expect  to  be 
itself  small  in  the  materials  under  consideration  7'hereiore,  equation 
(2  i)  will  remain  valid  over  a  relatively  large  range  of  n  ,  7,T\e 

narrower  the  bands  are,  the  larger  the  range  of  n  for  which  (2.  1) 
holds.  From  (2.  36),  the  value  of  p  is: 

P  =  ^=go/N'  (2.37 

Comparing  (2.  3?)  \vith  the  result  of  the  one-dimensional  calculation, 

(2.  2i),  shows  that  the  vaiu>»  of  ^  ic  the  same.  Note  that  in  the  limit 
of  narrow  bands,  a,  -  Uj  =  Pj  -  0,  and  (2,  37)  is  ■’xact  within  the 
validity  of  (2.  35), 


3.  Virtual  Crystal  Approximation 

We  shall  recalculate  &  using  an  entirely  different  model  but 
maintaining  the  spirit  of  the  Kartree-Fock  approximation.  We  con¬ 
sider  metallic  ions  with  spins  ordered  antiferromagnetically 
at  T  =0.  L;et  an  ion  with  spin  down  be  called  type  A,  a  spin-up  ion 
type  B.  Consider  the  sublattice  where  all  spins  are  down  when  perfect 
order  exists.  In  the  vicinity  of  an  A  ion.  an  electron  w'lth  spin  up  sees 
a  potential  vt .  whereas  a  spin-down  electron  sees  V*  S'milarlv 
Vg  and  Vg  are  the  potentials  in  the  vicinity  of  a  B  ion  seen  by  a  sp.,n-up 
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and  a  spin-down  electron  respectively,  ’t  is  clear  that  and 

V'a  ^  For  simplicity,  we  shall  take  the  case  of  one  3d  eiectron 

per  cation.  Since  N  is  the  total  density  of  cations,  the  density  on 
each  snblattics  is  N/Z. 

The  average  potential  seen  by  an  electron  with  spin  up  on  ihe 
sublattice  under  consideration  is: 

where  n,  is  the  number  of  A  ions  on  the  suolattice,  n*,  is  the  number  of 
A  B 

B  ions.  Similarly,  for  a  spin-down  electron: 


z,  J9l 


The  average  exchange  energy  is  then  just  the  difference  between  (2.  38i 
and  (2.  39),  or: 


<''ex>4<"A-  "B*  ^''a 


(Z.40i 


For  perfect  order,  n*  -  N/Z,  n_  ~  0,  and' 


^^ex/T=0  ~  ~  ' 


For  complete  disorder,  n,  -  n„  -  N/4,  and: 


/V  >  =  0. 

\  ex  / 

Since  n  +  n  =  N/Z,  and  the  number  of  intrinsic  carriers  is  just  the 
A  B 

number  of  ions  with  spin  up  on  the  sublattice,  (2.40)  can  be  written. 
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ex  /t 


4ti/m 


If  w€  as#imie  that  the  energy  gap  ia  proportional  to  the  average  exchange 
energy,  a*  was  found  In  (2.  18).  then  (2.41)  becomes: 

C2. 

Thus  p  ”  4E^^/N.  exactly  the  same  as  found  in  equations  (2.2?)  and 
(2.  37)  for  one  34  electron  per  cation.  2t  is  clear  that  if  there  are  z  Id 
electrons  per  cation  filling  the  valence  band,  the  result  for  8  in  the 
above  approximations  would  be: 

p  -  4  E  /  zN , 
go' 

D*  Crystalline  Structure  Distortion 

The  existence  of  antiferromagnetism  is  not  a  necessary  condition 
for  the  applicability  of  (2.  1).  The  relation  can  also  he  shown  to  he 
appropriate  when  an  energy  gap  is  caused  by  a  crystal  structure  dis¬ 
tortion  to  lower  symmetry.  This  type  of  gap  can  arise  from  an  energy 
‘hie  tO  heroieal  binding  —  the  lower  band  may  be  thought  of  as  a 
ha  the  upper  an  antihonding  band.  Excitation  of  an  electron 
across  the  t  rgy  gap  decreases  the  gap  because  the  excited  electron 
no  longer  contributes  to  the  chemical  binding*  Thus  the  situation  Is 
analogous  to  the  case  of  an  antiferromagnetic  crystal  dealt  with  in 
section  C.  However,  the  evaluation  of  0  is  more  difficult  when  the 
gap  is  caused  by  crystalline  distortion.  In  this  section,  we  shall  cal¬ 
culate  expressions  for  P,  using  simple  one -dimensional  models,  but 
employing  as  much  as  possible  the  physical  properties  of  the  vanadium 
oxideo  to  which  we  expect  the  theory  to  apply. 


very  •trongly  screened*  which  is  not  too  far  from  the  case  where 

an  extremely  high  dentity  of  free  electrons  exists.  When  the  a  .  are 

J 

independent  of  j,  VCx)  is  pericdic  and  we  can  solve  Schrdd  nger’s 
equation  exactly;  when  the  a  ^  depend  on  j,  Bloch  states  are  no  longer 
eigenstates,  and  the  problem  cannot  be  solved  exactly.  Consequently, 

let  us  assume  0^=0  for  all  j.  Then  Schrddinger's  equation  may  be 
written; 

jj2 

'^o  f  (2.43) 

wher  J  bs  (1  -  2€o)  a/2. 

Thus; 


=  A  sin  Y  -f  B  cos  y  (0  .'•'x  <  b) 

.^(x)  =  A«  siny^  +  B'  cos  Y^  (b  «  x  <  a) 

2  2 

where  y  s  2  m  E/fi  .  Bloch’s  theorem  gives: 

>!i{x)  =  e^^  p(kx) 

where: 

p<kx)  =  p[k{x  f  k.a)]  . 

Then: 

p{0) 

^(b)  =  e^^^p(kb). 

But  from  (2.44)  and  (2.46): 
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4*(0)  =  ^(0)  =  B 


4>{h)  -  A  sin  yh  +  B  cos  \  b  =  A’  sm  \  b  +  B'  cos 
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Integrating  (2.43)  from  0“  to  0^; 


4»*(0'*‘)  -  ij;'(0")  =  2k  ij/;0)  =  2  XB 
2 

where  X  --  m  V^/fi  .  Similarly,  integrating  (2.43)  from  b‘  to 

iji  (b  )  -  4<’{b  )  =  2Xii{b)  -  2X[A  sin  yb  t  B  cos  y  b]  .  (2.  49; 

From  (2.44),  we  see: 


vj<'(x)  -  Ay  cos  yx  -By  sin  yx  (0  <  b) 

(|i’(x)  =  A'  y  CQS  yx  -  B’  y  sin  yx  (b«  x  <  a) 

while  from  (2.  45): 

4>'(x)  =  e  4'*(x  4  a)  , 

Gomb.ning  (2.48),  (2.49),  (2.  50),  and  (2.51): 

y  A  -  e  [A'y  cos  ya  -  B‘y  sin  ya]  -  2  X  B 


(2.  50i 
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5*; 


A’y  cosyb  -  B*y  sin  yb  -  Ay  cusyb  +  By  slnyb 
=  2X[A  sm  yb  ~  B  cos  ybj  . 

The  secular  equation  is  obtained  from  (2.47)  and  (2.  52): 


sin  yb 

cos  y  b 

Sin  yb  -  cos  y  b 

0 

Jka 

w 

e 

sin  y  a  -  cos  y  a 

ika 

y  c 

Tka 

-ZX  e 

y  cos  y  a  y  sm  ya 

2X  sin  yb  +  y  cos  yb 

2X  cos  yb  -  y  sin  yb 

y  cos  yb  y  Sinyb 

. .  2  tka  .  ,  *  y 

-4X  e  sinvbsiny(a 

-b)+4y  Xe^^^cos  y  a+2y  ^ 

e  cosya-y  (i+c  ‘  i 

.  53 
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SoHdng  {Z,  53)  for  cos  ka; 


coska  s  cosva  +  ^  «in\b  sinv(a  -b) . 

Fo.  ,  teatcr  clarity,  let  K  =  ka.  y  =ya.  z  =  Xa,  recall 
®  ^  =  (a  -  2b)/2a.  Then  (2.  54)  becomes; 


C2.54) 


I  1 

2  siny  tc)  siny  ( ^  4  c  «  ) 


CO*  K  =  cos  y  4  Zz  --ip .y  -4  2z 

y 


(2.55) 

Equation  (2.  55)  give*  the  energy  band  structure  of  the  crystal.  We 
still  have  two  parameters  at  our  discretion,  the  strength  of  the  inter¬ 
action,  z,  and  the  amount  of  distortion,  €.  Since  we  wish  to  apply 
the  theory  to  narrow  band  materials,  we  must  choose  z  relatively 

large.  A  good  choice  for  €  is  c  ^  0.01,  the  distortion  present  in  V  O  . 

2  3 

The  energy  band  structure  at  T  =  0  (o  =  1)  is  plotted  in  Fig.  11- 1,  for 
z  =  -60,  c  =  0,009.  As  can  be  seen,  the  bands  are  narrow.  In  this 
case,  the  relevant  energy  gap  is  narrow  also,  as  we  might  expect 
^2^3*  ®  decreases,  the  gap  shrinks.  In  Fig.  II-2,  the  band 

structure  is  plotted  for  g  =  2/3;  Fig.  n-3  shows  E(k)  for  g  -  1/3. 


Clearly,  the  gap  must  vanish  as  c  goes  to  zero. 

It  is  important  to  demonstrate  that  c  decreases  with  increasing 
con»  entration.  In  the  distorted  state,  the  wave  functions  of 
electrons  in  the  valence  (bonding)  band  have  larger  amplitudes  in 
the  region  where  the  cations  are  closely  spaced,  whereas  those  of 
electrons  in  the  conduction  (antibonding)  band  have  larger  amplitudes 
in  the  region  of  wide  spacing.  Thus,  excitation  of  an  electron  across 
the  gap  removes  charge  d^jnsity  from  the  closely  spaced  region. 
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Consider  a  hydrogen  molecule.  With  both  electrons  in  bonding  orbi¬ 
tals,  the  equilibrium  distance  between  the  hydrogen  ions  is  0,  74  A, 
when  one  electron  is  removed,  we  are  left  with  the  hydrogen  molecular 
ion,  for  which  the  equilibrium  distance  is  1  06  A,  The  ions  thus 

move  apart  as  an  electron  is  removed  from  a  bonding  orbital.  Of 
course,  the  hydrogen  molecule  is  made  up  of  ions  and  electrons  which 
interact  by  means  of  a  Coulomb  potential,  rather  than  the  delta - 
function  potential  with  which  we  are  dealing.  Conceivably  cnis  could 
make  an  important  difference.  Consequently,  the  hydrogen  molecule 
analogy  is  analyzed  in  Appendix  C  using  delta-function  interactions.  It 
is  found  that  the  equilibrium  distance  of  is  0.  76  of  that  for 

which  is  roughly  the  same  as  the  0.70  ratio  for  the  real  molecules. 

We  conclude  that  a  decreases  with  increasing  concentration  of  free 
carriers.  It  is  clear  that  c  vanishes  after  N/2  carriers  are  excited, 
since  the  charge  density  between  all  the  cations  is  the  same,  and  thur 
the  cations  must  be  equally  spaced.  Thus  a  decreases  from  1  to  0  as 
n  increases  from  0  to  N/2.  We  shall  assume  the  simplest  possible 
behavior  for  a  ,  the  linear  relationship: 

a  =  I  -  2n/N  .  (2.  56) 

The  resul  of  an  involved  calculation,  to  first  order  in  n/N,  is: 

=  3.3  {J  -  3.  1  n/N). 

The  value  of  p  is  thus; 

^  =  3,1  ,  (2.  57) 


It  is  worthwhile  also  to  calculate  p  when  the  bands  are  narrow 


. . . . . . . . 


Compared  ts  the  energy  gap  For  this  purpose,  we  investigated  the 

case  ifehere  2  -  -6,  c  -  0.  i  b.  1  he  band  structure  is  shown  in 

Jigs.  11-4  through  11-7  for  q  -  i,  Z/3 ,  1/3,  and  0  ,  It  can  be 

seen  how  the  energy  gap  shtinks  to  zero  with  decreasing  a  .  while 

the  valence  and  conduction  bands  become  wider  and  move  towa-'  *s 

each  other  The  result  for  the  variation  si  £  is  * 

g 

E  =  15.  Z  (i  ^  3.  7  n/N)  , 


Hence: 


p  ^  3  7  E  /N. 

go^  {I 

In  general,  p  was  found  to  be  insensitive  to  changes  in  the  bandwidth 

for  a  constant  band  gap;  however,  p  was  sensitive  to  changes  in  the 

gap,  varying  from  1,6  E  /N  for  narrow  gaps  to  4,  9  E  M  for  the 

^  go' 

widest  gaps  investigated. 


Figures  11*4  through  II-7  can  give  us  some  insight  into  the 
nature  of  the  distortion,  f  igure  II-7  shows  the  band  structure  for 
the  undistorted  situation,  where  there  is  no  gap,  but  rather  a  half- 
filled  band  The  distortion  can  be  looked  at  as  the  band  generalization 
of  t.he  Jahn-Teller  effect.  When  the  bands  are  wide,  a  Jahn-Teiler 
distortion  would  just  introduce  a  small  gap  at  the  ends  01  the  reduced 
first  Brillouin  zone,  resulting  in  only  a  small  gain  in  electronic 
energy  This  gam  w'ould  quickly  be  overcome  by  the  loss  in  strain 
energy,  and  the  distortion  most  likely  would  be  sufficiently  small  so 
that  a  real  gap  in  three  dimensions  would  not  occur  Howe'^^er  -- 
narrow  band  case,  the  distortion  has  an  effect  on  all  points  in 
k-space,  lowering  the  energy  of  tne  entire  valence 
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COS  -  X  =  S  cos  X 

a  o  o  a  o 

t-  =  o'  E  . 

IT  ^ 

Then  {Z.61)  can  be  written; 

']i'*(z)  +  [C  +  ^  s  cos  2z  ]  44z)  =  0  .  {2.  62) 

Equation  (2.  62)  is  just  a  form  of  Mathieu’s  equation,  the  eigenvalues 
of  which  have  been  tabulated  [34]  .  Take  =  100,  and  for  simplicity 
assume  that  sufficient  electrons  are  present  to  fill  the  lowest  two 
bands,  leaving  all  higher  bands  completely  empty.  Set  £  -  0.01, 
the  correct  value  for  V^O  ,  A  long  cal'-ulation  shows; 

E  =  0.  30[  1  ^2^1  . 

g  t  N  J 

Hence; 

p  =  3.  5  E  /N.  (2.63) 

go 

For  a  Mathieu  interaction,  the  value  of  p  depends  som.ewhat 
on  the  number  of  electrons  present  and  strongly  on  the  amount  of 
distortion  e.  It  was  found  that  p  increases  for  decreasing  c  .  In 
many  ways,  the  Mathieu  interaction  is  inferior  to  the  delta -lunction 
interaction  for  our  purposes.  The  Mathieu  potential  is  a  very  weak 
attraction,  and  at  high  temperature  actually  vanishes,  leaving  the 
free  electron  Schrodinger ’ s  equation.  This  results  in  a  band  struc 
ture  where  the  bands  are  wide  and  the  gap  is  narrow,  just  the  opposite 
of  the  physical  situation  to  which  we  want  to  apply  these  results.  At 


'X 
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o 
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low  temperatures,  this  could  be  remedied  by  using  extroiuely  large 
values  for  s^;  howevei  ,  no  tables  of  Mathieu  functions  exist  for 
s  >  100,  so  such  calculations  are  no»  practical  at  present. 


Chapter  III 


CALCULATION  OF  CONDUCTIVITY  AS  A 
FUNCTION  OF  TEMPERATURE 


The  electrical  conductivity  of  a 
=  n  e  Vi 


st'miconductor  can  be 


expressed 


(3.  1) 


where  n  is  the  concentration  of  excited  carriers,  e  is  the  electronic 
charge,  and  p  is  the  mobility  of  the  carriers.  In  this  chapter,  we 
shall  assume  that  Eq.  (2.  1)  is  applicable,  and  calculate  n  as  a  function 
of  temperature.  We  shall  evaluate  n(T)  in  tw-o  opposing  limits,  the 
effective  mass  approximation  and  the  limit  of  narrow  bands. 

Let  us  take  the  bottom  of  the  conduction  band  as  the  zero  point 
of  energy.  The  probability  that  a  state  with  energy  E  will  be  occupied 
by  an  electron  is  given  by  the  Fermi  function: 


""  E-E, 


(3.2) 


kT 


4  1 


Cleaily,  the  probability  that  a  state  is  occupied  by  a  hole  is: 


f^(E)  .  1  -  f^(E) 


1 


(3. 


3) 


Let  p  (E)  be  the  numbe'  of  states  per  unit  volume  between  energies  E 
and  E  +  dE.  Then,  i  .ling  (3.2),  the  concentration  of  electrons  in  thr- 
conduction  band,  n,  is  given  by. 
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n 


dE  pJE) 


(3.4) 


Similarly,  defining  ^^{E)  as  the  density  of  states  in  the  valence  band, 
and  using  (3.  3): 


-E 


g 


aE  p^(E) 


'-00 


(3.5) 


where  p  is  the  concentration  of  holes  in  the  valence  band.  For  an 
intrinsic  semiconductor,  the  condition  of  neutrality  requires  n  =  p. 
This  condition  determines  the  Fermi  energy,  E_. 

if 

A.  Effective  Mass  Approximation 


We  first  consider  the  situation  where  the  valence  and 

conduction  band  structure  can  be  well-represented  by  ellipsoidal 

constant-energy  surfaces.  We  can  then  define  densitv-of-states 

effective  masses,  m  for  electrons,  m,  for  holes: 

’  e  ’  h 

m  =  M{mj  m^ 

where  m^,  m^.  m^  are  the  principal  values  of  the  effective  mass 
tensor,  and  M  is  the  number  of  equivalent  extrema.  The  densities 
of  states  are  now  given  by: 


j  2m  3/2  j 
p(E)=-L(— E*/^ 

2‘if^ 


(3.6) 


1  ^"^h  1/2 

p  (E)  ^  )  {-E  - 

''  2ir  ® 


(3.7) 


. . . . . . . . 
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ii.quations  (3.  6)  and  (3.  7)  are  always  true  near  the  band  extrema. 
However,  for  large  values  of  and  they  probably  do  not  re¬ 
main  valid  for  a  very  large  range  of  energy.  The  effective  mass 
approximation  assumes  that  these  relations  remain  applicable  for 
all  values  of  E  which  have  finite  probability  of  occupation  at  the 
temperature  under  consideration.  Thus,  the  effective  mass  approxi¬ 
mation  is  a  good  one  whenever  |Ep| »  kT,  and  ni  is  not  too  large. 

In  view  of  (3.6)  and  (3.7),  (3.4)  and  (3.5)  yield; 


,  2m  3/2  r 
1  ,  e  . 

2u 


Ji 


00 


dE 


E-E^ 

TT 


,  2m,  3/2 


dE 


-iQC- 


Ejr-E 


+  1 

1/2 


kT 


+  1 


(3.8) 


(3.9) 


We  first  consider  the  case  where  jE  -  E  |  »  kT,  and  Boltzmann 

r 

statistics  apply.  Then,  the  exponentials  in  the  denominators  are 
much  greater  than  unity,  and  (3.8)  and  (3,9)  become: 


n  = 


( 


in 


,  1/2  -X 

dx  X  e 


,  2m  kT  3/2 
1  .  e  . 


_  ^  /  e  ,  '  kT 

'4  ^ ^  ^ 
trfi 


(3.  10) 


Ep+E 

Zm^kTl/Z 


*■  i  \ 

P  =  4  *  ■  ,2"  > 


(3.11) 
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the  condition  of  neutrality  gives: 


Ej,  Ep  E 

3/2  FT  3/2  “FT  "F 
m  e  -  m,  e  e 
e  u 


IT  ^  x  3kT  ,  "'h 


(3. 12) 


Defining  m  as  the  geometric  mean  of  m^  and  m^,  (3.  10)  yields: 


_  1  .Zm^kT.S/Z  'zA 
^  J  i - T"  ^  e 


Letting  A  s  (2m  kX/wfi^)^/^  / 4,  (3.  13)  becomes: 


n  =  e 


£ 

'zA 


(3. 13) 


(3.  14) 


Substituting  (2.  1)  into  (3.  14): 


.„3/2  '2#r 

n  =  A  T  e  e 

Elqiiation  (3.  15)  can  be  solved  explicitly.  Let; 


(3.15) 


T  =  exp 


BAT 


(3.  16) 


i  E 

I 

n.  1  1  2kT 


(3.17) 


Substituting  (3.  16)  and  (3.  17)  into  (3.  15),  we  are  left  with  the  simple- 
equation: 


/Z  =  T  . 

The  solution  of  (3.  18)  can  be  written: 


(3.18) 


50 
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1) 


0 


') 


) 

) 


(3.  19) 


Equation  (3.  19)  giverf  n  as  a  function  of  temperature.  The  function 
■n(r),  which  can  be  called  the  infinitely  repeating  exponential,  is  a 
function  which  will  recur  often  during  this  analysis.  I*  is  not  dis¬ 
cussed  anywhere  in  applied  nnathematics  literature,  and  has  some 
unusual  properties  which  are  worth  detailing  at  this  time.  Its  be¬ 
havior  is  shown  in  Fig.  III-l.  i)  slowly  rises  from  0  to  e  as  x 

1 /e  1 /e 

increases  from  0  to  e  =  1.445.  However,  at  the  point  x  =  e  '  ' 

an  essential  singularity  exists  and  q  increases  without  limit.  Not 

only  does  r\  blow  up  at  this  point,  but  all  its  derivatives  are  infinite 

also.  Such  a  singularity  in  ii(T)  conjpletely  eliminates  the  band  gap, 

and  brings  about  a  first-order  semiconductor-to-rnetal  transition.  We 

can  determine  the  transition  temperature  easily.  At  the  point  of  singularity 
i  /e 

T  ~  e  '  ,  and  equat  on  (3.  16)  yields: 


M  t^/2  - 

2k  o 


(3.  20) 


Let: 


{ 


2k 

e  pA 


R=  e 


) 


Th€*n  (3.  20)  becomes; 
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R 

R 

Q  ^ 

Q  =  =  R  r  ti(R)  .  (3.21) 

JEiquation  (3,21)  gives  the  transition  temperature,  1’^,  in 
terms  of  the  band  parameters,  ^  A  transition  always 

occurs  in  this  case,  since  R  <  1,  and  thus  (3. 21)  always  has  a  solu¬ 
tion.  Using  Boltzmann  statistics,  the  number  of  carriers  grows 
without  limit  at  the  transition  temperature.  We  shall  be  interested 

in  the  ratio  E  /kT  ,  which  can  be  expressed: 
go '  o 

E  ^  ,  , 

=  Q  In  (^)  =  ri(R)  in(^)  .  (3.22) 

o 


In  Chapter  II  [see  equations  (2.27),  (2.37),  and  (2.42)],  we 
fr  und  that  for  an  energy  gap  arising  from  antiferromagnetism; 

j3  =  4  E  /N  .  (3.23) 

go' 

In  equations  (2.57),  (2.  58),  and  (2.63),  we  saw  that  (3.23)  is  also 
a  good  approximation  when  the  gap  is  due  to  a  crystalline  structure 
distortion.  Let  us  write’ 


p  ^  4(1  +  6)  E  /N 
go' 


(3.24) 


where  6  vanishes  in  the  c?.se  of  antiferromagnetism,  and  is  small, 
probably  between  -1/4  and  l/4,  in  the  case  of  crystalline  distortion. 
Then  the  parameter  In  (1/R),  which  determines  E  /kT  ,  becomes: 


In  (ij  =  z(i  ^  a:/  [|;  ( 


m 


wh 


E  3/2 

-^)  ] 


C  25) 


The  quantity  in  brackets  in  (3.25)  is  essentially  a  measure  of  the 
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r-irio  of  the  band  gap  to  the  band  width,  [in  fact,  using  (2,  18),  from 

the  one-dimensional  theory  of  antiferromagnetism:  ln{^)  s 

K. 

(1  +  6)  (  )  .  ]  E  /kT  is  plotted  as  a  function  of  ln(l/R)  in 

4TT^  Efc  J  go^  o 

Fig.  XII-2.  For  narrow  bands,  ln(l/R)  is  large,  and  E  /kT  is  a 

o 

slowly-varying  function  of  R.  Boltzmann  statistics  are  valid  for 

roughly  >  2. 

The  important  point  is  that  a  semiconductor-to-metal  transi¬ 
tion  can  occur  in  a  simple  way  in  a  band  theory.  The  resulting  be¬ 
havior  of  In  r  as  a  function  of  l/T  is  shown  in  Fig.  III-3  for  reason- 
able  values  of  E  and  m  . 

go 


Just  before  the  transition,  e.  Thus,  from  (3.  17),  at  this 


point' 


n(T  )  =  e  A  T^/^  e 
o  o 


(3.  26) 


But  it  can  be  seen  from  (3.  15)  that  (3.26)  implies: 


n(T^)  =  ZkT^P. 


(3.27) 


Substituting  (3.24)  into  (3.27): 


n(T  ) 
o 


2(1  +  6) 


(^)  . 


(3.28) 


From  Fig.  iXI-2,  we  see  that  for  the  narrow  band  situation  with 

which  we  are  nost  concerned,  E  /kT  >  5.  Equation  (3.  28)  then 

go  o 

shows  that  the  maximxim  value  of  n/N  before  the  transition  is: 


(&)  <  J 

'  N  max 


(3.29) 
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Equation  (3,29)  is  a  very  important  result,  since  it  shows  that  our 
use  of  (2.  1)  for  all  n  up  to  the  transition  was  valid.  Furthernnore, 
if  n/N  had  approached  unity  before  the  transition,  the  value  of  the 
effective  mass  approximation  in  arriving  at  (3.8)  and  (3.9)  would 
have  been  brought  into  question.  If  E^^/kT^  is  small,  and  (3.  29) 
does  not  hold,  Boltzmann  statistics  do  not  apply  and  Fermi  statistics 
must  be  used. 


where: 


In  an  approximation  valid  up  to  y  =  1-2,  or  just  before  the 
system  becomes  completely  degenerate,  Ehrenberg  [35]  has  shown 
that  ¥ ,  ^(y)  can  be  written: 

r,/2(y)  =  — (3.32) 

'  1  +  4  e’> 

Th.  s  approximation  will  be  valid  until  the  Fermi  energy  is  well  into 
the  conduction  band,  so  that  we  can  think  of  (3.32)  as  exact  in  all 
normal  s  tuations  involving  semiconductors.  Then  the  condition  of 
neutrality  gives,  from  (3.30),  (3.  31),  and  (3.  32): 


I 


.„+E  J  A  » 
I  +  4  e  ^ 


#  '»  The  Aolution  of  (3.33)  i«: 

E„7i^ 


(3.33) 


/ 


M  .  1  . . 

''T“  *  J 


.  M  -  1  .2  ^  ,  g' 

}  +  M  e  * 


(3.34) 


Ci^^ess TAtioti  of  electrons  in  the  conduction  band  ic  then: 


2m  kT  3/2 
N  =  f - -  ) 


1  +  4  e 


(3.  35) 


case  m  -  m,  =  m  ,  (3.  35)  becomeb; 
e  ft 


E 

/  --S- 

-  -.3/2  kT 
n  =  A  T  '  e 


1 

;  -E  /2kT 
1.1  K 
1  +  -T  e  * 

4 


(3.36) 


Th..*  differs  from  the  Boltzmann  expression,  (3.  14),  only  by  a  factor 
of  [I  e*p(“E  /2kT)/4]  which  is  very  nearly  unity  for  all  ordinary 
temperatures  unless  the  energy  gap  is  extremely  small. 

Combining  equations  (2.  1)  and  (3.36): 


3/2  "5^r 

'A  A 


n  =  AT 


1 


go  g 

~2PT  IkT 


(3.37) 


e 


Equation  (3.37)  can  be  solved  analytically  only  when  the  initial  energy 
gap  is  much  larger  than  kT  for  all  temperatures  under  consideration. 
In  that  case,  let: 

£ 


T  £  [ c( 1  - 


kT  2kT  TFT, ,  IFT  ^  ® 


)] 


(3.  38) 
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Equation  (3,37)  then  becomes; 


7l  =  V^  -  T  =??(t)  (3. 39) 

where  is  defined  as  in  equation  (3.  17).  Thus  a  transition  still  occurs 
with  Fermi  statistics. 

Just  before  the  transition,  /t  =  e.  Thus  the  concentration  of 

carriers  at  the  point  of  transition  is  the  same  as  the  Boltzmann  result 

in  terms  of  T  : 

o 


n(T^)  =  2kT^;p. 


(3.27) 


But  for  Fermi  statistics,  a  different  value  for  T  is  obtained.  The 

o 

-  l/e 

transition  temperature  occurs  when  T=e  ,  which  yields  the  condition: 


,1/2 


i  +  In  [ 1  - 


epAT*'"  -E  /2kT 
o  go'  o 

k  2k~  ^ 


2PAT  '  e 
o 


E  /2kT 

1  -  2k  _  go'  o 
o 


(3.40) 

Equation  (3.40)  can  be  solved  graphically  for  T^.  However,  it  can  be 

seen  that  Fermi  statistics  tend  to  raise  the  transition  temperature. 

Since  exp(-E  /ZkT  )  is  a  small  parameter,  an  expansion  of  (3.40)  to 
go  o 

first  order  gives: 


go 


e6AT^^^  ZkT 


~Zk 


=  e 


°  ^  I- 


(3.41) 


Ignoring  the  relatively  small  second  term  on  the  right  of  (3.41),  we 
recapture  the  Boltzmann  relation  for  T^,  equation  (3.20).  The  presence 
of  the  small  positive  addition  to  the  right-hand  side  increases  T^  by  a 
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smaii  amount.  Equation  (3.Z7)  then  shows  that  slightly  more  carriers 
have  been  excited  at  the  point  of  transition  than  had  been  when  Boltzmann 
statistics  were  used. 

In  the  case  where  the  electron  and  hole  effective  masses  differ, 
M#  1,  and  we  must  modify  equations  (3.38)  and  (3.41).  These  become: 


, _ ^  Jn 

„  _  .  «3/2  2Fr  2 
N  -  A  T  '  e  e 


go  7 

,^1-M  '2^  ri?f^(M-l)  "k 

‘  +  wm '  =  *  t^sm-  ' 


(3.42) 


and: 


epAr 


o  ^  e  (3  -  M) 

8  Jm 


(3.43) 


Equation  (3.43)  shows  that  the  transition  temperature  can  be  a 
sensitive  function  of  the  mass  asymmelry,  M.  When  the  electron 
effective  mass  is  about  twice  the  hole  effective  mass,  is  approxi¬ 
mately  the  same  as  for  the  Boltzmann  case.  For  greater  ratios  of 
m^  to  m^,  the  transition  temperature  is  lowered,  but  this  lowering  is 

small  for  normal  values  of  M,  As  an  example,  when  E  is  6kT,  and 

go 

the  electron  effective  mass  is  100  times  the  hole  effective  mass,  T 

o 

will  be  decreased  by  roughly  15%.  On  the  other  hand,  a  large  effective 

mass  for  holes  relative  to  that  for  electrons  tends  to  increase  T  some- 

o 

what  more  sharply.  When  E  is  6kT  .  and  m,  is  100  m  ,  T  in- 
^  ^  go  o  h  e  o 

c, teases  by  about  50%.  It  is  important  to  bear  in  mind  that  (3.43)  is 


valid  only  for  values  of  E^^  greater  than  approximately  4kT^. 

When  E  is  not  this  large  relative  to  kT  ,  we  must  return  to 
go  o 


1 


equation  (3.38),  We  cannot  solve  (3,38)  analytically,  but  we  can  solve 
for  the  transition  temperature.  Using  (3,23): 


,1 

j  ,  i  +  -T  e  ^ 

d  in  n  _  4 

dT"  ^  2n 

j _ kT  N 

1  +  n-  e  ® 


(3.44) 


From  (3.44),  we  see  that  the  transition  occurs  when: 


E 


_ £  (1  .  .in  ) 

o  2n  ,  ^  1  2kT  -  N  ' 
kT  TT  "  ^  ^  4  ^ 


(3.45) 


and’ 


Tf  we  let: 


-E  /2kT 
e  go'  o 

4*^ 


o 


then  equation  (3.45)  becomes: 


U  -  -  B  ®  ^71{B)  .  (3.46) 

Equation  (3.  46)  can  be  used  to  find  when  the  energy  gap  is  small. 

A  major  effect  of  Fermi  statistics  is  the  cutting  off  of  the  con¬ 
ductivity  anomaly  at  very  high  n.  This  can  easily  be  seen  from  (3.  35). 

The  maximum  v'alue  for  n  at  the  transition  temperature,  after  the 
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transition,  is: 


2m.  kT  3/2 

n={ - (3.47) 

whereas  with  Boltzmann  statistics,  n  grew  without  limit  after  the 
transition.  Such  a  cut-off  is  characteristic  of  Fermi  statistics. 

l/sing  (3.27)  and  (3.47),  we  obtain  the  jump  in  carrier  concentra¬ 
tion  at  the  transition  temperature: 


An  =  ( 


2m,_kT 
h  o 


fi' 


3/2 

) 


NkT 

o 

■  • 
go 


This  result,  as  well  as  (3.47),  is  physically  meaningful  only  if  the  value 

of  n  given  by  (3.47)  is  significantly  less  than  N,  since  the  effective 

mass  approxin -vtion  must  remain  valid  in  order  to  use  (3.35).  If  n  in 

(3.47)  is  greater  than  N,  cleaily  there  can  be  no  more  than  N  carriers 

per  unit  volume  in  the  conduction  band,  and  the  jump  in  n  at  T  is 

o 

lurji: 

An  :=  N(i  -  RT  /2E  )  , 

o'  go 

This  teiis  us  nothing  about  the  jump  in  conductivity,  since  the  mobility 
will  change  considerably  upon  transition  from  a  semiconducting  to  a 
metallic  state.  However,  this  result  will  enable  us  to  calculate  the 
mobilities  before  and  after  the  transition  from  experimental  measure¬ 
ments  of  the  conductivity. 


B.  Narrow  Band  Limit 

The  transition  metal  oxides  which  exhibit  semiconductor-to- 
metal  transitions  are  characterized  by  extremely  narrow  3d  bands. 
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For  such  materials,  the  effective  masses  of  electrons  and  holes  are 
so  large  that  it  is  not  meaningful  to  use  the  approximations  in  section  A. 
Thu  physical  situation  is  probably  closer  to  the  extreme  limit  of  delta 
function  bands.  We  shall  here  derive  an  expression  for  n(T)  in  this 
limit.  We  assume  one  3d  electron  per  cation,  although  any  number 
can  be  treated  analogously. 

We  return  to  equations  (3.4)  and  (3.5)  for  n  and  p,  but  now  use 
the  narrow  band  densities  of  states: 


p^(E)  ^  N6(E) 

p  (E)  =  N  6(E  +  E  }  . 
v'  g 

Substituting  (3.48)  and  (3.49)  in  (3.4)  and  (3.5),  respectively: 

N 


(3.48) 

(3.49) 

(3.50) 

(3.51) 


The  condition  of  neutrality  leads  to: 


He  nee,  the  concentration  of  free  carriers  is  given  by: 


N 


Putting  (Z.  i)  into  (3,52): 


(3.52) 


(3.53) 


Utlng  the  geneml  expression  for  p,  equation  (3,24),  {3,53)  becomes: 


(3.54) 


First,  we  ensume  »  kT  for  all  ordinary  temperatures.  This 


is  equivalent  to  the  Boltzmann  limit.  Then,  letting: 


X  s  n 


y  =  E  /2kT  , 

?  go' 

equation  (3,  54)  can  be  written: 


x.y  + 


With  the  substitutions: 


Z  s  X  e’ 


W.  exp[4(l  +  t)y  «*>■]  , 


this  equation  becomes: 


z  =  =  w  ^ 


=  ??(W)  , 


(3.55) 


Once  again,  we  find  a  transition  to  the  metallic  state. 

transition  occurs  at  W  =  e^^^,  where: 

o 


4(1  +  6)  e  y  =  e 
^  o 


.  ^  I  1  4(l  +  6>ey 
^o  _  ^Jl+6)e  I  ° 


(3.56) 


Equation  (3.  56)  gives  E^^/kT^  as  a  function  of  6,  The  solution 

is  {dotted  in  Fig.  ni-4.  It  can  be  seen  that  E  /kT  is  not  a  r£  -dlv 

go  o  ' 

varying  function  of  6  over  the  range  of  interest.  It  is  important  to 

note  that  6  is  a  constant  for  each  material.  Thus  E  /kT  will  also 

go'  o 


be  constant.  This  can  be  expressed  as; 


d  In  E 


d  In  T 


(3.57) 


where  X  is  any  external  parameter,  such  as  pressure.  Equation  (3.57) 
can  be  tested  experimentally,  and  is  an  important  prediction  of  this 
theory.  In  the  event  the  energy  gap  is  of  antiferromagnetic  origin, 

6  =  0,  Then  (3.56)  becomes: 


1  j  4ey 


.  i 

4e  y  =  {  e  i 
^  o  " 


(3.58) 


Hence: 


(3.59) 


Equation  (3.  59)  is  true  whatever  the  material,  provided  it  is  anti- 
ferromagnetism  which  was  responsible  for  the  energy  gap.  Just  before 
the  transition,  Z  =  e.  At  this  point,  =  0.067.  Thus  the 

concentration  of  carriers  in  the  conduction  oand  at  the  transition  point 


^  =  0.067  . 

This  small  a  value  of  n  justifies  the  use  of  (2.1)  throughoxit  the 
conducting  region. 


(3.60) 


semi- 


The  approximation  E  »  kT  need  not  have  been  made  in  the 

8 

narrow  band  case.  The  general  equation,  (3.54)  can  be  solved  exactly 
for  the  transition  temperature,  T^.  In  terms  of  x  and  y,  (3.  54)  be- 


c  jmes: 


X  =  [  +  ij  ^3 

where  we  have  set  6  =  0  to  investigate  the  antiferromagnetic  situation. 


(3.61) 


3iv  im  =i,  b  =  i  lor  v: 


1  1  -  X 


(3-  62, 


The  transition  occurs  when  dy/dx  ~  0,  or  when: 

i  -  4x  1  -  X 

o  ,  o 

n  —  • 

o  o  o 

The  solution  to  (3.63)  is  x  =  0,08,  which  corresponds  to  y 

C  '  o 

Hence,  the  complete  Fermi  solution  is; 


|3. 63 j 


3,  56. 


'.  :o  . 


(3,  64) 


As  in  the  eflective  mass  approximation,  Fermi  statistics  tend  to  raise 
the  transition  temperature  somewhat. 

The  solution  to  equation  (3.611  is  shown  in  Fig,  III-5.  The  part 
of  the  curve  for  x  >  l/4  is  unphysicai.  Since  the  energy  gap  becomes 
negative. 

Analysis  shows  that  the  final  result,  (3,64),  is  independent  of 
the  number  of  3d  electrons  per  cation. 

The  case  of  a  Gaussian  broadening  about  delta  function  bands 
can  be  solved  analytically  for  Boltzmann  statistics,  which  is  a  good 
approximation  if  the  spread  is  smai,.  Let  X.  be  the  parameter  which 
measures  the  root  mean-square  deviation  in  E,  Then,  for  one  3d 
electron  per  cation,  we  write  the  densit  =  e3  of  states: 


p  (E)  - - e 

^  Rk 


(:i.  65) 


(E4E  )' 


P  (El 

V 


(3. 66i 


clZ 


equation  (3*  70)  can  be  written! 


—  W 

_ .’7  ^w" 

Z  -  W  =  W  ” 


??{W)  . 


V/e 


The  transition  occurs  at  W  =  e  ■  ,  where! 

2 


4y  e 
'  o 


<2CT-' 

O 


I 

e 


The  solution  ol  equation  (3.72)  is: 


=  2  In  R  q(R) 


where: 


In  R  =  {4  e]  !  -i- 


2k 

a 

Tc  second-order  in  the  small  quantity,  X/2kT  ,  (3,73)  yields; 

C- 

E 


2  [  4  e  (  1  +  -^-2  )  ]  ^  {  exp[  4  e  ( 

o  4k  T 


=1 

+  }. 


4k  T 


3,  7J 


«3. 


(3.73) 


(3.  74) 


It  can  thus  be  seen  that  introducing  a  small  spread  to  the  bands  tends 

to  raise  E  /)f.T^  somewhat.  This  means  that  the  transition  occurs 

at  a  slightly  lower  temperature  than  in  the  zero  bandwidth  case,  as 

one  might  intuitively  expect.  For  X/E  s  0.0'.  E  /kT  =  7.40,  ?he 

go  go'  o 

same  value  as  when  X  -  Oj  for  X/E  *  0.05,  E  /kT  -  7.44. 

go  go'  o 

When  X/E  becomes  much  larger  than  0.  05,  not  only  does 
the  approximation  in  (3.  74)  begin  to  lose  validity,  but  also  we  must 
s'aft  to  take  into  account  the  effective  lowering  of  the  ►eal  band  gap 
due  to  the  finite  bandwidth.  The  parameter  X  is  approximately  1/5 


the  band  width,  so  that  we  may  write  the 
band  gap 


Oi  the  band  Width  to  the 


Band  Width 

~mb<r'c5p'  "  E  ■-  • 

^  go 


When  the  bands  are  about  half  as  wide  as  the  energy  gap,  then  E  /kTQ 

G  ■* 

7.  51  ,  When  the  band  width  and  energy  gap  are  roughly  the  same, 
solution  of  (3.73)  shows  that: 


{3.  7 


We  must  bear  in  mind  that  (3.  73)  is  valid  only  for  extremely  narrow 
bands,  so  that  we  cannot  allow  \  to  become  very  large.  Nevertheless, 
when  the  band  width  and  the  band  gap  are  both  small  in  absolute  mag¬ 
nitude  (e.g.  of  the  order  of  0.05  eVi  equa+^ion  (3,75)  could  have  meaning 
In  the  case  of  antiferromagnet. sm,  the  theory  presented  here 
gives  a  behavior  for  the  variation  oi  sublattice  magnetization  with 
temperature  which  differs  considerably  from  that  calculated  using  the 
ordinary  molecular  field  model  of  ant  1*  r  rmagnerism.  We  can  eva  uat^ 
the  subiattice  magnetization  in  the  narrow  band  limit,  using  equations 
(E.20),  (2,225,  and  (3.  55).  Equation  ’2.20),  which  gives  the  average 
value  of  sp.n  .n  the  valence  band,  provides  the  magnetization  at  T  "  0: 

M^(0)  -  gPg  N  <S^> 


r  "/a 


XT  ^  1 

I 


IT 


^0 


go 


(1  +  cos  ka)] 


2  ‘ 


I  M 

^  •-  g  kgN 

go 
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where  MpCO)  is  the  magnetization  at  T  “  0  on  the  eublatiice  where  the 
§pin  IS  primarily  down,  g  is  the  spectroscopic  splitting  factor,  and 
-  efi/2mc  is  the  Bohr  magneton.  In  the  limit  ==  0,  M^iO)  Is 
just  the  value  obtained  by  the  molecular  field  treatment  of  antiferro- 
magnetism.  For  all  finite  band  widths,  however,  MpfO)  is  reduced 
in  magnitude. 

Equations  (2.22)  and  (3.76)  give  the  sublattice  magnetization 
as  a  function  of  the  number  of  free  carriers: 


E.  j  _  E.  sin  ^ 

(  1  -  (-^  +  (3^77 

go  go 

Finally,  equations  (3.55)  and  (3.77)  can  be  used  to  calculate  M  as  a 
function  of  temperature: 


M 


r  ^  ?  E^  sinrTrx(y)l 

pCT)  ^  g  Jig  n{  (  1  ^  +  x(y)]  +  -g - - -  }  (3.  78 

go  go 

where  y  s  E^^/ZkT  as  before,  and  x(y:i  is  given  by  (3.55),  Equation 

(3.  78)  is  valid  only  when  E.  «  E 

O  go 

We  shall  investigate  the  subiattice  magnetization  for  two  cases, 
E^  0,  th  -icrcw  band  limit,  and  E^  E  /4,  just  abou^  the  upper 
limit  for  the  use  of  (2,20)  and  (2.22),  When  E^  —  0,  the  magnetization 
iS  given  by; 


M 


(T)  =  g  Pg  N  [  -^  +  x(y)]  , 


(3.  79 


When  =  E  /4,  the  magnetization  is: 

o  ' 


M  (T)  s  g  |q[  +  x(y)]  , 


(3,80 


From  (3.  79)  and  (S.SO)^  the  ratios  of  to  the  magnetization  at 

P 

T  -  0  are: 


M  (T) 

T:™  =1-2  xiy) 

P 

Mp(T) 


(3.  8i) 


-  1 


8 


x(y)  • 


Evaluating  (3.81),  we  find  that  the  ratios  remain  very  close  to  unity 
from  T  -  0  until  T  ~  0.  7  Tj^.  As  T  is  increased  above  this  point,  the 
ratio  begins  to  fall  off,  until  just  below  Tj^,  where: 


Mp(01 


Mp(T-4 


=  0.87 


0.82 


=:  0) 


(3.82) 


At  the  subiattice  magnetizations  drop  sharply  from  the  values 
given  by  (3.82)  to  zero.  This  behavior  for  M^T)/Mp(0)  contrasts 
with  the  Brillouin  curve  obtained  using  the  molecular  field  theory  of 
antiferromagnetism  in  that  it  remains  higher  at  low  temperatures  and 
then  drops  to  zero  at  much  more  drastically.  As  the  bands  get 
Wider,  the  ratio  Mp(T)/Mp{0)  approaches  the  molecular  field  behavior 
more  closely;  however,  as  we  have  seen,  the  magnitude  of  Mp(0)  is 
greatly  reduced  from  the  zero  band  width,  and  thus  the  molecular 


field  case. 


Chapter  IV 


COMPARISON  OF  THEORY  WITH  EXPERIMENTAL  RESULTS 


In  this  chapter,  we  shall  present  the  results  of  experiments 
performed  on  the  four  materials  known  to  exhibit  semiconductor-to- 
metal  transitions^  and  compare  the  results  to  the  predictions  of  the 
theory  given  in  Chapters  n  and  IH. 


A.  V,Oj 


V^O^  will  be  of  primary  interest  to  us  s^nce  it  is  the  material 
which  has  been  studied  experimentally  in  the  greatest  detail.  At  room 
temperature*  V^O^  has  corundum  structure,  with  rhombohedral 
symmetry.  The  lattice  parameters  are  most  easily  expressed  in  terms 
of  a  hexagonal  unit  cell,  consisting  of  6  molecules.  As  measured  by 
Warekols  [36]  at  300®K,  the  hexagonal  c-axis,  14.  00  A,  whereas 

-  a  OR  1 

H 

with  increasing  temperature  as; 
d  In  c, 

“dT 


the  basal  plane  lattice  parameter,  a„  "  4,95  A.  The  c-axis  contracts 


^  =  -5.8  X  10“^ 


while  the  a-axis  expands  as; 
d  In  a 


dT 


^  ^  22.9  X  10"^  (°K)"^ 


Thus  the  thermal  expansion  coefficient  is; 

^  40,0  X  10"^  (°K)‘^  . 

in  the  vicinity  of  ISO^K,  a  phase  transformation  occurs.  At 
lower  temperatures,  has  monoclinic  symmetry.  As  also 


measured  by  Warekois  [36]  at  7?^K,  the  monocUnic  unit  cell*  con¬ 
sisting  of  12  rrsolecules,  can  be  expressed  as: 
a  =  8. 57  A 
b  =  4.98  A 
c  ’  13.88  A 
p  =  91.  6°. 

A  simplified  diagram  of  both  structures  of  'V’20j  is  shown  in  Fig.  IV- 1 . 

The  monoclinic  distortion  can  be  thought  of  as  a  shifting  of 
pairs  of  cations  in  the  hexagonal  basal  plane  towards  one  another*  re¬ 
sulting  in  an  effective  tilting  of  the  c-axis.  From  the  above  data  at 
77’^K,  the  parameter  €,  defined  in  Chapter  II,  can  be  evaluated  as: 

c  sin  (P  -  :r/2) 

IZ  0.0323  A 

^  a  '  2.86  A 

3 


-  0.011  . 

Dilatometric  data  by  Foex  [  16]  showed  that  there  is  a  sharp 
contraction  of  volume  at  the  transition  temperature.  Recent  experi¬ 
ments  by  Minomoru  and  Nagasaki  [37]  have  shown  that  this  contraction 
is  about  3.  5%. 

The  presence  of  antiferromagnetism  has  never  been  established 
in  Carr  and  Foner  [38]  measured  Xg  and  Xj  with  respect  to 

the  c-axis,  and  found  both  components  the  same  order  of  magnitude 
and  essentially  constant  down  to  4*^K.  Paoletti  and  Pickart  [39]  per¬ 
formed  neutron  diffraction  experiments  which  placed  an  upper  limit 
on  the  antiferromagnetic  moment  of  0.5  Bohr  magneton  per  cation. 
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Abrahams  [2iJ  was  able  to  detect  antifer romagnetism  in  Ti^O  by 
means  of  neutron  diffraction,  and  found  a  moment  of  0,  2  Bohr 
magneton  per  cation,  so  antiferromagnetism  is  certainly  not  as 
yet  ruled  out.  Wucher  [40]  and  later  Teramshi  and  Terama  [4  ] 
have  interpreted  their  magnetic  susceptibility  measurements  to  infer 
that  anliferromagnetic  clear  through  the  transition,  up 

until  the  vicinity  of  525°K;  we  shall  refer  to  this  point  as  the  high- 
temperature  transition. 

Transport  measurements  on  have  been  done  by  manv 

worker^.  Foex  [16]  found  that  at  the  same  temperature  as  the  dilate- 
metric  anomoiy,  there  was  a  semiconductor-to-metal  transaction.  Be¬ 
low  the  transition  temperature,  the  resistivity  increased  exponent  ally 
With  j/T,  the  activation  energy  being  0.20  eV.  Above  the  transit  on 
and  up  to  the  high-temperature  anom.aly  at  525°K,  the  resistivity 
increased  linearly  with  T,  and  was  metal -like.  There  was  a  jump 
in  conductivity  of  a  factor  of  10^  at  the  transition.  Foex  found  that 
above  525  K,  the  resistivity  decreased  exponentiailv  once  agam  wi»h 
an  activation  energy  of  0.  04  eV. 

Monn  s  experiments  [20]  (see  Fig,  J-1)  confirmed  most  of 
Foex’s  results.  In  Morin’s  samples,  the  jump  in  conductivity  was 
alt...  a  factor  of  10  ,  and  his  measured  activation  energy  in  t  he  semi¬ 
conducting  region  was  0.  17  eV.  MacMillan  [42]  prepared  a  number 

of  samples  of  stoichiometry  varying  ^rom  V  0_  to  V.,0-  and 

2  2.  yU  2  3. 05 

found  virtually  identical  conductivity  behavior  for  all  samples.  The 
activation  energy  was  0.  13  eV,  independent  of  the  composition 
Recently,  Feinleib  [21]  performed  measurements  on  single  crystals 


of  found  to  be  within  0.3%  of  stoic hiOjUetry.  He  found  a  jump 

£r  J 

8 

of  10  in  conductivity  at  the  transition  temperaturej  measured  to 
be: 

T  ^  152"‘K  .  (4.  1 

o 

The  magnitude  of  the  activation  energy  varied  from  0*  12  to  0,  18  eV. 

Other  electrical  conductivity  work  has  been  done  on  by 

Austin  [43]  ,  and  by  Goodman  [44]  .  All  experiments  confirm  that  a 

semiconductor-to-metal  transition  occurs  with  near  or  at  152°K, 

6  8 

that  it  is  a  sharp  transition  with  jumps  in  s’  of  10  -10  ,  that  above 
the  behavior  is  metallic,  and  that  below  the  conductivity  is 
thermally  activated,  with  activfvtion  energy  in  the  range  0.  ll  to 
0.  20  eV.  These  results  show  that  ^2^3  *  material  to  which  the 

theory  developed  in  Chapters  II  and  III  may  be  applicable,  and  they 
provide  the  value  of  T^,  equation  (4.  1).  Where  there  is  a  choice,  we 
shall  use  the  results  of  Feinleib  [21],  since  his  data  were  taken  on 
single  crystals,  and  his  material  was  closest  to  stoichiometric.  The 
values  obtained  for  the  activation  energy  conceivably  could  give  us 
some  information  about  the  size  of  the  energy  gap.  In  particular,  if 
the  mobility  were  independent  of  temperature,  our  previous  results 
would  give: 


r  ”  n  e  p 


(4.2) 


Since  we  ’  ve  shown  that  n/N  does  not  become  very  large 


before  the  vransition,  (4,2)  would  imply: 


However,  it  is  unrealistic  to  assume  a  temperature -‘independent 

mobility.  For  example,  if  polar  scattering  is  important,  then  the 

mobility  would  be  thermally  activated,  giwing  a  negai,^ve  contribution 

to  the  measured  activation  energy  of  approximately  k  where  © 

D  D 

is  the  Debye  temperature.  Since  the  Debye  temperature  is  in  the 

vicinity  of  500°K  [l5] ,  this  would  modify  (4.  3)  to  giv^e? 

E.  -  E  /2  -  0.  04  eV. 

A  go' 

Other  types  of  scattering  may  be  present,  such  as  spin-disorder 
scattering  or  impurity  effects.  Also,  if  polaron  formation  is  im¬ 
portant  as  we  expect  in  the  case  of  narrow  bands,  there  may  be  a 
polaron  contribution  to  the  apparent  act^vatiOn  energy.  In  genera, 
the  most  we  can  say  is: 


E  /2  4  E 
go 


X 


(4. 4 


where  E^  contains  all  contributions  besides  that  of  the  energy  gap, 
and  may  depend  on  the  particular  sample  being  investigated.  Only  ^f 


lattice  scattering  dominates  can  be  negative,  all  other  con^r  bu  = 
Uons  are  positive  or  zero.  li  we  as.i.ume  E  pcsiiive.  the  meas- 

X 

ured  range  of  activation  energies  of  0.  j  ^  eV  to  0,20  eV  gives  us  an 


upoer  Kmit  for  E  : 

go' 

E  <0,22  eV  . 
go 


*  4 .  5 ' 


In  order  to  find  the  exact  value  of  E  .  Feinieib  [24]  per- 
formed  optical  experiments  on  V^O^  at  temperatures  above  and  b* 
the  transition.  No  transmission  was  observed  m  the  room  temperature 


ICiOW 


5U 


metallic  region  betweeri  u,  05  eV  and  6  eV.  However,  m  the  semi¬ 
conducting  region,  at  77“K,  transmission  was  obtained  between 
0.  1  eV  and  0,  4  eV.  The  transmission  changed  by  several  orders  of 
magnitude  in  this  range,  and  clearly  indicated  the  presence  of  an 
absorption  edge  at  approximately  0,  I  eV.  Since  at  77®K,  T  -  l/E  T 
and  we  have  shown  that  n/N  is  negligible  at  this  point,  this  value 


can  be  taken  as  E  for  V.O*: 

go  2  3 

E  ^0,  lOeV. 

go 

From  {4.  1)  and  (4.  6),  we  find: 
E 


^2^3’  expect  the  narrow  band  analysis  of  section  ITfB 
to  be  appropriate.  If  the  energy  gap  m  were  of  antiferromag¬ 

netic  origin,  equation  (3.59)  should  apply.  The  result,  (4.7),  agrees 
with  this  well  within  experimental  error.  If  the  gap  were  due  to  a 
crystalline  distortion,  then  as  we  have  seen,  (3.59)  still  remains 
approximately  true.  We  can  then  use  the  experimental  result  (4.  7) 
to  evaluate  p.  This  can  be  done  directly  from  (3,56),  or  from 
Fig.  ni-4.  The  result  isr 


P  =  4.4Eg^/N. 

There  is  no  evidence  that  is  antiferromagnetic.  The  presence 

of  a  phase  transformation  of  the  correct  type  at  T^  is  strong  evidence 
that  the  gap  arises  from  a  crystalline  distortion.  If  so  the  theory 
of  Chapters  II  and  IH  would  account  for  the  phase  transformation  as 
well  as  for  the  semiconductor-to-metal  transition.  We  can  gain 


further  insight  into  the  nature  of  transition  from  the  resuk  of 
the  thermodynamtc  argument  in  section  II  B: 


2  as 

Feinleib  [2ij  has  measured  the  total  pressure  coefficient-  ■'j 


dE 


F  d 


3  E  3  E 

/  8=  /  g  =  <hi 

T  "  H  n  -P,  T  dP  *  ^ 

Assuming  Boltzmann  statistics  apply,  the  concentration  of 
carriers  is; 


n  =  n  e 
o 


E 

_ ^ 


(4.  "0) 


Recall; 


-  V  P  • 

g  go  ^ 

Substituting  (2.9)  into  (4.  10): 


(2.9) 


n  =  n  e 
o 


go  pn  V  P 
2kT  2Fr  2kT 


H.  M 


Thus; 


\  ,  P  dn 

'SP  ”  2kT  dP  ’ 


H.  2/ 


dn 


Solving  (4,  12)  for  ^  . 


dn 


I3L 


dP  Zk'f  ~  np 


(4.  13) 


But  from  (2.  9l‘ 
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‘  “  V  -  P  C4.  14) 

Substituting  (4.  13)  and  (4,  14)  into  (4.  9): 

‘  (4.13) 

At  T  =  132°K,  Feinleib  rfieasurea  ^  =  2.2  X  10“^  eV  bar"^  Hence, 
from  (4.  15): 

y  *  1.  8  X  10“^  eV  bar"*  .  (4.  I5) 

Unfortunately,  the  isothermal  compressibility  has  not  been 
measured  for  V^O^.  But  the  value  X  =  0.  58  X  10"^  bar"*  has  been 
found  for  the  very  similar  materials  Fe^O^  and  FeTiO^  [45] ,  and 
the  compressibility  of  will  undoubtedly  not  be  much  different. 

Using  this  value  for  "X  in  equation  (2.  1')  yields: 

3 .  6  E  a  jc.* 

N  "  *  B  n  *V,  T  {4.17) 

Thus  the  measured  value  for  S,  equation  (4,8),  is  almost  entirely 
accounted  for  by  the  variation  of  E  due  to  explicit  volume  changes. 

In  fact,  with  just  a  very  slight  band  spread,  such  as  that  due  to  spin 
disorder  broadening  (see  Chapter  Vi),  this  variation  alone  is  suf¬ 
ficient  to  obtain  the  experimental  value  (4.7).  This  result,  taken 
together  with  the  measurement  of  a  finite  volume  change  at  T  by 
Minomura  and  Nagasaki  [37],  is  another  indication  that  the  transition 
is  due  to  crystailme  structure  effects  rather  than  antlferromagnetism. 

No  matter  what  the  cause  of  the  gap,  the  theory  predicts 
that  the  ratio  E^^/kT^  will  remain  constant.  We  have  expressed 
this  as: 


d  In  E 


d  In  T 


In  order  to  test  equation  {3*  57),  Fei.nleib  [21]  varied  both 
^go  ^PP^yi^g  hydrostatic  pressure  and  uniaxial  stress. 

His  results  were  as  follows.  The  transition  temperature  varied  w^th 
pressure  as: 


d  In  T 


2.  6  X  10“^  bar“^ 


(4.  18} 


where  P  is  the  hydrostatic  pressure.  With  uniaxial  stress  applied 
parallel  to  the  corundum  structure’s  b-axis,  the  relation  between 
transition  temperature  and  stress  was  determined  to  be: 


d  in  T 


=  -2.8  X  10"^  bar"^ 


(4.  19 


where  S  is  the  stress.  With  stress  applied  along  the  c-axis,  the 
uniaxial  stress  coefficient  was  at  least  an  ord^r  of  magnitude  smaller: 


d  In  T 


“•  <  0.  3  X’  10  ^  bar 


H.20) 


The  pressure  and  uniaxial  stress  coeffic.ents  for  the  energy  gap  were 
determined  from  the  stress  dependence  of  the  activation  energy  for 
electrical  conductivity.  These  results  can  be  expressed,  in  the 
case  of  hydrostatic  pressure,  as: 


d  In  E 


^  -  2.  2  X  10“^  bar' 


( 4 .  2  i  } 


whereas  with  uniaxial  stress  along  the  h-axis* 


d  In  £ 


=  -  3.  0  X  10"^  bar""  . 


(4.22i 


In  Chapter  VI.  we  show  that  the  pressure  coefficient  of  conir .butions 


to  the  activation  energy  from  polaron  effects  and  non* stoichiometry 
are  of  tne  order  of  10  ^  bar  and  thus  these  contributions  do  not 
affect  the  results  (4.21)  and  (4.22), 

A  com|>arison  of  (4,21)  with  (4.  18)  and  (4.22)  with  (4,19) 
shows  that  equation  (3.  5")  is  indee  satisfied  within  experimental 
error.  Thus  the  two  major  predictions  of  the  theory,  equations  (3,  S’?) 
and  (3.59)  are  verified  by  Feinlcib's  results.  There  is  no  a  prion 
reason  for  the  validity  of  (3.  57)  and  (3.  59),  and  no  other  model  here¬ 
tofore  suggested  predicts  these  relations.  Therefore  the  agreement 
with  experiment  must  be  considered  as  good  evidence  for  the  applica¬ 
bility  of  the  model  of  Chapters  II  and  in. 

The  anisotropy  in  variation  of  T^  with  uniaxial  stress  found 
in  (4.  19)  and  (4,  20)  provides  still  another  indication  that  it  is  crystal 

structure  changes  which  bring  about  the  energy  gap.  The  distortion 
which  doubles  the  number  of  cations  in  each  unit  cell  is  entirely 
in  the  basal  plane  of  the  corundum  structure.  Therefore  the  varia¬ 
tions  of  and  with  stress  applied  along  the  c-axis  should  not 
be  very  great,  as  is  borne  out  by  (4.  20).  However,  the  changes  in 
and  T^  with  stress  applied  along  the  a -axis  or  b-axis  are 
strikingly  large,  as  we  would  expect. 

Feinleib  [21]  also  investigated  the  high-temperature  anomaly, 
measuring  resistivity  vs.  temperature  from  300°K  to  800°K.  He 
found  that  resistivity  increases  V  nearly  with  temperature  with  the 
same  slope  both  below  and  abo  e  the  high-temperature  transition. 
However,  in  the  vicinity  of  550*"  K,  the  resistivity  undergoes  a  rather 
sharp  anomalous  increase.  In  the  region  T^  <  T  <  500°K,  the 


measured  resistivity  can  be  expressed: 

Pl  (T)  =  4,  3  X  10“^  ohm-cm  {0.  51  +  0, 49  T/T  ).  {4.  23 

If  we  write  for  the  total  resistivity  up  to  800°K: 

p{T)  =  pj^(T)  4  p^(T)  (4.24 

where  P^CT)  is  an  extrapolation  of  (4.23),  and  p^(T)  is  the  anomalous 
resistivity,  then  P^{T)  is  a  function  which  is  zero  up  to  500®K,  then 
sharply  rises  to  a  value  of  12  X  lO***  ohm-cm  by  600°K,  above  which 
temperature  it  remains  constant.  Such  behavior  for  p^(T)  bears  a 
striking  resemblance  to  the  spin-d’sorder  resistivity  calculated  by 
De  Gennes  and  Friedel  [46] ,  and  therefore  suggests  a  magnetic  order¬ 
ing  temperature  ot  600^K.  As  we  have  mentioned,  Wucher  [40]  and 
Teranishi  and  Terama  [41]  have  concluded  that  this  is  the  Neel  tem¬ 
perature  of  V20^.  We  shall  return  to  this  pmnt  later. 

B.  VO 

The  experimental  results  for  VO  are  quite  similar  to  those  for 
^2^3’  salt  structure  above  T  =  126®K.  At  thus 

temperature,  there  is  a  crystalline  structure  distortion  to  orthorhombic 
sym.metry.  The  exact  low  temperature  crystal  structure  has  not  as 
yet  been  reported.  Little  is  known  about  the  magnetic  properties  of 
VO  as  well.  Neutron  diffraction  and  magnetic  susceptibility  measure¬ 
ments  have  not  been  performed. 

What  we  do  Hhow  about  VO  is  its  electrical  properties.  Morm  [20] 
isee  Fig.  I-i)  found  that  there  is  a  sharp  semiconductor-to-metai 


transition  at: 


86 


=  126'’k  .  (4.  25) 

The  Jump  in  conductivity  *t  T  was  meaeured  to  be  a  factor  of  10°. 

Below  the  activation  energy  was  “  0. 14  eV.  Above  T^,  the 
resistivity  increases  linearly  with  T  and  is  metal-like.  It  can  be  seen 
from  Fig.  I-l  that  the  electrical  properties  of  VO  greatly  resemble 
those  of 

Austin  [43]  has  performed  pressure  experiments  on  VO,  and 
has  observed  the  effects  of  quasi-hydrostatic  pressure  on  the  electrical 
properties  of  single  crystals  of  variable  composition.  From  meas¬ 
urements  of  resistivity  as  a  function  of  pressure  at  94®K,  Austin  ob¬ 
tained; 

dE  . 

-gpS  s  .2.9  X  10“  eV  bar"  .  (4.26) 

Since  no  direct  measurement  of  the  energy  gap  hi»  s  been  made  in  the 
case  of  VO,  we  cannot  do  more  than  assume  the  gap  is  the  same  as 
that  of  V^O^  as  a  first  approximation,  and  perform  a  self-consistent 
calculation  using  (2. 13)  and  (3.  56)  to  obtain  a  better  value.  Thus,  we 
begin  with: 

(4.27^ 

We  use  this  value  to  calculate  the  contribution  to  p  of  the  volume  term 
in  {2.  13).  In  the  case  of  VgO^,  this  term  was  sufficient  to  give  the 
correct  ratio  01  ^gQ  ^  VO  is  completely  analogous  to  ^20^, 

as  Fig,  I-l  suggests,  we  should  obtain  good  results  in  this  manner. 

From  (4.26): 

2.9  X  10“^  eV  bar’^  . 
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Using  (4,  15); 

7  =2.7xlO-‘eVbar-^  ,4.  ,g, 

Fo  -  the  isothermal  compressibiliiY*  the  best  value  available 
is  that  for  ZnO  [45] : 

"K  =  0.  78  X  10'^  bar*^  .  29': 

Substituting  (4,28)  and  (4.29)  into  (2*  13),  we  can  evaluate  the 
contribution  due  to  explicit  volume  changes  as: 

6.  9  E 


P  - 


N 


(4.  30j 


Using  this  value  for  p  in  the  general  expression  for  narrow  bands, 
equation  (3.  56),  we  find: 


o 


,’4.  3 


But  (4.341  implies  that  E^^  0,  10  eV,  as  was  assumed  in  (4.27#.  Thus 

we  see  that  this  assumption  is  se.f-consistent,  and  that,  just  as  for 

^2^3'  volume  term  in  (2.  13/  and  narrow  band  theory  give  the 

correct  ratio  of  E  to  T 

go  o 

However,  now  we  can  use  Austin's  results  for  dT^/dP  to  cherk 
the  vaxidity  of  (3.57).  This  relationship  provides  an  important  test 
of  the  .hecry.  From  (4.  26)  and  (4.  27#^  we  can  express  Austin's 
measurements  on  the  pressure  variation  of  the  energy  gap  as: 
d  In  E 


“dF 


=  -29  X  10'^  ba 


6  ,  -1 


(4.32) 


As  we  mentioned  before,  we  shall  show  m  Chapter  VI  that  the  pressure 
coefficients  of  contributions  tc  E^  from  pciaron  effects  and  non- 
stoichiometry  are  of  the  order  of  10”^  bar  '. 


so  that  (4.  32)  probahiV 


does  represent  only  the  variation  of  the  energy  gap. 

At  94®K,  Austin  found: 
dT 

-  -  3  X  10"^  bar“*. 

Hence: 

d  in  T 

— gp —  =  -  32  X  lU  ^  bar**  . 

Compartaon  of  {4.  32)  with  (4. 33)  demonstrates  thai  (3.5?)  is  satisifec 
Within  experimental  error. 

3  j  conclude  that  it  appears  likely  that  the  energy  gap  in  VO 
arises  from  a  crystalline  structure  distortion  and  has  a  value  of  about 
0.  lO  eV,  approximately  th  same  as  tha^  of  V^O^.  The  lower  transi¬ 
tion  temperature  of  VO  seems  to  be  due  to  a  greater  change  in  volume 
at  T^.  This  is  reflected  in  a  larger  va^ut-  for  p  in  equation  (2.  I). 
Outside  of  the  slightly  lower  transition  temperature,  the  behavior  of 
VO  and  seems  identical. 

c.  VO^ 

The  crystalline  structure  of  VO^  is  relatively  simple  and  well- 
known.  Above  T  =  340°K,  the  structure  ^s  tha^  o*  rwtiJe,  the  cations 
occupying  the  positions  of  a  body-centered  tetragonal  lattice.  Below 
340  K,  the  symmetry  is  monociimc,  the  structure  being  that  of  MoO 
I  he  low  temperature  phase  is  just  a  distorted  rutile  structure  — 
the  cations  which  in  the  rutile  phase  were  collinear  and  separated  by 
2.8?  A,  are  slightly  non-collinear  and  spaced  alternately  2.65  A  and 

ft 

3.  12  A  apart  [47l .  This  is  an  almost  classic  example  of  the  model 
of  section  HD  —  the  umt  cell  is  doubled  by  a  distortion  of  the  lattice 


33) 


in  one  dimension,  alternate  cations  pairing.  The  parameter  e  ,  defined 
in  II D,  is: 

1  2. 65  A 

€  -j - T 

5.  77  A 

=  0.04  . 


The  details  of  the  transformation  have  been  s*-ud<ed  by  Mmomura  and 
Nagasaki  [37]  .  They  find  that  at  T^  the  n  ^noe,  me  a=ax  .8  ccn'ra.ts 
from  5,77  A  to  5.70  A,  twice  the  rutOe  sti  cture's  c-axiS.  simn  - 

■3  9 

taneously,  the  monoclimc  b-axis  expands  from  4.50  A,  to  4,  54  A,  the 
rutile  a-axis,  while  the  monocllnic  c-axis  contracts  from  5.39  A  to 
5.24  A  Z/ fS  of  the  rutile  a-axis.  They  found  no  measurable  total  vcium’? 
change  at  the  transition  pc^nt,  just  a  change  ‘n  ’he  slope  thermal  expan¬ 
sion  coefficient).  They  concluded  that  the  t  ransformaticn  was  of  sec  ond 
order,  Jr/fray  and  Dumas  [  18]  measured  a  vo  ume  ce.ntiaction  of 
0.08%  on  powdered  VO2.  which  is  consisten’  with  the  la^’  re  constant 
measurements  of  Minomura  and  Nagasaki.  In  any  event,  the  volume 
change  .s  more  ♦han  an  order  of  magn  tude  smaller  than  that  whlc  h  occur 
in 

Magnet’c  susceptibility  data  have  been  taken  on  by  Rudorff 

et  al  [48]  and  by  Kawakubo  and  Nakags  [49]  .  Both  sets  of  experiments 
found  susceptibility  independent  of  temperature  both  above  and  below 
340®K,  with  a  jump  in  at  340°K  of  a  factor  of  rougluy  7.  They  con¬ 
clude  that  the  low  temperature  susceptibility  can  be  accounted  for  by 
temperature  independent  paramagnetism,  find  that  there  is  no  evidence 
for  antiferromagnetism.  Similarly,  Kasper  [50]  found  nc  magne’iC 
scattering  in  neutron  diffraction  measurements  below  340°K. 
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The  electrical  conductivity  behavior  of  VO,  was  also  studied 
by  Morin  [20]  (see  Fig.  I-l),  who  found  a  sharp  semiconductor -to - 
metal  transition  at: 

Tq  =  340®K  (4,34) 


the  temperature  at  which  the  phase  transformation  takes  place.  Below 

T^,  the  activation  energy  was  0. 13  eV,  The  jump  in  conductivity 

2 

was  a  factor  of  10  ,  Neuman,  ^awson,  and  Brown  [51]  found 
tbr  same  value  for  T^,  but  obtained  a  jump  in  s  at  T^  of  10^  and  an 
activation  energy  of  0.44  eV. 

The  energy  gap  has  not  as  yet  been  measured.  The  narrow 

band  limit  appears  to  have  given  good  results  in  the  cases  of  V>0,  and 

VOj  we  expect  it  shovdd  be  at  least  as  appropriate  for  VO^,  since  the 

cations  are  V  ,  which  are  smaller  in  spatial  extent  than  are  V 

and  It  can  be  concluded  that  the  cause  for  the  higher  transition 

temperature  in  VO,  is  most  likely  a  larger  value  for  E  .  This  can 

i  go 

be  easily  explained  by  the  four  times  larger  alue  obtained  for  the 
distortion  parameter  €  In  the  case  of  VO^,  which  should  certainly 
produce  a  much  larger  gap.  In  our  analysis  of  section  II  D,  we  found 


that  using  a  delta  function  potential,  the  larger  the  initial  gap,  the 
larger  the  value  obtained  for  p.  The  most  reasonable  conclusion  Is 


that  the  final  E^^/kT^  ratio  for  VO2  is  probably  somewhat  larger  than 
it  is  for  ^20^.  As  a  rough  estimate,  we  shall  take: 


10  kT 

o 


0.3  eV  . 


Pressure  measurements  on  VO^  have  been  carried  out  by 


Minomura  and  Nagasaki  [37]  and  by  Nenman  et  [51]  -  Minomura  and 
Nagasaki  measured  the  variation  of  T  with  pressure  and  found: 


d  In  T 

o 


1.4  X  10 


.6  .  -1 
bar 


14.35) 


This  is  a  very  small  value  and  contrasts  strikingly  with  the  pressure 
coefficients  of  ^^2^3  VO  [see  (4.  18),  (4.19),  and  (4.33)],  which 
are  an  order  of  magnitude  larger.  Neuman  et  al  could  not  find  a  shift 
of  with  pressures  up  to  6  kilobars  within  the  0.  scatter  in 
itself.  Thus  they  find  a  still  smaller  value  for  3  In  T^/dPthan  iS  given 
in  (4.  35). 

Neuman  et  al  also  measured  the  change  in  activation  energy 
with  pressure.  Their  experiments  show; 


dE. 

A 

“dp- 


-  5.  0  X  i0~^  eV  bar 


f4.  36i 


From  their  measured  value  of  =  0.44  eV,  (4,  36)  yields: 

d  In  E  .  , 

^  r  -  1.  1  X  10“*'  var”"  . 


"dP 


This  is  also  an  order  of  magnitude  smaiier  than  the  values  for  V^O^ 

and  VO,  given  in  equations  (4,  21),  (4.  22),  and  (4,  33).  However,  in 

this  case  we  cannot  use  (4.36)  to  evaluate  d  In  E  /dP  in  order  to 

go 

compare  with  (4.  35)  and  check  the  valid  ty  of  (3.57).  For  equation 


(4.4)  gives: 

■JP" 


(4.  37) 


In  the  cases  of  V^O^  and  VO.  the  measured  values  of  dE^/dPwerc 
so  large  that  we  could  neglect  dE^/dP  in  (4.  36).  As  we  ment-iOned 
previously,  we  show  in  Chapter  VI  that  dE^/dP  is  normally  of  tht 


Order  of  10  eV/bar.  But,  in  the  case  of  this  is  just  the 

order  of  magnitude  measured  for  dfl^/dP  in  (4.  3tt.  A  further  com¬ 
plication  in  the  results  of  Neuman  et  al  is  the  extremely  high  value 

found  for  E  .  ^nce  we  expect  E  -  0.3  eV,  (4.4)  shows  that 
*  go 

\  much  larger  than  in  V^O.  or  VO.  Therefore,  for  VO  , 

we  can  only  use  (3.57)  and  (4.35)  to  conclude: 
d  In  E 

dP*®  =  -  1.4  X  10“  bar"^  .  (4,38 

Using  the  approximate  value,  E  ~  0.3  eV,  (4.38)  shows: 

& 

dE 

-  .  4  y  0“^  eV  bar"*  . 

We  can  use  this  approximation  and  equation  (2. 13)  to  calculate  the 
contribution  to  ^  from  volume  changes.  Since  the  volume  change  at 
the  transition  is  so  small  in  VO^,  we  expect  this  term  also  to  be  very 
Small.  Taking  for  the  compressibility  the  value  measured  [45]  for 
the  rather  similar  material.  TiO,,  X  -  0.  59  X  lO’^  bar"\  equation 
(2,  13)  gives: 

P  =  0.  05  E  /N  -  O  E  /Bn\^  ^  . 

go'  go'  V,  T 

the  contribution  due  to  explicit  changes  in  volume  is 
negligible,  as  we  expected.  The  transition  in  VO^  seems  to  be  due  to 
a  crystalline  distortion  with  little  or  no  accompanying  volume  change. 
Thus  hydrostatic  pressure  is  actually  a  poor  variable  to  use  in  the 
study  of  this  material.  B;  should  be  much  more  advisable  to  study  the 
changes  in  electrical  conductivity  with  uniaxial  stress. 


In  crystal  structure,  Ti^O^  appears  to  resemble  con¬ 
siderably,  Pearson  [19]  studied  from  300^K  to  650^K  and  found 

that  the  structure  was  that  of  corundum,  as  in  V_0,.  At  300°K,  the 


lattice  parameters  are  c, 


13.64  A.  a„  =  5.  15  A,  The  c-axiS  ex- 
M. 


pands  and  the  a-axis  contracts  with  increasing  temperature,  slowjy. 
except  in  the  neighborhood  of  450-550*^K,  At  b50°K,  the  lattice 
parameters  of  Ti^O^  and  are  virtually  the  same.  Abrahams  [23] 

examined  the  structure  of  Ti^O^  down  to  4.2°K  and  found  the  lattice 
constants  to  be  within  0.  5%  of  those  at  300°K.  Thus  there  is  no  phase 
transformation  in  Ti^O^,  and  the  symmetry  is  rhombohedral  at  all 
temperatures.  The  only  region  resembling  an  anomalous  one  ts  the 
range  450-550*^K,  where  the  thermal  expansion  parameters  sharp  \ 


increase. 


The  magnetic  structure  of  Ti20j  is  known  unambiguously. 
Abrahams  [23]  performed  neutron  diffraction  experiments  over  the 
♦emperature  range  * .  4'^K.  to  lll^K  and  found  that  antiferromagne^ism 
was  present  until  a  Neel  temperature  in  the  vicinity  of  600° K.  The 
magnetic  structute  was  monoclinic,  strikingly  resembi'ng  the  lew 
temperature  crystal  structure  of  V^O^,  Abrahams  determined  that 
the  c-axiS  pairs  were  antiferromagnetlcally  aligned,  whereas  the  basal 
plane  spins  were  all  parallel.  The  spins  were  perpendic  ilar  to  the 
c-axis,  thus  reducing  the  symmetry  to  monoclinic.  The  magnitude 
of  the  antiferromagnetic  moment  was  about  0.2  Bohr  magneton  per 


c  aiion. 


Pearson  [19]  and  Foex  and  Wut  her  [52]  have  measured  magnetic 
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found  a  jump  in  conductivity  of  a  factor  cf  40  at  about  550  K,  approxi- 

mately  the  experimental  Neel  temperature.  The  variation  in 

is  unique  to  Ti^O^,  and  prevents  us  from  deciding  on  an  exact  value 

for  T  .  The  most  we  can  say  is: 
o 

T  =  450-550°K.  (4,41) 

o 

Since  there  is  no  crystalline  distortion  in  the  case  of 
it  must  be  the  antiferromagnetism  which  is  responsible  for  the  energy 
gap.  Thus  the  theory  of  section  n  C  applicable.  We  must  also 
decide  whether  the  effective  mass  approxima^^ion  or  the  narrow  band 
limit  is  more  appropriate.  Theoretically,  we  expect  the  bands  to  be 
wider  in  Ti^O^  tlian  in  since  the  overlap  between  ions 

should  be  significantly  greater  than  that  between  the  smaller  ions. 

Furthermore,  the  c>axis  pairs  are  Sfocioser  in  Ti^O^,  which  wouid 
tend  to  increase  the  overlap,  and  hence  the  band  width,  in  the  band, 
which  we  shall  show  in  Chapter  V  is  the  band  of  interest  m  the  case  of 
Ti-O-.  Experimentally,  the  presence  cf  antiferromagnetism  with 

w  .J 

such  a  small  moment  as  0.  Z  pg  per  cation  is  an  indicaticn  of  some 
degree  of  band  width,  as  we  noted  in  the  discussion  following  equations 
(2.20).  Application  of  the  exact  one -dimensional  theory  of 
'‘ectlon  II  C  shows  that  a  band  width  equal  to  three  times 
the  energy  gap  leads  to  a  sublattice  magnetization  of  0,2  pg 
per  cation.  [We  would  expect  In  three  dimensions  that 
a  smaller  band  width  to  band  ^p  ratio  could  suffice  for 
the  same  sublattice  magnetization,  since  the  density  of 
states  In  the  center  of  the  Brlllouln  zone  (low  spin  region) 
Is  far  greater  compared  to  that  at  the  zone  edges  than 
In  the  one-dlmenslonal  case,]  a  band  width  considerably 
wider  than  the  energy  gap  would  seem  to  preclude  use  of 
the  narrow  band  limit. 
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An  additicmal  point  i#  the  Jew  value  of  found  in  equation 
(4,  401,  as  compared  with  the  E  >  0.  6  eV  in  V,0  We  can  take  ^hls 
as  eviaence  that  polaron  effects  are  not  as  important  m  Ti^O^  as  m 
V2O3,  as  we  shall  discuss  Jn  Chapter  VI.  But  this  is  .till  one  more 
indication  that  the  bands  are  relatively  wide  in  Finally, 

we  note  from  our  discussion  of  the  Jahn- Teller  effect  in 
section  II  D  that  a  lack  of  observable  crystalline 
distortion  is  a  wide  band  characteristic,  Ihus  we  con- 

Ci.ude  that  it  is  probably  more  accurate  to  use  the  effective  mass 
approximation  than  the  narrow  band  limit  for 

Unfortunately,  the  presence  of  an  extra  parameter,  m^,  affords 


am  M  greater  degree  of  freedom  in  the  calculation  of  the  transition 
temperature  in  this  case  than  in  the  narrow  band  limit*  Frederikse  [53 
measured  the  effective  mass  of  the  bottom  of  the  lowest  3d  band  in 
TiOg,  and  obtained: 


{4.42: 


We  shall  use  this  value  as  the  best  available  approximation,  noting 
that  the  ratios  of  nearest  neighbor  cation>cation  and  cation-anien 
dlsUnces  to  cationic  radii,  probably  the  best  single  measures  of  over¬ 
lap  available,  are  virtually  the  same  in  TiO^  and  Ti^O^.  Using  (4.  39) 
^go’  solution  of  equations  (3.22)  and  (3.23)  is: 


0.6. 


ThiS  shows  that  Boltzmann  statistics  are  invalid,  ae  might  be  ex¬ 
pected  With  such  a  small  energy  gap.  The  general  Fermi  equations, 
(3.37)  and  (3.46)  must  be  used,  in  the  case  of  Ti<,0,.  thn  solution 
of  (3.37)  and  (3.46),  using  (4,39)  and  (4.42),  is: 


Equations  (4,39)  and  (4,43)  predict  for  the  transition  temperature 


=  490°K  . 

A  comparison  of  (4.44)  with  (4.41)  shows  that  this  value  is  within 
the  experimental  range. 

For  this  solution,  the  jump  In  carrier  concentration  at  the 


(4.44; 


transition  is  much  smaller  than  those  for  the  oxides  of  vanadium. 


Furthermores  the  wider  bands  of  Ti„0  should  give  a  larger  value 
for  the  mobility  in  the  semiconducting  state=  These  points  should 
be  reflected  inthe  larger  conductivity  in  the  semiconducting  region 
and  a  srnaLar  jump  in  conductivity  at  the  transition  in  the  case  of 

As  can  he  seen  in  Fig.  I-l,  this  is  indeed  the  experimental 
Situation.  Connected  with  this,  however,  is  the  fact  that,  for  this 
solution,  a  large  percentage  of  the  available  carriers  have  been  ex¬ 
cited  before  the  transition.  This  probably  means  that  the  effective 
mass  approximation  begins  to  break  down  below  T  ,  and  that  (3.  6)  and 
(3.7)  should  be  modified  by  introducing  the  details  of  the  band  struc- 
fure.  One  possible  way  around  this  is  to  allow  the  effective  mass  in 
35)  to  depend  on  temperature.  If  we  add  to  the  effective  mass  a 
smali  negative  contr  .bution  iinear  in  T,  we  find  that  the  sharpness  cf 
the  transition  disappears,  and  the  transition  to  the  metallic  sta^e 
takes  piace  over  a  small  range  of  temperature.  But  th,s  is  just  the 

peculiar  characteristic  of  the  transi*  on  in  T1,,0,  which  we  noted 

c  3 

previously.  Furthermore,  small  deviations  from  stoichiometry  could 
now  cause  large  variations  tn  and  this  may  be  the  explanation  cf 
the  Wide  range  cf  experimental  values  of  T, ,  as  given  in  (4.41}. 


Chapter  V 

MODELS  FOR  BAND  STRUCTURE 

In  this  chapter,  we  take  the  result 9  of  Chapter  togethe 
what  is  known  about  the  structure  of  the  various  materials,  and  try 
to  piece  together  the  approximate  band  structure  of  each.  Since  some 
critical  information  is  missing  about  every  one  of  these  materials, 
the  band  models  presented  here  contain  some  guesswork.  But  they 
arc  all  consistent  with  both  the  symmetry  and  the  electrical  properties 
of  the  crystals. 

As  we  discuss  in  detail  in  Appendix  B,  little  information  can  be 
gained  from  actual  band  structure  calcu.at  one.  For  the  sake  of  this 
chapter j  we  are  prli  larily  mteresied  .n  the  bands  ai.s^ng  from  the  over¬ 
lap  of  the  3d  wave  functions  of  the  magnetic  ions.  We  shall  use  the 
point  of  View  presented  in  Appendix  B  to  discuss  these  bands.  In  short, 
wt;  begin  With  the  tight  b'nding  approximation,  taking  the  toniC  wave 
functi'ns  to  approximate  the  Wannier  functions  of  the  crystai.  We  a  so 
adopt  the  outlook  of  Anderson  [8] ,  in  which  the  Wannier  func*^'on3 
iCc.aiized  around  the  cations  contain  cont:  fusions  due  to  cverlap  between 
the  cal  ons  and  all  the  Ugands.  These  functions  are  pr  mar^iy  3d  lOnic 
electron  orbitals,  but  they  contain  finite  amplitudes  on  the  surrounding 
anions.  For  example,  for  a  cation  at  the  origin: 

.mt 

R,  are  the  pcs  t^ons  of  ’^he  oxygen  ons.  For  oar  purposes, 
t  s  enablos  us  to  proceed  by  ignoring  the  anions,  assuming  that  they  are 


taken  into  account  in  some  effective  way  ty  using  functions  such  as 
(5*  1),  Functions  of  the  form  (5.  1)  can  have  finite  overlap  even  when 
the  cation  wave  functions  *1^,  ,(r  -  R  )  have  negligible  overlap*,  sin-'e 
they  can  have  finite  amplitude  on  the  same  anion.  We  shall  label  the 
functions  (5.  1)  by  the  type  of  ionic  orbital  of  which  they  mainly  consist. 

We  first  note  that  since  the  ligands  are  large  oxygen  ions  and 
are  doubly  chargedj  the  octahedral  crystal  field  is  undoubtedly  very 
laige,  and  accordingly  we  expect  a  large  separation  of  the  and  e 
bands.  Experimentally,  this  splitting  seems  to  be  of  the  order  of 
1.5  eV  or  »nore.  In  these  materials,  where  the  3d  bands  are  rela¬ 
tively  very  narrow,  this  means  that  we  can  neglect  the  interaction 

between  and  e  bands.  For  all  the  transition  metal  oxides,  each 

8 

cat..on  IS  surrounded  by  an  octahedral  array  cf  0  ligands.  The 

negative  charge  clouds  of  the  e  orbitals  arc  directed  at  these  ncga- 

8 

t..ve  ligands,  and  therefore  the  e  orbitals  have  a  higher  electrostatic 

8 

energy  than  do  the  orbitals,  which  are  directed  between  the  ligands. 

Thus  the  t„  bands  are  lower  than  the  e  bands.  Since  none  of  the 
2g  g 

oxides  of  titanium  and  vanadium  have  more  than  6  3d  electrons,  we 

need  not  consider  the  e  bands  at  all. 

g 

The  bands  contain  6N  states  per  uni*  volume,  which  are 
2g 

not  split  any  further  by  a  pure  cubic  fie’d.  Tetragonal  and  trigonal 
fields  split  the  t.,  band  into  two  sub-bands,  one  with  a  concentration 

2g 

of  2  states  per  cation,  the  other  with  4  stages  per  cation.  Momclinic 
fields  produce  three  sub-bands  with  2  states  per  cation  in  each.  The 
presence  of  antiferromagnetism  or  of  a  uni*  cell  wliich  consists  of 
two  cations  can  bring  about  a  splitting  in  half  of  all  the  bands. 


For  the  corundum  phase  of  ^2^3*  trigonal  field  will  split 

the  t_  band  into  a  t  and  a  t  band.  The  t  band  is  associated  with 
Eg  o  ±  o 

orbitals  directed  al  3  the  c-axis,  the  t  band  with  orbitals  pri- 

i 

manly  in  the  basal  plane.  Below  150°K..  the  crystal  has  monoclinic 
symmetry,  and  the  t^  band  is  split  into  two  sub-bands,  which  wc  shall 
call  the  t^  band  and  the  tj^  band.  The  c-axis  jiairs  form  the  closest 
caUon-cation  distances,  and  these  cations  are  in  octahedral  arrays 
of  anions  which  share  a  common  face.  We  should  therefore  expect 
a  relatively  large  bonding -antibonding  splitting  of  the  t^  band.  The 
monoclinic  distortion,  as  we  have  demonstrated  in  Chapter  IV  leads  td 
the  pairing  of  cations  in  the  basal  plane,  albeit  wi^h  a  somewhat 
larger  cation-cation  distance  than  along  ii  e  c-axis.  Thus  the  bonding¬ 


antibonding  splitting  of  the  t^  and  bands  is  relatively  small.  We 
infer  that  the  t^  bonding  band  is  lowest  for  V^O^,  and  that  it  is  sepa¬ 
rated  from  the  t^  antibonding  band  by  an  energy  Eg2‘ 
band  13  next  lowest,  and  it  is  separated  from  the  t  antibonding  band 
by  an  energy  E  j.  We  assume  Ih-  situation  is  as  given  in  Fig. 

V-Ha).  Since  the  narrow  band  limit  appear f^d  to  give  good  results  for 


V^Oj,  we  have  drawn  the  bands  as  quite  narrow,  and  thus  there  is  no 
overlap.  We  shall  give  an  order-of-magnitude  estimate  of  the  band 


width  below.  Each  of  the  six  sub-bands  contains  one  state  per  cation. 
Thus,  for  with  two  3d  electrons  per  ion,  the  bottom  two 

bands  are  exactly  filled  at  T  =  0,  while  the  top  four  bands  are  com¬ 


pletely  empty.  The  material  is  thus  a  narrow  band  semiconductor  with 


[Of 


SUGGESTED  ( 
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an  energy  gap  We  determined  in  Chapter  IV  that  E  ,  s  0.  10  eV. 

The  theory  presented  there  shows  that  at  =  150^K,  a  phase  trans¬ 
formation  occurs  due  to  the  breaking  up  of  the  basal  plane  pairs, 
changing  the  crystal  structure  to  tae  higher  symmetry  rhomhohedral 
phase. 

As  a  very  rough  estimate,  we  can  use  equations  (2.  18;  to 
determine  the  band  widths  for  ^2^3*  Since  these  equations  were  ob¬ 
tained  from  a  one-dimensional  model,  v,fe  must  consider  these  esti¬ 


mates  as  order-of-magnitude  approximations  at  best.  From  (2,  18): 


(5.2) 


For  the  t^  bands,  the  lattice  parameter  of  interest  is  the  c-axis  pair 

s 

d. stance  of  2.  70  A,  which  yields: 


(E^)^  ~0.03eV.  (5.3) 

o 

For  the  t^  bands,  a  "•2.88  A,  giving: 

(E^)^  ~0.02ey.  (5.4) 


We  have  used  for  the  effective  mass  the  value  calculated  in  sec- 

♦ 

•  On  V%E  for  the  semiconductinc  state  of  V_0, :  2L_  ;;  70. 

2  3m 

From  150*^K  to  600°K,  the  situati  on  is  now  as  shown  in  P'ig. 
V-l(b/.  In  this  temperature  range,  ^  metal,  the  band  being 

1/4  filled.  Thus  is  also  the  temperature  at  which  a  semiconductor- 
to-metal  transition  occurs,  as  is  experimentally  observed.  However, 

^he  theory  of  Chapters  II  and  III  now  can  be  applied  to  the  energy  gap 
E  2  between  the  bonding  and  the  t^  an*^it©nding  bands.  The  application 


of  the  final  results  of  Chapter  III  is  not  immediately  obvious  m  this 

case,  since  there  are  now  electrons  present  in  the  gap  between  the 

two  bands  of  interest.  The  situation  is  analogous  to  that  of  a  normal 

semiconductor  containing  a  large  percentage  of  impurities.  However 

note  that  we  are  concerned  with  the  thermal  excitation  of  electrons 

from  the  t  -  band  to  the  t  -  c  band,  for  it  is  this  excitation  which 
o  o 

reduces  the  bonding  and  thus  tends  to  close  down  th  =  gap.  The  cor 

centration  of  electrons  in  the  t^  -  c  band  is  still  given  by  (5.50), 

and  the  concentration  of  holes  in  the  -  P  band  is  given  by  (3.  51), 

with  E  ~  E  But  now  the  position  of  the  Fermi  level  depends  on 

the  intermediate  partially-filled  band.  If  this  band  is  approximately 

half-way  between  the  t  -  j3  and  the  t  -  o  bands,  E„  =  -E  /2,  and 

o  o  F 

the  final  result,  (3.64),  follows  exactly  as  in  Chapter  III.  If  the 
intermediate  band  is  nearer  one  or  the  other  t^  band,  the  theory  mu'?! 
be  modified  somewnat.  Equation  (2.  should  be  generalized  to: 


E  -E  -pn-dp 
g  go  '  n  "^p  ^ 

since  n  and  pare  no  longer  equal.  The  results  of  Chapter  II  indicate 

that  ~  Pp  ~  "2  *  ^ith  thiS  modification,  it  can  be  shown  that  (3.64) 

remains  approximately  true,  even  when  the  intermediate  band  is 

very  close  to  one  of  the  t  bands. 

o 

Thus  another  transition  is  predicted  at  a  temperature  T^  - 
which  point  ^g->  shrinks  to  zero,  leaving  one  t^^  band  1/3 
filled.  The  behavior,  of  course,  is  still  metallic,  and  no  striking 
change  in  conductivity  will  be  observed  at  T^.  However,  this  transitioi 

V 

changes  the  shape  of  the  Fermi  surface  and  could  account  for  an 


104 


anomalous  resistivity  change  such  as  the  high  temperature  transition 
discussed  in  Chapter  IV.  From  the  temperature  of  this  anomaly, 

=  &00°K»  we  can  estimate; 

E  ,  “  0.4  eV  . 

g2 

Recall  that  the  anomaly  in  resistivity  at  resembled  the 
additional  resistivity  due  to  spin-disorder  scattering  [see  the  dis¬ 
cussion  following  equation  (4.23)]  .  If  is  antiferromagnetic  at 

T  -  0,  then  antiferromagnetisrn  would  contribute  to  the  t^  band  split¬ 
ting.  In  particular,  the  antiparallel  spin  arrangement  of  c-axis 
pairs  is  just  the  magnetic  configuration  found  by  Abrahams  [23]  in 
Ti20^-  We  would  expect  the  same  spin  configuration  in  as  in 

the  structurally  similar  Ti^O^,  so  this  is  a  reasonable  hypothesis. 
Then  the  theory  predicts  the  breakdown  of  long-range  order  (i.e.  the 
Neel  point)  to  occur  at  -  600*^K,  and  the  high  temperature  anomaly 
of  ^2^)^  1®  easily  explained.  We  shall  return  to  this  point  when  we 
consi<^er  spin-disorder  scattering  resistivity  in  Chapter  VI. 


B.  Ti203 


In  the  case  of  Ti-O,,  wlxich  has  corundum  structure  at  all 


temperatures,  the  trigonal  field  results  in  a  splitting  of  the  band 
into  and  t^  sub-bands.  Experimentally  as  well  as  theoretically, 
the  band  turns  out  to  have  lower  energy  than  the  t^  band  [54].  Since 
the  structure  of  Ti^O^  below  500°K  is  the  same  as  that  of  between 


*  o  o 

150  K  and  600  K,  we  might  expect  similar  band  schemes  in  these 
temperature  ranges.  However,  as  we  demonstrated  in  Chapder  IV, 
the  bands  in  Ti^O^  are  much  wider  than  those  in  V^O^,  The  same 


estimate  which  led  to  (5,  3)  and  (S.'i)  can  be  used  with  the  values 

» 

=  25j  a  =  2.  59  A  to  calculate  the  band  widths  in 


(Ej^)^  ~  0.  10  eV. 
o 


5) 


We  know  the  energy  gap  due  to  the  antiferromagnetism  of  c-axis 
pairs  is  0.  06  eV.  Thus  the  ratio  of  band  w’dth  to  band  gap  is  drasti¬ 
cally  different  in  Ti^O^  and  V^O,.  Note  that  the  cation  separation  in 
the  basal  plane,  which  determines  the  bandwudth  in  (5.2)  is  con¬ 
siderably  larger  in  Ti^O^  (3.02  A)  than  in  (2.88  A),  whereas  the 

separation  of  the  c-axf  s  pairs  is  in  the  opposite  order  for  the  two 

materials.  Thus  (5.2)  estimates  the  t  bandwidth  for  Ti-C-  as: 

3 

(Ej^)t^  ~  0.  06  eV  .  (5.6) 


The  band  scheme  for  Ti  O.  below  T  is  given  in  Fig.  V-2(a). 

Since  there  is  only  one  3d  electron  presen*  per  cation,  the  ”  P  banc  . 

whcch  corresponds  to  electrons  primari.y  with  spins  on  their  own 

5ubla*^tice  (analogous  to  the  bonding  bands  ..n  V.,0-),  is  exactly  filled 

d  5 

while  all  higher  t^^  bands  are  compiete-y  empty.  In  this  temperature 

range,  Ti-O-  is  therefore  a  semiconductor,  with  a  gap  of  E  =  0.06  eV. 
^  ^  S 

brought  about  by  c-axis  antiferromagnetism.  The  theory  then  pre¬ 
dicts  a  semiccnductor-to-metal  transition  at  about  T  =  500°K.  Above 

o 

500  K,  the  band  situation  is  given  in  Fig.  V-2(b).  In  this  range,  the 

bands  are  l/6  filled,  and  the  material  is  metallic. 

2g 

The  details  of  the  band  structure  presented  here  offer  an 
alternative  explanation  for  the  spread-out  nature  of  the  transition  in 
Ti^O^.  This  is  the  effect  on  the  transition  of  the  probable  overlap  of 


the  I  band  with  the  t  -  e  (conduction,  band.  Since  the  transition 
a  o 

in  Ti,jO  does  not  occur  until  the  conduction  band  is  relatively  w-eil 
thermally  populated,  the  significant  change  in  density-of-states  due 
to  the  t^  band  will  clearly  be  important. 

It  should  be  noted  that  the  main  difference  between  Ti^O^  and 
IS  the  number  of  3d  electrons  present  per  calion. 
be  semiconducting  up  until  the  high-temperature  transition  were  it  not 
for  the  additional  3d  electron  on  each  ion.  Then  the  conduct  v  ’y 

would  greatly  resemble  that  of  Ti^O^,  despite  the  smaller  gap  and 
Wider  bands  of  the  latter.  Further  evidence  for  this  picture  will  bo 
given  in  Chapter  VT. 

C.  VO  and  TiO 

Beicw  L26*^K,  VO  undergoes  a  d.  stort  on  to  orthorhombic 

symmetry,  so  that  the  t.,  bonding  and  antibonding  bands  are  each 

*•8 

sp.it  into  three  sub-bands,  which  we  cai.  t  ,  t  ,  t  .  The  exact  low 

X  y  z 

*emperature  structure  is  as  yet  unknown,  so  that  we  cannot  make 
any  statement  about  the  relative  positions  of  these  sub-bands.  We 
arb  trarily  take  ^he  orthorhombic  lattice  parameters  to  be  smallest 
along  the  x-axiS,  largest  along  the  z-ax^s.  In  the  case  of  VO,  we 
expect  the  bands  to  be  wider  than  in  ^2^3’  probably  not  so  wide 
as  in  Ti,,0^,  since  the  V  ic  n  is  smaller  in  extent  than  is  the  Ti^  ^ 
iOn.  Whether  or  not  there  is  overlap  between  the  sub-bands  is 
r relevant  in  VO,  as  long  a*,  there  is  a  real  gap  between  the  t  bonding 
and  the  antibonding  bands-  From  Ciiperiment,  we  know  that  this 
gap  exists  and  is  about  0.  10  eV.  The  band  structure  of  VO  below 


126°K  is  sketched  in  Fig.  V-3(a).  There  are  three  3d  electrons  per 
cation,  and  so  the  lower  three  sub-bands  are  just  filled,  while  the 
upper  three  sub-bands  are  completely  empty.  Thus  VO  is  a  semi¬ 
conductor,  whose  gap  of  0.  10  eV  is  brought  about  by  a  crystalline 
distortion.  The  theory  predicts  a  transition  at  T^  =  126°K  to  the 
situation  as  shown  In  Fig.  V-3(b).  Above  T^,  VO  has  a  l/2  filled 

t.  band  and  is  metallic. 

2g 

Note  that  TiO  retains  r'.ck  salt  structure  down  to  at  least 
4°K.  The  presence  of  antiferromagnetism  has  not  been  established, 
but  since  TlO  co.ntalns  only  two  3d  electrons  per  Ti^^  ion,  it  must 
exhibit  metallic  behavior  at  all  temperatures.  As  can  be  seen  from 
the  results  of  Morin's  work  [20]  (see  Fig.  I-l),  this  is  the  case  down 
to  at  least  1.  5°K. 

D.  VO^ 

The  rutile  phase  of  VO^  has  tetragonal  symmetry,  and  splits 
the  t-  band  into  a  1  sub-band  with  4  Jtates  per  cation  and  a  t  band 
with  2  states  per  cation.  Since  c  is  considerably  smaller  than  a 
in  rutile  VO-,  we  expect  the  t  band  to  be  well  below  the  t  band.  At 
low  temperatures,  the  crystal  undergoes  a  distortion  to  monoclinic 
symmetry,  the  c-axis  cations  of  the  rutile  structure  pairing  and 
puckering  somewhat.  This  monoclinic  distortion  has  two  effects. 

The  t  band  is  split  into  sub-bands,  which  we  shall  call  the  t  ,  and 
the  t^2  bands.  But,  more  important,  the  t^  band  is  split  into  bonding 
and  antibonding  sub-bands  by  the  pairing  along  the  c-axis,  which  is 
a  classic  example  of  the  theory  of  crystalline  distortion  given  in 


section  II  D.  The  distortion  parameter  €  is  rather  large,  which  tends 

to  bring  about  a  large  energy  gap,  estimated  in  Chapter  FV  to  be 

0.3  eV.  The  band  structure  below  T  =  340°K  is  outline'!  in  Fig. 

o 

V-4(a).  Since  there  is  just  one  3d  electron  per  ion,  the  ~  P 

band  is  completely  filled,  and  all  higher  bands  are  empty.  The  theory 
predicts  a  transformation  at  T^  to  rutile  structure.  Above  340'^K, 
the  band  structure  is  given  in  Fig.  V-4(b).  It  can  be  seen  that  the 
band  is  l/Z  filled,  resulting  in  a  scmiconductor-to-metal  transition 
at  340°K,  as  is  observed. 
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FIG.  ¥-4  SUGGESTED  BAND  SCHEME  FOR  VO 


Chapter  VI 


DISCUSSION 


A.  The  Metallic  State  of  V^O, 


Primarily  because  of  the  exhaustive  work  of  Fe^nleib  [2l], 
much  more  is  known  about  than  about  any  of  the  other  materials 

which  exhibit  semiconductor -to-metal  transitions.  Therefore,  much 
can  be  gained  by  considering  hi  more  deta‘1. 

In  cdew  of  the  recent  accurate  measurements  of  the  lattice 
parameter  changes  at  the  transition  [37]  ,  we  shall  recalculate  the 
latent  heat  of  transformation  and  the  Fermi  energy  from  Feinleib's 
pressure  measurements  [21]  .  The  Clausius-Clapey ron  equation  gives: 


L  =  ^  X  V  X  T^  / 


3 

.  035  X  30.  8  X  152°K 

T" 


3. 86  X  10 


•  9  K  cm‘ 
Syne" 


=  1020 


cal  A  A  eV 


moie 


=  0.044 


molecule 


(6. 


We  can  use  the  value  of  latent  heat  to  approximate  the  Fermi  energy, 
as  was  done  by  Feinleib  [55]  .  u  assume  no  spin  changes  occur  at 
the  transition.  The  difference  in  free  energy  between  the  semi¬ 
conducting  and  metallic  states  at  the  transition  point  can  then  be 
expressed: 


0  AG  -  AE  -  T  {AS  ,  ,  .  +  AS,  ) 

o  electronic  lattice 

2  .2  'I 

,  TT  k  T 

«  AE  -  i - -  kT  . 

I  o 

F 


{6.2  5 


According  to  the  model  we  presented  in  Chapter  V,  two  3d  electrons 
per  molecule  (one  3d  electron  per  cation)  contribute  to  th»'» 

metallic  band.  Setting  AE  in  (6.2)  equal  to  half  the  value  of  the 
latent  heat  per  molecule  found  in  (6.  1),  we  obtain: 

2  2  2 

tr  k 

tr  -  O 

L  -  2k  T 

o 

ss  0.  09  eV  .  (6.  3 

This  is  a  reasonable  magnitude  for  Ep,  consistent  w.dh  the 

orders  of  magnitude  assumed  in  Chapter  V  for  the  band  parameters 

in  With  this  value  for  Ep,  the  band  width  in  *-he  metal, *c  sta’e 

can  be  estimated  as  E.  ~  0.  18  eV,  Note  that  E„  =  7  kT  ,  so  that 

b  F  o 

usual  statistical  approximations  are  valid. 

We  can  now  use  (6.  3)  to  calculate  the  effective  mass  of  eiec  - 
Irons  in  the  metallic  state.  In  the  effects. ve  mass  approximation  for 
metals  [56]: 


E 


F 


,,  2...2/3 
- 5  (3tr  N)  '  . 

2m 


c  ,  4 


Hence: 

m  *"  2m  Ep 

=  50  .  ,6.  5  > 
This  order  of  magnitude,  effective  mass  is  roughly  who.t  we  would  ex¬ 
pect  for  V^O^.  We  can  approximate  the  effective  mass  in  the 
metallic  state  in  a"  cher  way,  to  check  (6  5).  Feinleib  [57|  has  estimaie 
the  plas  .na  frequency  from  optical  data  oji  the  metallic  state  of  \ 
ujp  =  1. 2  eV  . 


ci  S  I 


Bui  irom  the  usual  expression  for  the  plasma  frequencry  [58]: 


■*  5 

m  4ttN 


This  is  in  agreement  with  the  approximation  (6.  5). 

We  can  now  use  this  value  for  the  effective  mass  to  calculat« 
the  relaxation  time  in  the  metallic  state.  From  (4.22),  we  get  for 
the  conductivity  just  above  the  transition: 


=  4.6  X  10^  ohm'^  cm"^  . 


Hence,  the  mobility  is: 


(6.6/ 


2 

~  0.  72  cm  /volt  -sec. 

The  relaxation  time,  in  view  of  (6.5),  is; 


r  =  X 

e  m 


=  2.0x10  sec. 


(6.  7 1 


(6.8) 


This  IS  a  typical  magnitude  of  relaxation  time  for  a  good  metal.  It 
IS  clear  that  the  difference  in  mobility  between  metallic  and 

metals  such  as  copper  is  solely  due  to  the  large  effective  mass  in 
^2^3*  means  that  the  carriers  are  moving  much  more  slowly 

than  they  are  in  copper.  Thus,  although  the  average  time  between 
collisions,  (6.8),  is  the  same  for  and  copper,  the  average 

distance  between  collisions,  or  mean  free  path,  is  much  lower  for 


1.5 


V  0-.  Using  (6.  3),  we  evaluate  the  velocity  of  carriers  at  the 
Fermi  surface  as; 

Up  -  ^ 

rn 

8  / 

=  0.  03  X  10  cm/sec  .  (t ,  9) 

From  (6»  8)  and  (6.  9),  w-c  find  for  the  mean  free  path; 

A  =  T  Up 

-  6  A  .  (6.  10, 

Therefore,  the  mean  free  path  is  about  2  ~  lattice  constants  in 

V^Oj,  as  opposed  to  about  100  lattice  constants  in  copper.  All  these 

values  are  consistent  with  the  picture  of  a  partially  filled  narrow 

band  resulting  in  the  observed  metallic  conductivity.  However,  we 

are  approaching  a  borderline  situation,  since  if  the  mean  free  path 

gets  much  below  one  lattice  constant,  the  carriers  are  essentially 

localized  and  the  bands  can  be  considered  to  be  washed  out. 

As  the  temperature  is  increased  above  T^,  kT  increases, 

until  at  600*^K,  it  is  almost  lialf  the  Fermi  energy,  as  given  by 

(6.  3).  At  this  point,  the  statistical  variations  in  Fermi  energy 

begin  to  become  important,  and  we  are  entering  another  borderline 

region.  But  here  the  high  temperature  transition  helps  retain  metalTc 

behavior.  For  at  =  600°K,  the  B  and  the  t^-  o  bands  merge, 

and  above  this  temperature,  there  is  essentially  one  relatively  wide 

t-  band,  1/3  filled.  This  means  that  the  number  of  free  electrons 

2g 

contributing  to  the  metallic  conductivity  has  doubled.  Assuming  the 


effective  mass  remains  the  rams,  the  doubling  of  N  serves  to  in-  . 
crease  the  Fermi  energy  to: 

Ep  =  0,  14  eV.  (6. 

We  conclude  that  the  experimental  data  on  above  150°K  can  be 

satisfactorily  explained  by  the  model  presented  in  Ciiaptere  IV  and 
V,  However,  the  bands  can  probably  not  get  much  narrower  and 
sti)l  maintain  metallic  conductivity-  We  shall  add  one  more  cupport- 
ing  argument  in  section  B,  where  we  consider  spin-d”sorder  scat¬ 
tering. 

B.  Spin-Disorder  Scattering 

The  theory  of  spin-disorder  scattering  has  been  given  by 

De  Gennes  and  Friedei  [46]  .  There  are  two  main  effects  of  this 

type  of  scattering,  which  we  completely  ignored  in  Chapter  III. 

Firstly,  there  is  a  contribution  to  the  resistivity,  which  is  small 

when  the  spins  are  highly  ordered  {i.e.  below  Tj^),  but  adds  a  term 

independent  of  temperature  after  long  range  order  disappears. 

Secondly,  there  is  the  effect  of  the  broadening  of  the  bands  due  to 

spin  disorder.  This  broadening  will  affect  the  relations  developed 

between  E  and  T  , 
go  o 

If  our  hypothesis  about  presented  in  Chapter  V  is 

correct,  we  can  calculate  the  spin-disorder  resistivity  above  Tj^  - 
600°K,  and  compare  it  to  the  experimental  jump  in  resistivity  at 

measured  by  Feinleib  [2.1]  [see  discussion  below  equation  (4.E4)]  . 
We  first  apply  the  theory  of  sections  II  C  and  III  B  to  esti¬ 
mate  the  antiferromagnetic  energy  gap  from  Tj^.  Since  we  expect 


narrow  band  theory  to  be  applicable,  equation  (3.  ?>9)  is  the  appropriate 


relation  for  E  ,  and  we  find; 
go 

E  “  0.  4  eV 

gc> 


.2' 


as  we  already  estirnateo'  in  Chapter  V  (where  we  referred  to  this  gap 
as  ^g2^*  gsp  is  due  to  the  f  oheresit  scattering  arising 

from  the  exchange  potential.  But  it  is  the  incoherent  scattering 
caused  by  tbiis  aame  potential  which  gives  rise  to  the  spin-disorder 
resistivity.  Using  the  Born  approximation  and  the  extreme  simpj.t  ra¬ 
tions  of  spherical  energy  surfaces  and  quasi-free  electrons,  ♦^ho 
spin-disorder  resistivity  above  Tj^  can  easily  be  calculated  in  terms 
of  the  energy  gap. 

Under  the  assumption  of  quasi-free  electrons,  the  energy 
gap  can  be  expressed  as  twice  the  Fourier  component  of  the  exchange 
potential; 


;6.  3 


where  k  -  k'  is  tw'i.ce  the  Fermi  momen*^um,  k^.  For  the  case  of 

V.,0  ,  using  the  value  for  E  calculated  in  equation  (6.  12): 

3  go 


Vjj,  =  0.  2  ev  . 
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For  spherical  energy  surfaces,  the  resistivity  can  be  expressed  [59] 


Pso  =  - ^2-^  I  f  dn  dfi'  i  ^  (I  -  cos  6)  . 

16it  fie  V  N 


Performing  the  integrations,  (6.  15)  can  be  written: 


^so  "  T 


3»  ,  "■''"Kg-' 

-  .  (4  ) - — _ 

°  ^  fie  N  E 


c  ^  ^  m 


But  for  the  high  temperature  phaee  of  using  (6.5),  (6.  ll)^ 

and  {6.  I4)j,  equation  {ht  16)  gives  for  the  *pin*di#order  resistivity: 


118 


^SO 


IZ  X  i0“^ 


ohm -cm. 


(6.  17) 


This  is  just  the  value  oX  measured  by  Feinlcib  for  the  anomalous 

resistivity  increase  in  above  600’^K,  Although  the  values  of 

Fermi  energy  and  of  the  effective  mass  used  in  (6, 16)  are  just  rough 
approximations,  they  can  both  be  arrived  at  by  two  different  methf"!* 
and  ws  thus  do  not  expect  too  much  variation.  Therefoie,  the  agree¬ 
ment  between  theory  and  experiment  can  be  considered  quite  good. 
However,  we  cannot  prove  the  existence  of  antiferromagnetism  in 
by  such  arguments,  since  the  agreement  could  be  accidental. 

If  ^2^3  were  non-magnetic  at  all  temperatures,  there  would 
be  a  contribution  to  resistivity  from  the  spin-disorder  scattering. 

In  order  to  estimate  this  contribution  from  {6.  15),  we  must  use  the 

ijt  o  o 

appropriate  values  of  Ep,  N,  and  m  for  the  range  150  K  <  T  <  600  K, 

and  we  tiust  approximate  in  some  way.  One  possibility^  is  to 

ass'.me  that  the  exchange  energy  is  the  same  for  V..O  as  for  Ti^O.. 

Then  we  can  use  the  value  of  E  in  to  obtain* 

go  2  3 


0.03 


ev 


Evaluating  {b.  15)  in  this  case  gives: 

*4 

PgQ  “  i.4  K  10  ohm -cm. 


(6.  18) 


But  the  temperature  independent  part  of  the  metallic  resistivity  as 
measured  by  Feinleib  j  2l|  and  given  in  (4.  23)  is: 
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Thus  (6.  18)  shows  that  spin-disorder  atteTing  could  account  ii-t 

almost  n/3  of  the  temperature  independent  part  of  the  resist  Vitv 

above  T  .  The  remaming  l/3  could  be  due  to  the  effects  of  im- 
o 

purities,  imperfections,  or,  as  we  shall  see  in  section  C,  polaron 
formation.  Clearly,  spin-disorder  resistivit'/,  being  of  the  order 
given  by  {6.  18),  is  completely  negligible  in  the  semiconducting 
region. 

It  i®  also  irr.portant  to  estimate  the  amount  of  broadening 
of  the  bands  brought  about  by  spin-disorder.  Us^ng  the  same  assump¬ 
tions  as  in  equation  {6.  13),  the  perturbing  potential  due  to  spin 
disorder  can  be  expressed  as  half  the  change  in  energy 


/V  ;>  =  1  (E  -  E  )  . 
\  pert/  2  go  g 


i‘'  = 


Using  equation  (2.42)  for  the  change  m  energy  gap,  we  obta  n  (?Gm 
9): 
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Thus,  the  mean  square  deviation  in  energy  due  !o  spin  disorder 


(AE)  ”  E 


‘o.  2 


Equaiiun  (b.2i)  shows  that  the  etfec*  cf  sp’n-disorder  broaden'ng  s 
similar  to  the  modifications  brought  about  by  a  Gaussian  spread 
around  delta  function  bands  in  the  na?  row  band  limit,  except  that 
the  spread  is  no  longer  cons’^ant,  but  depends  or.  the  amount  of  spin 
disorder,  and  thus  on  n.  We  have  already  worked  out  the  case  '•■f  a 
Gaussian  spread  in  deta  l  in  section  IK  B,  Using  (6.2’?  as  the  vaiuc 


for  X  n  equation  (3*  70), 


wc  obtain: 


(6,22) 


y|U4x) 

X  -  e  ^  e 
Solving  (6,22)  give;s; 

E 

~  f»6  .  (6.23) 

o 

Thus*  spin-disorder  broadening  lowers  the  transition  temperature 
somewhat,  as  wc  would  expect. 


C.  Polar  on  Effects 

A  second  effect  which  has  not  been  considered  thus  far  is 
polaron  formation.  It  is  clear  from  the  extremely  narrow  band 
widths  that  only  '*8mali”  polarons  are  involved  in  these  materials. 
Holstein  [10]  has  calculated  the  conductivity  in  the  two  cases  where 
transport  occur  s  in  a  polaron  band  and  where  an  electron  hops  to  a 
neighboring  site.  Cond  ction  in  a  polaron  band,  which  should  domi¬ 
nate  at  very  low  temperatures,  is  characterized  by  an  exponential 
decrease  with  increasing  temperature.  This  is  certainly  not  the 
situation  in  the  oxidej^  of  vanadium  and  titanium.  For  thermally 
activated  hopping,  the  experimentally  observed  temperatu-  e  de¬ 
pendence  is  obtained.  The  activation  energy  in  the  temperature 
range  where  polaron  hopping  is  the  dominant  mode  of  conduction 
is  [  10]  : 


K 

w 


dk(l  -  cos  k) 
- —2 - 


(6.24) 


where  K  is  a  constant  depending  on  the  mat  s  of  the  ions  and  the 
strength  of  the  electron-phonon  interaction,  and  are  the  optical 


\1 


phonon  frequencies.  Assuming  a  narrow  -band  vibrational  spectrum 
(Einstein  model),  which  is  a  good  approximation  for  optical  phonons: 


w,  =  w 
k  o 


we  obtain  from  (6.  24): 
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(6.25;- 


Although  the  mean  phonon  frequency  can  be  expected  to  increase 
somewhat  with  pressure,  the  likelihood  is  smal.'  that  it  would  change 
sufficiently  to  account  for  the  decrease  in  activation  energy  by  a  factor 
of  2  at  pressures  of  20  kxlobars,  as  measured  by  Feinleib  [21]  for 
and  by  Austin  [43]  for  and  VO,  Thus  it  is  unlikely  that 

thermally  activated  hopping  of  polarons  Is  the  major  contributor  to  the 
(  onductivity . 

However,  we  can  estimate  the  con:  r  button  of  polaron  effects  to 
the  pressure  cocftici#;nt  of  the  total  act-vat  on  energy  from  (6.25). 

Fic] II avion  (6.2  5)  implies: 


d  In  (E  .  )  ,  2  d  In  u> 
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But  Gruneisen's  relation  [60]  shows; 
d  In  cj 

u  c 

■JITV  ' 


(6.26) 


p  K  C, 
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where  a  is  the  thermal  expansion  coefficient,  p  is  the  density,  x  is  the 
isothermal  compressibility,  and  c^  is  the  specific  heat  at  constant 
vol'ime.  Solving  (6.27)  for  the  pressure,  coefficient  of 


d  In  u) 

“dT^ 
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(b.28) 


nz 


The  value*  appropriate  to  V,0,  are  a  =  40  X  10"^  c  -  6.1  K 

6  *  ^  V 

fiO  dyne-cm/g«'  K,  and  p  ^  S.O  g/cm^.  Subftituting  in  {6.Z9)t  we  find; 
d  In 

*  1,2X10“^  bar*^  (6.29) 

Combining  16.26)  and  (6.29): 


— 


.2,4  X  10’^  bar’“ 


(6.301 


Comparing  (6.30)  to  (4.21)  and  (4.22),  we  see  tixat  the  contribu¬ 
tion  to  the  pressure  coefficient  from  polaron  effects  is  an  order  of  mag¬ 
nitude  smaller  than  the  total  observed  value.  We  were  therefore  justified 
in  neglecting  polaron  effects  in  the  case  of  V^O^.  A  similar  comparison 
for  VO  shows  that  the  polaron  contribution  is  a  negligible  part  of  (4.32). 
However,  for  ^02*  this  is  not  the  case.  As  we  pointed  out  in  section 
IV  C,  the  polaron  contribution  to  |4.36)  *s  quite  significant.  In  (4.36), 

We  found; 


dE 

-5^  *  -5.0  Xl0‘^eVbar’"  . 


(6.31) 


Us*ng  C3.5?l,  (4.35),  and  our  estimate  of  E  -0.3  eV,  we  conclude; 


i  ~  -  2,  1  X  10“^  eV  bar’"  . 

Putting  (6.  31)  and  (6.  32)  into  (4.  37); 
dE 

'^[P  ^  -  2. 9  X  10  eV  bar"*  . 


Since  E^  —  0.3  eV,  (6.33)  shows; 


(6.32) 


(6.33) 


d  In  £ 


IF 


-  1.0  X  10"^  lar' 


(6.34) 


Evaluating  (6. 28)  for  VO^#  we  obtain: 
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A  comparison  of  (6.34)  and  (6.  3S)  shows  that  polaron  effects  could  be 
responsible  entirely  for  the  obtained  by  Neumann  et  al  [51]  in  VO,. 

We  expected  VO.^  to  be  the  material  with  the  narrowest  bands  of  those 
we  considered.  Consequently,  we  would  expect  the  largest  polaron 
effects  in  VO^,  It  is  significant  that  the  value  of  ~  0, 3  eV  for 
was  indeed  largest.  For  V^O^,  which  should  have  slightly  wider  bands 
than  VO_  has,  E  •“  0.  07  eV  to  0,  15  eV.  For  VO,  E  0.  1  eV,  similar 
to  ^2^3*  Finally,  for  Ti^O^,  where  the  effective  mass  approximation 
was  found  to  give  better  resxdts  than  the  narrov/  band  limit,  (4,40)  shows 
that  £^"0,01  eV.  It  seems  reasonable  to  conclude  that  polaron  effects 
are  responsible  for  at  least  a  large  pirt  of  in  these  materials. 

It  is  worth  noting  that  the  decrease  in  resistivity  as  a  function 
of  pressure  in  NiO,  CoO,  CuO,  and  MnG,  as  measured  by  Minomura 
and  Drickamer  [61] ,  are  considerably  less  sharp  than  they  are  in 
and  VO.  For  the  former  materials,  the  calculated  changes  in  with 
pressure  by  means  of  (6.28)  can  completely  account  for  the  experi¬ 
mental  pressure  coefficients  of  E^.  Thus,  for  these  transition  metal 
oxides,  the  model  of  conduction  by  means  of  polaron  hopping  may  indeed 
apply.  Since  we  would  expect  the  bands  to  be  still  narrower  in  these 
materials  than  in  which  was,  as  we  have  seen,  almost  a  border¬ 

line  case,  it  is  not  surprising  that  they  are  non-conducting  at  all 
tempe  ratures. 

Polaron  effects  can  also  manifest  themselves  in  the  metallic 
region.  For  example,  if  there  were  a  weak  activation  energ'y,  so  that 


1E4 


■esiitlvity  were  of  the  form: 


P  *  T  e 


FT 


T  + 


then  the  activation  energy  would  contribute  to  the  temperature  inde¬ 
pendent  term  in  (4.23),  A*  wo  thowed  in  (6.  18),  spin-disorder  scat¬ 
tering  could  account  for  at  roost  2/3  of  this  term.  It  is  possible  that 
polaron  effects  are  responsible  for  the  remaining  l/3, 

D.  Effects  of  Non-ftoichiometry 

It  is  also  important  to  consider  the  effects  of  deviations  from 
stoichiometry,  both  ia  the  form  of  impurities  and  imperfections,  par¬ 
ticularly  in  the  materials  under  consideration,  wluch  are  difficult  to 
obtain  in  their  pure  form.  For  simplicity,  we  shall  consider  the  case 
where  there  is  a  given  number,  N^,  of  donor  levels  per  unit  volume. 
The  case  where  acceptor  levels  are  present  follows  analogously.  Let 

n.  be  the  concentration  of  intrinsic  carriers  excited  into  the  conduction 

1 

band.  At  extremely  low  temperatures,  it  is  possible  to  have  electrons 
partly  trapped  on  the  donor  sites;  however,  such  temperatures  are  not 
reached  in  practice,  so  we  shall  ignore  this  region.  At  higher  tem¬ 
peratures,  the  donors  become  completely  ioniaed,  and  saturated  ex¬ 
trinsic  conductivity  dominates.  Here  the  concentration  of  carriers  in 
the  conduction  band  is  given  by; 


,  -E  /2kT 
..  ,  e  * 

"  '  '’d  *  - - - 


(6.36) 


In  this  temperature  region,  the  conductivity  will  not  vary  •tr'mgly  with 
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temperature,  but  will  vary  with  composition.  As  the  temperature  is 
increased,  becomes  comparable  with  and  the  general  expres¬ 
sion  for  n  must  be  used; 

4n^  1/2 

n=|N  [l  +  n^—i)  ].  (6,37) 

^d 

At  still  higher  temperatures,  intrinsic  conductivity  dominates,  and  n 
can  be  approximated  by: 

-E  /2kT  j 

n  =  N  e  ^  ^  .  {6. 38) 

Ii*  this  region,  conductivity  v/^U  exhibit  an  exponential  temperature  de¬ 
pendence,  but  be  nearly  independent  of  composition.  In  order  to 
determine  the  extent  of  the  intrinsic  region,  the  experiments  of 
MacMillan  [42]  are  particularly  useful.  MacMillan  obtained  V  O. 
samples  by  reducing  V20^  in  hydrogen  at  different  temperatures.  He 
was  thus  able  to  vary  the  composition  from  to  V^O^  In 

measurements  down  to  100°K,  it  was  found  that  conductivity,  activa¬ 
tion  energy,  and  transition  temperature  were  all  independent  of  com- 
positior  This  is  strong  evidence  that  conductivity  is  primarily 
intrinsic  in  this  temperature  range,  and  that  equation  (6.38)  should  be 
used  to  evaluate  n. 

The  effect  of  donor  and  acceptor  levels  on  the  theory  can  now 

be  determined.  Equation  (3.62)  must  be  modified  in  the  following  way: 

y(2  -  4x)  -1  , 

+1]  +  2 

where  =  N^/N.  We  are  dealing  with  concentrations  of  impurities 
such  that  «  N,  and  thus  «  1.  Solution  of  equation  (6,  39)  shows 


IZb 


that  the  traiiSitioa  is  moved  only  slightly:  x  is  decreased  and  y  is 

'  o  o 

increased.  The  fractional  Increase  in  y  Is  Z.B  x..  or  2,8  N,/N. 

o  d  d- 

Ihua  donor  levels  decrease  the  transition  temperature,  A  similar  cal¬ 
culation  performed  in  the  next  section,  shows  that  acceptor 
levels  also  tend  to  decrease  in  the  same  proportion  as  do  donors. 

Since  the  experimental  results  quoted  in  Chapter  IV  seem  to 
be  in  the  intrinsic  conductivity  region,  the  observed  activation  energies 


must  still  contain  contributions  of  1/2  E  .  Then,  since  the  s^olution 

S 

of  (6,39)  maintains  the  proportionality  of  to  T^,  equation  (3,  57) 
remains  valid  for  small  amounts  of  non-stoichic  etry. 


A  major  effect  of  small  deviations  from  stoichiometry  is  the 
contribution  of  the  corresponding  scattering  to  the  mobility.  The 
mobility  due  to  scattering  by  ionised  impurities  is  given  by  the  Conwell- 
Weisskopf  formula  [62]  t 


(6.40) 


where  K  and  Tj,  are  constants  which  depend  on  the  number  of  impurities 
and  the  dielectric  constant  of  the  material.  Thus,  if  impurity  scat¬ 
tering  dominates,  (6,  40)  shows  that  the  contribution  to  the  observed 
activation  energy  is: 

<=A>i  “  I ''■I' *  r— TT-r?  • 

jin  (i^  )  I 


Since  the  number  of  impurities  is  small.  Tj  T,  and  the  first 
term  of  (6.41)  is  the  important  one.  The  activation  energy  due  to 


ionizcd  impurities  is  thus  small  and  somev^’hal  temperature  dependents 


Since  there  is  no  pressure  def  endence  of  the  actlva^iOn  energy  given 

by  (6.41),  ionized  impurity  scattering  cannot  account  for  the  pressure 

variation  of  the  observed  activation  energies  in  V  O*,  VO,  or  VO*. 

i  i  £. 

Furthermore,  we  were  justified  in  neglecting  in  Chapter  IV  the 
pressure  dependence  of  the  contributions  to  from  this  type  of 
scattering. 

Scattering  by  neutral  impurities  or  by  dislocations  are  tem¬ 
perature  independent,  and  are  small  except  at  very  low  tempera¬ 
tures  [63] .  Thus  there  is  no  contribution  to  from  these  sources. 

As  we  already  pointed  out,  the  transition  temperatures  observed 
in  the  oxides  of  vanadium  do  not  seem  to  change  with  small  deviations 

from  stoichiometry.  However,  the  measured  T  i.n  Ti.,0  does  seem 

o  ^3 

te  vary,  although  the  fact  that  the  transition  is  not  very  sharp  makes 
it  difficult  to  decide  on  an  exact  value  for  in  this  case.  It  appears 
that  the  effects  of  non-stoichiorrietr>'^  can  be  important  only  in  Ti^O,. 
Some  supporting  evidence  for  such  a  conclusion  comes  from  the  work 
of  Yahia  and  Frederikse  [22]  ,  who  found  that  the  Hall  mobility  in 
Ti.,0-  obeys  a  Conwell- Weisskopf  law  below-  ’56?K.  Between  200°K 
and  450°K,  polar  scattering  appears  to  predominate.  All  this  is 


consistent  with  the  idea  that  the  relatively  wide  bands  of  Ti^O^  cause 
polaron  effects  to  be  negligible,  whereas  in  the  narrow  bands  of  ^^^3* 
VO,  and  VO^,  polaron  effects  are  quite  important. 


E,  The  V  O  -Ti  O  System 
J  iL  5 

The  models  for  V^O,  and  Ti^O^  presented  in  Chapters  TV  and 
V  taken  together  with  the  discussion  of  non- stoichiometry  in  section 
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VI- D  enable  n#  to  make  prediction®  abovt  the  electrical  conductivity 
properties  of  the  ayste;  .  Subatitution  of  for  =  + 

introduces  koXea  in  the  valence  band,  since  a  ion  has  only  one 

3d  electron  and  consequently  acts  like  an  acceptor* 

Let  bo  the  fractional  amount  of  Ti  introduced  into  a 
crystal  of  V^O^*  Thus  =  N^/N,  where  is  the  concentration  of 
acceptors,  la  general,  we  know; 
p  -  n  = 

2 

np  =  n. 

so  that: 

,  4„i' 

In  the  narrow  band  limit,  (6.42)  becomes; 


i  r  1  /  ,  4  -^:y2kT  1/2 

x=^x^[l  +  (l+-^e  «  )  ] 


(6.43) 


where  x  £  p/N  is  the  fractional  amount  of  holes  in  the  valence  band. 
For  small  x  «  I,  equation  (6. 43)  becomes; 


X  ^  e 


-E  /2kT  , 

S  *ix 

z  a 


(6.44) 


wWch  is  the  relation  for  acceptors  analogous  to  (6.  38).  Introduction 
of  (2. into  (6.  44)  gives; 


-y(l-4x)  1 


(6.45) 


where  y  =  E^^2kT  as  before.  The  solution  of  (6.  45)  shows  that  a 
small  percentage  of  acceptors  lowers  the  transition  temperature  b>'  a 


small  amount.  For  example,  when  1%  Ti  is  introduced  into  V^O^, 

the  transition  temperature  determined  by  (6.45)  decreases  by  3%. 

The  approximation  leading  to  (6.  45)  is  va^id  in  the  case  of  up 

until  2^c  Ti^^^  conc^  itration. 

As  X  gets  large,  we  must  return  to  the  general  relation, 

■'6.43).  When  N  »  n,,  (6.43)  can  be  written: 
a  1 


X  -  X  +  e 

a  X 
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,  -E  /kT 


(6.46) 


Introduction  of  (2.  1)  into  (6.46)  yields: 


X  ~  X  +  —  e 
a  X 

a 


1  -2y[l-4(x-x  )j 


(6.47) 


Using  the  notation: 


X  ix  -  X  ) 
V  -  21  a 

^  2ir” 

e  ' 


Y  s  exp[|^  e^^] 


equation  (6.47)  becomes; 


X  Y  ^ 

X  =  Y  =^(Y) 


(b.  48) 


The  transition  occurs  at; 


X  -X  /4e 
a  ,  a  . 

Vo  "  ' 


(6.49) 


For  values  of  x  'vher'^  (6.47)  is  valid,  (6.49)  shows  that  the  transution 

temperature  is  extremely  low.  For  x  =  0.  I,  the  transition  occurs  at 

\  =3  3,  or,  in  the  case  of  at  a  temperature,  T  =  I8*^K.  Thus, 

'  o  w  3  o 

4  4'4- 

substitution  of  10%  or  more  Ti  for  V  in  lowers  the  transition 
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temperature  sufficiently  to  produce  metallic  behavior  at  a'mc=st  all 
normal  temperatures.  We  found  that  for  2%  or  less  in  the  crvaial, 

the  transition  temperature  is  changed  only  sUghtlv.  WTien  between 
2%  and  10%  is  Introduced  into  intermediate 

region,  and  the  general  expression,  {6.43),  must  be  used  without 
further  appvoximaticns.  It  can  be  shown  that  the  predicted  transition 
temperature  begins  to  drop  sharply  when  the  concentration  ex¬ 

ceeds  8%. 

Since  we  have  considered  only  the  extreme  narrow  band  limit, 
we  can  solve  only  tor  the  transition  temperature,  at  which  point  the 
two  bands  overlap.  We  cannot  find  out  if  the  large  number  of  holes 
present  in  the  valence  band  at  low  temperature  are  sufficient  to  pro¬ 
duce  a  metal  even  below  the  transition  point,  since  the  zero  width 
bands  cannot  give  metallic  behavior.  The  Fermi  energy  in  this  limit 
approaches  the  valence  band  as: 

Ep  =  E  +  kT  In  ~ 

F  V  X 

a 

for  0. 1  <  <  1,  Since  cannot  get  lower  than  the  top  of  the  valence 

band,  the  holes  can  never  become  degenerate.  But  this  is  only  be¬ 
cause  of  the  zero  bandwidth  assumption.  Since  we  have  some  idea  of 
the  effective  mass  of  carriers  in  from  the  discussion  in  sections 

A  and  B,  we  can  look  at  the  effective  mass  approximation  to  investigate 
the  degeneracy  of  the  yal^hice  band  holes. 

In  the  effective  mass  approximation,  metallic  behavior  sets 


2m,  kf  3/E 
f 


in  when: 


For  greater  acceptor  ..  oncenirmtion.  tt  r  ormi  enoigy  is  beio’*  the 

is  i 

top  of  the  vaience  band.  Using  the  value  "t  nt  /rr-  bU,  estiriia!=ar! 
in  section  A  for  met-ailic  ''^2^3' 

N  ’  ,  50  X  9  X  :0‘^  r  .3/2,,  N 

»  =  4  * - — ‘  '  '  '  17  '  ■ .  STO 

=  .  ^0  {  ^  }  ■  N  .  ib. 50 

o 

'^  +  + 

Thus,  for  Ti  concentrations  of  -0%  cr  morej  metallic  behavior 
should  occur  at  all  temperatures.  At  T  =  lOO^'K,  the  holes  in  the 
valence  band  will  be  degenerate  lor  concentrations  of  greater 

than  6%. 

MacMillan  [42j  studied  le  x^^3  system  in  detail. 

He  found  that  substiti  lion  of  Ti"^  ^  in  V,0  n  small  amounts  none- 

J 

lonicaily  depressed  the  transition  lempei dture.  This  effect  was 
measured  for  concentrations  up  to  2%,  and  i.*'  just  as  predicted  above. 
MacMillan  also  measured  conductivity  as  a  function  of  temperature 

^  '4' 

for  10%  and  50%  Ti  concentrations.  For  10%  Ti  ’  ‘ ,  the  con- 
ductivity  v/as  metallic  below  iOO  K,  decreased  slowly  about  a  factor 
of  *0^  as  the  temperature  was  raised  to  ^40^Ki  and  then  increased 
to  its  previous  value  somewhat  more  sharply.  This  *s  what  is  pre¬ 
dicted  by  (6.  50)  for  a  concentration  of  6%,  and  thus  indicates  that 
an  effective  mass  of  70  free  electron  masses  for  holes  in  ^2^3 
in  the  semiconducting  state  is  more  appropriate  than  the  value  of 
50  used  in  {6.  50).  Since  m  /m  -  50  was  the  value  estimated  for 
metallic  -^y  should  expect  a  iarg.  r  m  /m  for  the  narrower 
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bands  ot  the  semiconducting  states 

For  S0%  Ti  »^encentrationj  MacMillan  fo\uid  a  temperature - 

independent  metallie  conduetiviiy,  just  as  we  would  predict  from  the 

above  theory.  In  this  case,  the  t  -  d  band  i.i  half-filled. 

ft 

F.  Thermally  Activated  Hopping  Theories 

Som  s  authors  [64]  have  suggested  that  the  origin  of  the  con- 
ductirdty  in  the  semiconducting  state  is  thermally  activated  electron 
hopping.  Then  the  observed  activation  energy  should  be  totally  associ¬ 
ated  with  the  mobility,  and  the  concentration  of  carriers  should  be 
temperature  independent.  There  Is  now  much  experimental  evidence 
against  this  hypothesis.  The  direct  measurement  of  the  energy  gap 
in  by  Feinleib  [E4]  and  the  measurement  of  the  Hall  coefficient 

in  Ti^O,  by  Yahia  and  Frederikse  [Z2]  strongly  indicate  the  presence 
of  a  normal  semiconducting  energy  gap  with  thermal  activation  of 
free  carriers.  Furthermoie,  a  theory  which  ascribes  conduction 
to  a  diffusion  of  charge  predicts  an  increase  in  activation  energy  with 
pressure  [5,  65]  .  Such  an  increase  has  been  observed  by  Young  et  a|.[6S] 
m  CoO,  CuO,  and  Gu^O.  However,  in  the  cases  of  V^O^,  VO, 
and  Ti^Oj,  the  activation  energy  decreases  with  pressure  [21,  43,  5i]  . 

G.  Theories  which  Postulate  a  Critical  Lattice  Parameter 

The  theory  of  Mott  [4,  26,  27]  assumes  that  there  exists  a  criti¬ 
cal  lattice  constant,  R^,  above  which  there  is  insufficient  overlap  to 
obtain  physically  meaningful  bands.  Mott  suggested  that  materials 
whose  lattice  parameters  exceed  must  be  non-conducting,  and  all 
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the  electrons  have  to  be  vonsiderod  as  lor*lized  on  individual  ions 
rather  than  spread  throughout  the  crystal*  One  immediate  difiiculty 
with  theories  based  on  this  hypothesis  is  the  explanation  of  the  con¬ 
ductivity  belov/  the  transition*  If  thermaily  activated  hopping  is  con¬ 
sidered  to  be  the  conduction  mechanism  below  T  »  all  the  objoctions 
of  the  previous  section  apply*  It  is  hard  to  think  oi  any  other  model 
employing  localized  electrons  which  can  account  for  the  considerable 
(c  10  ohm-cm)  conductivity  in  the  insulating  states  of  and  Ti^O^, 
Other  objections  emerge  from  the  details  of  the  transition  in  Ti^O^ 
and  in  VO,.  At  in  Ti^O^,  the  c-axis  actually  expands.  Thus, 
the  nearest  neighbor  cation  distance  is  larger  above  than  it  is 

below  T  .  This  is  difficult  to  reconcile  with  the  idea  of  a  critical 
o 

lattice  spacing  being  exceeded  at  low  temperatures  but  not  at  high 
temperatures*  Furthermore,  below  T^  in  VO^,  the  c-axis  cations 
pair  up.  Thus,  the  nearest  neighbor  cation  distance  is  much  smaller 
in  the  non-conducting  state  than  it  is  in  the  conducting  state.  But  also, 
the  Lnear  thermal  expansion  is  such  that  at  very  h'gh  temperatures 
all  cation-cation  distances  are  greater  than  the  distance  betv  een 
the  more-separated  c-axis  cations  at  low  ^empet atures*  Thus,  even 
a  theory  based  on  an  average  also  seems  to  fa-*. 

As  a  final  point,  the  transport  properties  of  the  ^ -x)^3 

system,  discussed  in  section  E,  are  not  consistent  with  a  theory  which 
ascribes  the  transition  to  a  critical  value  of  the  lattice  parameter. 

In  particular,  TiVO^,  which  has  a  structure  w;th  lattice  constants 
very  close  to  V,0-^,  is  metallic  even  below  ISO^K,  where  both 
and  Ti^Oj  are  semiconducting. 
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H,  Theories  which  Ascribe  L-ack  ef  Conductivity  to  Direct 
Cation -Cation  Iste ’-action 

Goodenough  [28,  29,  30]  his  postulated  that  direct  interactions 
between  cations  lead  to  the  formation  of  homopolar  bonds,  and  there¬ 
fore  a  non-conducting  state  can  occur  if  all  the  3d  electrons  are  used 
up  in  such  bonds.  As  with  Mott*s  theory,  models  based  on  this 
principle  have  great  difficulty  in  explaining  the  low  temperature  con¬ 
ductivity  and  its  behavior  under  pressure.  There  is  also  some  inde¬ 
pendent  experimental  evidence  which  seems  to  contradict  this  particular 
hypothesis.  In  the  pairing  of  the  c-axis  cations  is  not  in  a 

straight  line,  but  the  pairs  are  puckered  below  T^;  yet  in  the  high 
temperature  rutile  structure*  all  these  cations  are  collinear.  If 
cation-cation  interaction  was  the  dominant  reason  for  the  pairing, 
it  is  hard  to  see  why  the  ions  would  not  move  directly  towards 

one  another.  A  more  serious  objection  is  the  lack  of  anisotropy  in 
conductivity  of  between  iSO^Kand  600°K.  Goodenough’s 

theory  predicts  that  the  conductivity  should  be  alm'^St  entirely  in 
the  basal  plane  in  this  temperature  range,  the  t^  electrons  being  tied 
up  in  homopolar  bonds  along  the  c-axis».  As  Feinleib  [21]  has  shown, 
this  is  not  in  agreement  with  experiment. 
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CORRECTiON  FOR  THE  cT  CONTFIBUTION 
TO  THE  ENERGY  GAP 


In  section  IIA,  we  wrote  for  the  energy  gap- 

g  go  ij 

where  E  is  the  gap  at  T  =  0^  /BT)  ,  and  P  s-(S  E  /Bn) 

s  g  n  g  T‘ 

In  our  fundamental  equation,  (2.  i),  we  dropped  the  term  linear  in 

T*  ^  this  appendix,  we  shaii  show  that  this  omission  does  not  ieac 

to  any  serious  errors  in  the  e-^aluation  of  conductivity  as  a  function 

of  temperature,  which  was  carried  out  in  Chapter  HI, 

We  shall  calculate  the  effects  of  the  term  -c  T  in  the  case  ol 

V,0^.  Consequently,  we  can  specialize  to  the  regions  where  the 

narrow  band  limit  is  applicable  and  Boltzmann  statistics  are  valid. 

Then,  equation  (3»  52)  in  the  Boltzmann  limit  becomes: 

-E  /2kT 

n  N  e  ^  . 

(A. 

From  equation  (A.  1): 


fA.  I 


dE 

1  _  q  dT 

dn  ^  dn 

dE 

- S  -  -  o  _  8 

dT  ‘  ^  • 

But,  equation  (A.  2)  show's: 


(A.  3) 


(A. 4) 


dn  f  £ 

dT  "  I  ^  ‘ 


,  dE 


(A.  5) 


Using  equation  (A.  4)  in  (A.  5): 
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Solving  equation  (A,  6)  for  dn/dX: 

E 

dn  n  ^ 

3T  "  7”Tir 
^  *iFr 

Substituting  (A.  7}  into  (A.  3),  we  find: 
dE 

g  ^  .  a  2kT  - 

-  »  P  “  —  g  — * —  . 

+  a 


We  may  rewrite  (A.  8)  as: 


dE 

— j  ^  ”  9  r  16  3 

dn  ^  t  pn 


2kT  -  pn  1 

^”+TT  i 

g 


Thus,  the  fractional  error  in  equation  {2.  1)  is  just: 

a  T  r  2kT  , . 

E  6  c  T  J  • 

g 


(A,  6) 


iA.l) 


(A.  8) 


(A.  9) 


For  the  cases  of  V2OJ  and  VO,  when  we  evaluated  p  from 
thermodynamics,  we  found  that  virtually  the  entire  contribution  to  p 
came  from  the  term  which  gave  the  change  in  energy  gap  due  to 
changing  volume.  We  can  assume  that  the  same  will  be  true  in 
evaluating  a  .  A  thermodynamic  argument  [66]  ,  similar  to  that  which 
led  to  equation  {2.  13),  gives: 


where  X  is  the  thermal  coefficient  of  volume  expansion,  and  k  is  the 
isothermal  compressibility.  Since  we  are  taking  the  entire  contri¬ 
bution  to  a  to  be  the  first  term,  we  have: 


a  s 


where  y  =  .  (8  E^/d  as  defined  in  equation  (2,  10), 

experimental  values  given  in  section 

IV  Ai 

K  =  40.  0  X  10*^ 

K  =  0,58  X  10'^  bar“^ 

V  1.  8  X  10  ^  eV  bar*^  , 

Then,  (A.  =0)  yields? 

c  =  -  1.2  X  10"^  eV/°K  .  ■ 

Thus,  from  (A.  9),  the  correction  to  0  near  the  transition  temperature 
is  of  the  order  of  2%  .  Thus,  the  dropping  of  the  term  in  c  T  in 
equation  {2,  1)  introduces  an  error  of  the  order  of  a  few  per  cent, 
which  is  certainly  negligible.  Note  that  smce  c  is  negative  in  V,0.^, 
equation  (A.  9)  shows  that  6  is  slightly  reduced.  This  implies  that 
the  calculated  ratio  of  to  kT^  obtained  in  Chapter  IV  should  be  a 
few  per  cent  lower.  Thus  the  transition  temperature  is  increased 
a  small  amount,  as  we  would  expect,  having  introduced  a  slight  in- 
^tease,  linear  in  temperature,  to  the  gap. 
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Appendix  B 

ENERGY  BAND  CALCULATIONS 

The  epre&ding  of  localized  ionic  energy  levels  into  energy 
bands  is  a  consequ^sncc  of  the  overlap  between  the  wave  functions 
corresponding  to  the  ions  which  make  up  the  crystal.  We  can  easily 
show  by  perturbation  theory  that  as  long  as  there  is  some  overlap, 
there  must  be  some  spreading  of  the  localized  levels,  and  thus  a 
finite  band  width.  It  would  be  most  desirable  to  be  able  to  determine 
the  energy  band  structure  of  given  materials  from  first  principles, 
but  to  date  such  calculations  are  long,  laborious,  and  grossly  in 
error.  This  has  been  particularly  true  in  the  case  of  transition 
metal  compounds,  where,  for  example,  careful  energy  band  calcula¬ 
tions  on  NiO  have  led  to  conclusions  that  one  of  the  best  known  in¬ 
sulators  is  a  good  metal  [67]  or  that  it  forms  a  crystal  consisting  of 
neutral  a*oms  [68] .  One  of  the  reasons  for  the  general  worthlessness 
of  such  calculations  seems  to  be  the  inability  to  obtain  correctly  eien 
order  of  magnitude  estimation  of  overlap.  Better  results  have  been 
acquired  for  the  oand  structure  of  crystals  when  overlap  is  treated 
as  an  experimental  parameter. 

What  we  are  really  interested  in  finding  out  is  whether,  in 
the  real  three -dimensional  case,  the  exchange  splitting  in  antiferro¬ 
magnetic  crystals  or  the  splitting  due  to  lower  symmetry  in  distorted 
crystals  is  sufficient  to  produce  true  energy  gaps  in  the  densities  of 
states.  In  an  attempt  to  ascertain  some  information  about  this,  we  . 


-shall  Outline  an  energv  band  calculation  for  VO.  We  shall  consider 
the  rock  sa’t  phasr  oi  VO,  assume  ihe  crystal  is  an  antiferromagnet , 
and  tr\  to  determine  how  strong  the  exchange  energs?^  must  be  to 
obtain  a  true  energy  gap.  The  reason  for  choosing  VO,  although 
we  do  not  beKevc  VO  to  be  antiferromagnetic,  i-s  purely  the  simplic  :y 
of  its  structure;  a  similar  calculation  for  corundum  Ti202  would  be 
much  more  difficult.  However,  we  expect  that  the  orders  of  magni¬ 
tude  involved  will  be  about  the  same  in  ail  the  oxides  of  titanium  and 
vanadium. 

We  shall  adopt  the  point  of  view  of  the  tight  binding  approxi¬ 
mation,  in  which  the  Wannier  functions  of  the  crystal  are  approxi¬ 
mated  by  the  iomc  functions  of  the  component  cations  and  anions. 
However,  we  shall  make  a  further  approximation,  consistent  with 
the  spirit  of  the  perturbation  theory  described  by  Anderson  [8]  .  VV^e 
shall  take  as  the  Wannier  functions  of  the  cations  the  wave  func¬ 
tions  plus  a  small  part  of  the  wave  functions  at  the  surrounding  oxygen 
sites.  Then  we  shaV  use  these  wave  functions  to  calculate  the 
catiOn-catxcn  overlap.  By  this  method,  we  can  neglect  the  anions 
after  it  is  dec  dr-d  what  percentage  of  O  wave  function  to  add  to  the 
V^^  function  at  each  cation  site.  As  has  been  observed  by  Amdeison  [8] 
this  method  eliminates  the  distinction  between  direct  cation-cation 
Overlap  and  indirect  cation-anion-cation  overlap. 

In  an  attempt  to  learn  something  about  the  real  problem,  we 
first  look  at  a  very  Simple  model,  a  linear  ciiain  of  alternating  V 
and  O"”  ions.  For  further  simplicity,  we  shall  deal  with  only  one 
state  from  each  ion,  Actu?Jly,  there  should  be  three  states  making 
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up  the  triplet.  We  place  ions  with  spin  down  at  x  s  (4n-l)4, 

and  ions  with  spin  up  at  x  C4n  +  l)a.  The  0“‘  ions  are  then  p,it 
at  X  ¥  Ena, 


First,  we  neglect  the  O  ions  completely.  There  are 
ions  at  x  —  a.  Let  p  s  (a  00).  Th,"'  wave  functions  are; 

^ ?  -  p)  =  y  *  F(  IF  -"p  |)c 

(B, 

C  r  +  p)  -  yz  F(  }r  +  p  })  j5 

where  F(r)  is  the  appropriate  radial  function.  Applying  Schrodinger's 
equation,  ws  obtain  the  secular  equation: 


0  ST 


e  •  E  20  cos  2ka 

20  cos  2ka  c  -  E 


where: 


(B.2) 


a  e 


/  dr  (r  H  (r) 

c  a  Hi 


“  p) 


e  -  /  dF  4,*  (  ?  +  r  +  p) 

Solving  CB.  2): 


E  = 


c  +  a 


40^  cos’ 


2ka 


CB.3) 


This  solves  the  problem  in  terms  of  the  physical  quantities  a  and  a , 
and  the  non-physical  quantity,  0.  which  is  the  effective  cation-cation 
Overlap.  S*nce  e.  and  2  are  intraionic  integrals,  they  can  be  evaluated 
from  tne  Hartree-Fock  ionic  wave  functions.  However,  we  must  find 


a  way  of  approximating  6^ 

-iOw  vhe  same  problem,  only  taking  this  time  as  the 

Wave  ,ii.iCt*ons  the  above  functions  plus  a  smai.i  part  of  the  O”” 
functions.  Then  w©  can  use  perturbation  theory  to  second  order* 
neglecting  all  out  the  physical  V+-^.0“"  overlap,  to  evaluate  both 
the  mixing  parameter  and  the  effective  cation-cation  overlap,  p. 


.  Il)  .-s  -  Bi  r  £[4,  (r)  +  4,  (?.2p)1 

O’M  o”*2 

a^p(?hp)  ?)  +  €[<!.  __  (r)  +  4-  (r  +  2p)l. 


O  ,  1 


o  *2 


The  secwlar  equation  for  this  problem  's,  to  first  order  in 


j  G  +  2e(v  +  y  i/  +  <  p  g  ») 


Z[€  /'/  h  V  COS  2ka 


CB.4i 


2|c(V  T]  :  2ka  S  -t  2ctY  ’)  +  €^(a  '+0^1 


IB,  5 


wnere; 


y  -  /  dr  4»^  (f  -  p)  Hiri  it] 

O  ,I 

/  d?  4)*  (f  -  p)  H(^:i  4  (r-2pl 

0”.2 

c*-/dr4>^  Cr)  H(T)  4>  f?) 

o“M  O’M 

a'  =  /dr4)  (r  -  2p)  H(r|  4)  {r  -  2p)  . 

O  1 2  0”,2 

Solving  equation  (B,  5),  we  obtain; 


E  =  < 


2«2 

a  +  2cy  +  cos^  2ka 

a  -  a 


fB.  61 


e  +  Zfv  +  €  . 


cos  2ka 


14Z 


where  a  »  a'  +  Comparing  (B.6)  to  (B,  3),  we  see  that 

with  the  choice: 


c  ’  +  t  ‘ 


we  obtain: 


{B.l) 


The  sign  of  0  can  usually  be  determined  by  physical  considerations. 

The  result,  (B.  7),  Is  quite  reasonable.  It  states  that  the  effective 
cation-cation  overlap  is  proportional  to  the  product  of  two  cation- 
anion-cation  overlaps,  which  is  what  we  would  physically  expect. 

Note  that  there  is  always  a  splitting  of  the  two  bands,  equal  to  |a  -  a  |, 
the  exchange  energy,  in  this  simple  model.  This  is  characteristic 
of  one-dimensicnai  models,  as  we  saw  in  Chapter  II. 

With  the  result  (B.  7),  we  can  calculate  the  important  overlap 
parameter  ^  from  the  Hartree-Fock  wave  functions  for  the  ion  [69] 
and  the  O  ion  [70].  These  Hartree-Fof  uuictions  are: 

F 

O 


r)  s  [0. 078  e  4  8.  5?.  e  +  1,  66  e  j  r 


T?  I 

(B.8i 

The  potential  V  was  obtained  by  extrapolation  from  Brown's 

o'“ 

[71]  calculations  for  Ne  and  F’,  tne  potentials  being  adjusted  by 
scaling  both  ionic  radii  to  the  ionic  radius  of  O*".  From  Pauling  [72]  ; 


f  ? .  74 


1.  iSr 


4  19. 8 e 


■2.33r 


t  37.  9e 


•4.42t 


19.3 


■  8.3lr-j  3 


,  76-  A 


1.  36 


12  A  . 


The  extrapc.at  cn  was  linear,  by  necess  ry, 
expression  for  the  potential  obtained  is; 


Thf-  best  analytic 


''q-J'-*  +  4.(i8  c-''«’''  +  0.90c-‘'''®*' 


(b.9; 


From  Watson  [69],  we  get  for  the  average  3d  electron  energy  in 
=  -2.23  Ry.  .  -30.  3  eV  . 

Watson  [70]  give«  as  the  2p  electron  energy  m  O"': 


c  a  = 


163  Ry.  “  -2.2  X  eV  . 


Callaway  [73]  gives  as  the  expression  for  the  change  in  th, 
Hartrec-Fcck  energy  when  the  spin  of  one  electron  is  reversed: 

^,-x  '  7  F^C3d,  3d)]  . 

Using  the  values  catcuiated  by  Watson  [69]  for 

^ex  ""  7  +  0.417)  Ry. 

=  0.  5  Ry.  rt  1. 9  eV  . 


(B.  iZi 


From  (B,  lO;  and  (B.  12),  we  get: 

a  =2.31  Ry.  E  31  eV  , 
£  ==  2.  '5  Ry.  ^29  eY. 

Only  y  rema  to  he  evaluated. 


(B.  13) 


(B.  14 


kns%  irom  the  usual  tight  binding  approximations  [74]  ; 

Y  ~  /  dr  f  W{  r)  I  ?  -  ) 

V  O' "  ^ 

where  W{f)  is  the  difference  between  the  one-eleetron  crystalline 
potential  and  the  potential  of  an  isolated  ion,  and  =  (OOl)a  is 
the  posi,E:iOn  of  a  nearest  neighbor  oxygen  ion.  We  shall  make  the 
following  assumptions  about  W(  F): 

a)  WCri^O  {r<a/2) 

Ui  W(r}^V(F)  Cr>a/2) 

where  V{?/  is  the  amon  potential  centered  on  the  nearest  neighbors. 
Hence; 
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i  hen  IB*  4;  becomes: 
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shall  us 


a  pjCisTfc  spUtfuidal  coordi 


naies,  with  tcsci  at  r  =  Ch  o 


Then  s.Bi  ihi  can  be  written? 
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where: 


Thus: 
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dtr 


0 


fl.r, 

d4>  cos  4>  I  di)  d|  [same  integrand] 

Jo  J.l  Ji 


r  + 


T1  = 


r  -  r' 


«*  f  >•  4 

t  =  —  i£  ri: 

a  ' 


^  ^  /a 


.2  2 


-  T|);  r  -  r***  = 

Cn  •^he  ease  c4  VO,  a  -  2.  04  A  =  4.  0  a.u.  The  integrals  are  all  eie- 
mer>ary.  be  rg  of  ’he  form  /  dx  x"  e**"*^.  After  a  long,  laborious 
calculation,  we  ob*ai.n  from  (B.  17); 

V  -  n  06  Ry. 

Thus,  from  (B.  7)  and  (B,  ^  1): 

p  =  ±0.0"  38  Ry.  =  ±0,  19  eV  .  -g\ 

Equanon  (B.  =8?  gives  the  effective  value  of  cation-cation  over¬ 
lap  in  \0.  We  must  still  perform  the  band  calculation  for  the  VO 
structure.  This  is  extremely  lengthy,  and  we  shall  just  outline  it.  In 
obtaining  (B.  iS,,  wc  le^ed  only  one  state,  i-j  =  fl5/Zn  y  z  Ffr),  of  the 
^2g  must  also  consider  the  other  two  degenerate  states: 
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^  /i|  XI  Fit) 

♦3 '  jTi  *y  ■ 


Introduce  t'-.e  notation: 

“mn*^  W(r-R.)»^(f-R.) 

^m„<*ij>  ^  I  if  Or-  R.)  W(r  -  R.)  .  R.)  . 


CB.  19) 


All  off*<iiagonai  elements  o  =  0,  The  diagonal  elements  are  equal, 
and  we  have  called  them  Cj^  ~**22  ~**33~°'  Bestrieting  the  over¬ 
lap  integrals  to  nearest  neighbors  only,  we  find: 


=  Pu' V’  '  ^23'\z'  *  P" 
Pn<«x.>  =  P22(V  '  ^33*V  "  ^ 


),,(R  )  =  0 

23  yz 


p 


1  l<«yx>  '  Pj  =  P22<^«*  “  ^33*  V  '  P33<>^x.l  ^  P’ 


Since  the  potential  W(r)  is  negative  everywhere,  we  can  conclude 
from  the  definitions,  (B.  19).  that  p  >  0,  ^‘  <  0.  We  now  introduce 
the  spin  dependence^  which  implies  a  distinction  between  the  potentials 
for  a  sptn-up  and  a  spin-down  electron,  W®  and  W^,  This  leads  to 
the  introduction  o£  the  quantities  c,  B,  p’,  for  the  ‘•excited”  band^ 
The  secular  equation  is  a  6  X  6  determinant,  which  is  quite  com¬ 
plicated  and  not  worth  reproducing  here.  Since  we  are  really  inter¬ 
ested  tn  Just  the  antiferromagnet.c  splitting,  we  set  3“  -  0,  which 
essentially  factors  the  secular  determinant  into  •‘hree  equivalent  2X2 
sub-determinants.  For  Himplicity,  we  also  make  the  assumptions 


*y 

^  ~  P*  P  “  P  j  whiCii  are  undoubtedly  quite  accurate  within  our  rarige 
of  error.  The  secular  equation  is  now  simply: 


jc  +2pcos(k^-k^)a  +  2p*[cos(k  +k 

f 

+  cos(k  +k  )a1  -  E 
j  ‘  X  y  J 


Noting  from  (B.  12)  and  (B.  18)  that  1^1, 


)a  20cos(k  +k  )a  +  2Bjeos  k  -k  a 
X  z  '  *■  V  r 

+  cos(k^-k^)a] 

Q  +  2d  cos(k  -  k  )a 

£ 

+  2d'[co8(k  -k  )arcos;s 
y  z  X 

-  -  E 

(B.20' 

Id’  !  «  !c  -  e  I  ,  the  solu¬ 


tion  of  {B.  20)  can  be  expressed; 


2pcosfk  -  k  )a  +  2d,  [cos(k  +  k  ja  +  co$(k  +  k  )al 

^  X  ^  Y  2,  ^ 

+  jC  -  c  j  -[dp  cos  (k  +  k  )a  -t  4p*^[cos{k  -  k  )a 

X  z  I  X  y 

+  cos(k^  -  ^  +  4p  P'  cos(k^  k^)a  [cos  (k^  -  k  )a 

-t  cos  :ky  -  k^la]  } 

=  a  r  Zp  cosCk  -  k  )a  +  2p'[co8{k^  +  k  )a  +  cos(k  +  k  )al 

"  X  y  y  2  * 

-  (c  -  a)  '  {same  expression  as  m  curly  brackets  above} 

(B.  2  =  ) 

with  both  the  valence  and  conduction  bands  three -fold  degenerate. 

Since  VO  has  three  3d  electrons  per  cation,  the  valence  band  in  tius 
approximation  will  be  completely  filled  and  the  conduction  band  com¬ 
pletely  empty  provided  that  E_  and  E^  do  not  overlap. 

From  (B.  Zi),  the  top  of  the  valence  band  occurs  at  the  two 
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poiiits  in  the  first  Brillouin  zone: 


4a 


(l.-i.  1) 


{B.22) 


It  ti  I  *i| 

4a  ^  • 


At  these  pc4nti 


*  3  -t  2p  4  4  f3’  I  * 


V  max 


(B.23) 

B  is  much  more  difficult  to  find  the  bottom  of  the  conduction 

hand,  since  there  is  no  point  in  k-space  where  the  terms  in  E  are 

c 

Individually  minimized. 

Up  until  this  point,  we  have  been  completely  general.  It  is 
now  necessary  to  put  in  values  for  3',  o,  2  in  order  to  proceed. 
However,  note  that  the  bottom  of  the  conduction  band  must  be  >igher 
than: 

(E  )  .  >  ,  .2P.4IPM  - 

C  min  ^  'f'  ‘ 


c  -  c 


Thus,  there  is  always  an  energy  gap  if; 

(c  -  a  )  >  43  4  8}3M  +  . 


(B.24) 


e  -  c 


In  order  to  evaluate  3',  we  use  the  relationship  found  by 
Fletcher  [75]  fcr  the  ratio  of  the  same  overlap  integrals  in  the  case 
of  the  3d  electrons  of  metallic  nickel: 


.0.297 


Equation  (B.  18)  then  yields: 


3»  ”  -0,041  Ry  =  -0.06  eV 


(B.25) 


oi!i:!n<|i'!!!!i'iiiiituiiiiiiii'iiiiiiiiniiiiiiuiiiiuii>t  ..in'iiiin 


8 


e  can  now  use  iB.  iZi,  (B.  !8),  and  to  test  eqaation  (B.24i, 

We  find; 

e  ~  Q  =  1.  9  ev 


4p  +  8  i0M  +  - 


4  eV  . 


Thus  equation  (3.24)  Is  satisfied,  and  the  true  energy  gap  is  at 
least  0.  S  cV. 

We  can  also  easily  obtain  an  upper  limit  for  the  band  width. 
From  {B.2i;-,  it  is  clear  that: 

<  4p  +  8  |0* !  +  4  2  ip*  i)^ 

o  -  c 

^  4  eV  . 

In  order  to  find  the  real  values  cf  and  E^,  we  must  locate 
the  conduction  band  minimum.  Since  p  >  !0*|,  the  minimum  occurs 
near  the  pianes,  =  k^±  w/a.  Minimizmg  the  other  terms  subjcc" 
to  this  constiamt,  it  can  be  seen  that  the  bottom  of  the  conduction 
band  is  ’n  the  v^cmity  of  the  twelve  points  equivalent  to: 

1,-3). 

The  conduction  band  minimum  is  approximately: 

-  -I  !»' I  +  (B.26) 

(3,  Zbi  to  (B,  23l,  we  conclude: 

E  =  (s  ,  £  )  -  4p  -  4  jp»  I  + 


=  0.9  eV  . 


(B.Z71 


ISO 


It  is  ilear  from  (B.Zl); 


s  4P  +  4p* 


40*)^ 


=  1.0  eV  . 


{B.E8) 


it  iS  also  useful  to  calculate  the  effective  ma  ses.  The  effective 
mass  of  holes  in  the  valence  band  is  easily  determined  from  (B,  El), 
it  iS  found  that  the  band  is  extremely  anisotropic  near  the  maximutnj 
the  diagonalized  effective  mass  tensor  can  be  expressed,  in  units 
of  the  free  electron  mass: 


m,  =4;  ^  10;  m*  =  90.  {B.E9) 

Defining  a  "density-of-states "  effective  mass  as; 

m  =  M{m*  m^  m* 

(B.  29)  shows: 


iB.30) 

The  final  results,  {B,27)  and  (B.  28),  are  each  about  a  factor 
of  5  greater  than  we  expected  from  the  experimental  observations 
ill  Chapters  IV  and  V,  This  is  undoubtedly  due  to  proportional  over¬ 
estimates  of  the  import<ait  band  parameters  (a  -  o  ),  0 ,  and  0*.  Such 
overestimates  are  a  general  characteristic  of  band  calculations.  The 
explanation  may  be  in  covalency  effects,  most  of  which  are  neg¬ 
lected  in  this  type  of  calculation,  or  it  may  be  due  to  a  serious  deriva¬ 
tion  of  the  ionic  wave  functions  in  the  crysi.  i  from  those  of  the  iso¬ 
lated  ions.  Note  that  the  effective  mass  value,  given  by  (B.  30)  is 


quite  reasonable  for  VO. 


Appendix  C 


HYDROGEN  MOLECULE  AND  MOLECULAR  ION 


WITH  DELTA- FUNCTION  INTERACTIONS 

In  section  ii  D,  we  solved  the  problem  of  a  one -dimensional 
crystal  with  two  cations  per  unit  cell,  the  cations  attracting  the  3d 
electrons  by  means  c-f  a  delta -function  interaction.  We  demonstrated 
that,  when  there  is  one  3d  electron  per  cation,  such  a  crystal  will 
be  semiconductmg  due  to  the  doubling  of  the  periodicity  arising  from 
the  pairing  of  cations.  As  we  showed  in  section  U  D,  excitation  of 
an  electron  across  this  energy  gap  tends  to  decrease  the  gap.  How¬ 
ever,  It  is  not  Ofawcus  that  such  excitation  also  lends  to  increase 
the  average  distance  be*ween  the  clcsey-spaced  cations. 


m  ordet  demonstrate  the  decrease  of  crystalLne  distortion 
as  electrons  move  from  bonding  to  antibonding  bands,  we  used  an 
analogy  relatmg  to  the  hydrogen  molecule  and  the  hydrogen  moiecuiar 
iOn.  With  both  electrons  tn  an  H^  motecule  in  bonding  orbitals,  the 
equilibrium  distance  between  the  two  ions  is  0,  74  A  [76]  .  If  one 
of  the  two  electrons  placed  in  an  antUbondmg  orbital,  the  molecule 
iS  no  longer  bound,  and  the  equilibrium  distance  becomes  infinite, 
in  the  band  situaLcn,  it  is  more  reasonable  to  look  at  one  electron 


being  taken  out  of  a  bonding  orbital  and  put  into  an  antibonding  orbital 
m  another  unit  cell.  Thus,  we  should  compare  the  H^  equilibrium 
distance  with  that  cf  an  molecular  ion,  whose  H'*'  tons  are  L  06  A 
aprtxt  [f7j  ,  Theri-fcre,  the  ions  move  apart  a  factor  of  1,43  when 
an  electron  is  removed  from  a  bonding  orbital. 
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Howevgi",  tile  first  cftSc  out  in  section  II  I)  is  for 

electrons  attracted  to  the  cations  by  delta -functioft  interactions,  which 
are  essentially  iminiteiy  screened  Csulotnb  interactions,  and  are 
very  different  from  pure  Coniomb  potentials  as  can  be  seen  from  the 
work  of  Kohn  |3i|  discussed  in  Chapter  L  Therefore,  it  is  imoortant 
to  see  if  the  hydrogen  tr olecaie  analogy  holds  up  if  the  electrons  are 
attracted  to  the  H  ions  by  means  of  a  delta -function  potential.  Such 
a  calculat'-on  ts  the  purpose  of  this  appendix. 

Consider  an  '’H2”  molecular  ion,  with  the  following  Hamiltonian: 

2 

(C.  1) 


H 


«1  .  V  6M  -  V  6{x  .  X)  4  ^ 

dx 


where  X  >  0  is  the  distance  between  the  ions.  Take  the  zero  of 


energy  at  e  /X,  Then  Schrodmger’s  equation  can  be  written: 

H  i(x)  =  41  (x)  . 

The  true  energy  is; 


(C.2; 


Let: 


E, 


k  s 


(C.5) 


For  X  -  (-00,  0),  equations  (C.  1)  and  (0.2)  give: 


(x)  -  k  4>(x)  K  0  . 

Solving  (C.4),  and  applying  the  boundary  condition  at  x  =  -co: 

-  .  .  kx 

ip{}i)  ~  A  e 

Similarly,  for  x  =  (0,  X),  the  solution  of  equation  (C.2)  is: 


(C.4) 


(C.5) 
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^  B  +  C  . 

Finally,  in  the  region  x  -  (X. «),  the  solution  of  (C.E),  using  the 
boundary  condition  at  x  =  ce,  is: 

4»(x)  =  D  e'^  . 

From  (C.  5)  and  (C.6),  continuity  at  x  =  0  implies? 

A  =  B  +  C  . 


Similarly,  from  (C.  6)  and  (C.  7),  continuity  at  x  =  X  yields: 


n  2kX  „  ^ 

Be  ~  D  +  C  . 


Integrating  (C.4)  from  x  =  0“  to  x  =  0^* 


<|>*(0")  -  =  --q  <|>(0} 


where: 


r\  s 


EmV 


ti 


Z 


Similarly,  integration  of  (C.4)  from  x  =  X“  to  x  =  implies: 
t’(x')  -  o^CX”)  =  -  ti(j>{X}  . 

Using  {C.5),  (C.  6),  and  (C.  7),  equations  (C.  10)  and  (C.  11)  yield? 

k 

D  -  C  +  B  ^  D  . 

k 

Solving  (C.  8).  (C,9),  (C.  12).  and  (C.  13): 

B  =  i  e-2>“D 
2k 


B 


(2k  -  n) 

ZkT) 


.2 


D  * 


(C.  6) 

(C.7) 

(C.8) 

{C.9) 

(C.  *0) 


(C.  il) 

(C.  12) 
(C.  13) 


(C.  14) 
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nations  iC*  14)  kave  a  non-triviai  solution  only  If: 


1  .2kK  .2 

o  s  t2k  -  ti)  , 


Let  ^  -  Zk  1|,  Then  (C.  15)  can  be  written: 


|C.  15) 


-  2 


Let  e  =  e 


1 

^t|X 


Then  {C.  16)  becomes: 


{C.  16) 


<C.  17) 


1  1_ 

iJL  = 

‘  e 

Let  Ms  ^/t|c  ;  N  =  c  ,  Equation  (C.  17)  is  now  expressed: 

iM  =  N^.  (C.  18) 

There  are  thus  two  solutions  to  Schr6dinger*s  equation.  Taking  the 
positive  sign  of  tC.  18)*  we  find: 


M  =  =  N  ^ 


N  • 


n(N)  . 


The  solution  (C,  19)  corresponds  to: 


k  =|[1  -t  3  '  q{N)] 


CC.  19) 


-  K[  I  +  ap(e°)]^  . 


(C.ZO) 


2  2 

where  K  s  m  V  /2h  .  The  negative  sign  in  (C,  18)  yields: 
C-M)  -  , 


The  solution  of  CC.21)  is  clearly: 
M  '  -  1 


CC.21) 


which  corresponds  to= 


K  ^  u 


Hence,  this  solution  does  not  produce  a  bound  state.  There  is  one 
and  only  one  bound  state,  given  by  (C.  20).  From  (C.  1),  ihe  total 
energy  of  this  s*ate  is; 


E  -  -Kf  i  +  c  T|{c“)] 


®  n  2  .  e^ 


X  ■ 


?r  2  ? 


The  equilibrium  distance,  X_^,  is  relatively  close  lO  X  '  0.  Note  tha 
m  the  vicinity  of  X  =  0,  (C.  20)  implies: 


«  -  4K  I  ^  qX). 

New,  cens  der  the  m-  ecule,  with  the  Hamiltonian: 


-2  „  2 

H  = 

G  2m  ^  2 

<fx. 


V  6{x,)  .  V  6(x,  -X) 


^C, 23, 


Vdfx^)  -  V6{x2  -X) 


V6(x,  ‘X^) 


Schrodinger 's  equation  is: 


^o^Cxp  x^)  =  EOilxj.x^i 


(C.24) 


(C.251 


Since  the  ’  mcleculai  ion  had  only  one  bound  state,  we  know  that 
the  solution  of  (C.  24|  and  \C,Z5]  neglecting  the  electron-electron 


interaction,  is; 


=  4>(x^)  4»{x2)  -  c^Pi  ]  I 

From  (C.  5),  C.  6;,  and  (C.  ■"),  using  equations  (C.  8),  (C.  9),  (C.  "2;, 


IC.26) 


Illl:  .ml  l(ll|l|l|l|l  il 


and  CC»  13),  we  find: 


kx 

n  ® 


{— ae,  0) 


X  J  (2k  -  «)  kx  .  Ek  -  n 

"  1  SloT"  ^ 


(O.X) 


(C.2?) 


(X.  oo)  . 


Taking  into  account  the  electron -electron  interaction,  the  solution  of 
(C,  24)  and  {C.25)  is: 


where: 


E  s  2E^  +  ^  +  E’ 


V  /  dx  i  i>(x){  ^ 


[  /  dx  !#(x)r  ] 


(C.28) 


(C,29) 


Near  X  =  0,  (C.27)  and  {C.29)  give: 

E's  2K(1  -  I  nX)  . 

For  iarg»*  X,  {0.27)  and  (C.29)  show: 


(C.30) 


E’  ZK  , 


(C.3I) 


Analogous  to  the  case  of  the  real  hydrogen  molecule,  the 
electron-electron  interaction  is  relatively  independent  of  the  distance 
between  the  H  ions,  and  hence  doesn't  affect  the  calculation  of 
equilibrium  distance  in  a  major  way.  Using  {C.30)  and  {C.3i),  (C.2S) 
yields' 


tX  0) 


IB 


-  6K{1 


(X  -  oc)  . 


(C.  32} 


Comparing  (0,23)  with  (C.  32),  we  see  that  for  the  ‘‘H,”  molecule, 
the  equilibrium  diatance  is: 


X 


oo 


=  0.  ?6  X  . 
o 


{C.  33i 


Equation  {C.  33)  ahowg  that  removal  of  an  electron  from  a  bonding 
orbital  in  increases  the  separation  of  the  ions  by  a  factor  of 

1.32,  very  close  to  the  value  in  the  real  case.  Hence,  the  hydrogen 
molecule  analogy  used  in  section  II  D  remains  valid  for  delta- 
function  interactions. 
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